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REVISED DECLARATION BY AYMERIC PERCHANT UNDER 37 C.F.R. S 1.131 

In connection with the Applicant's Response to the Final Office Action issued on April 
7, 2009, this declaration sets forth the pertinent facts proving conception of the claimed 
invention prior to December 2001 . Further, this declaration also sets forth the pertinent facts 
proving due diligence of the inventors from conception of the claimed invention until the filing 
date, July 18, 2002, of the priority application related to the present application (constructive 
reduction to practice). 

I, Aymeric Perchant, hereby declare that: 

1 . I am a co-inventor of the subject matter described and claimed in the above-identified 
application, which relates to a method for processing an image acquired by means of 
a guide consisting of a plurality of optical fibres. 
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2. As evidenced by the attached copy of laboratory notebook, pages and the 
corresponding English translations of the relevant notebook pages, we performed 
various testing of the claimed subject matter as early as July 2001, which is prior to 
December 6, 2001 , the priority date of Shankar el al. (U.S. 6,885,801 ). 

3. We had been diligent in working on this invention from prior to December 6, 2001 to 
Ihe date of reducing it to practice. Evidence of diligence is attached [in French]. To 
accompany the diligence documents, the following spreadsheet contains an English 
description of each diligence document and how each diligence document relates to 
the image processing technology of the present application. The spreadsheet below 
also includes the dale of each diligence document. The dates range from September 
4, 2001 to November 4, 2002. 



4. 1 hereby attest that the below English descriptions and dates corresponding to each of 
the diligence evidence documents are accurate and correct. 



Filename 


Date 


Title 


Authors 


Description/link with image 
processing 


FMODIF_ittterface_ 
eellViZio.ps 


18.Mar.20O2 


Functional 
Requirement of 
Software Interface 
for Software 
prototype 


Georges l.c 
Goualher, 
Aymeric 
Perchant 


Requirements of the first software 
prototype controlling the tomoscope, 
former name of Cellivizo prototype 
device, and including image 
processing capabilities. 


MEMO acquisition 
NOV2001.ps 


27.Nov.200! 


linage Acquisition 
in November 2001 


Georges Lc 
Goualher 


Image Acquisition protocol with the 
tomoscope prototype to acquire proper 
raw image for image processing 


MEMO_devlog.ps 


I.Fcb.2002 


Software 


( imialher 


So It ware architecture, including image 
processing. 


ManuelUiilisation/ 


5.Aug.2002 


CellvizioUser 


Georges 1 e 
Goualher. 
Aymeric 
Perchant 


Software guide as a web page. 
Screenshots show the fiberdetect 
functionality, the background 
subtraction. 


ReferencelmageCell- 
2002-02-07.ps 


7.Feb.2002 


Reference Guide 
oflmageCell 


Aymeric 
Perchant 


ImageCell if the name of the image 
processing included in the Cellvizio 
{ or tomoscope at the time) software. 
Drawings show the implementation of 
the filler removal algorithm, and 
image reconstruction. Results from the 
Nov. 200 1 database and algorithm 
block-diagram are shown. First 
release. 
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ReferencelmageCcII- 


12.Mw.2002 


Second release of 
.Reference Guide ... 


Aymeric 
Jferchanl 


Second Release. 


2ffl)2.-03-12,K 

RcfcrcnceimageCell- 


24.Apr.2002 


Third release of 
Reference G.ukk — 


Aymeric 


Third Release. 


Cahicrdcseharges 

PROio.ps 


13.Sep.2001 


requirements of 
the software 
prototype of the 
tomoscope 


Georges Le 
Goualher, 
Aymcric 
Perchant 


Includes image processing module and 
technical choices are addressed. 


im igeCcUSA_cdc.ps 


4.Nov.2002 


Stand-alone 
interface for 
ImageCel! 


Aymeric 
Perchant 


Command line specification for raw 
image post processing. 


traitementimage.ps 


4.Sep.2001 


Image Processing 


Aymcric 
Perchant 


Optima! estimation of photon signal in 
the fiber on a raw image. Noise 
analysis. 



We filed an application in France on July 18, 2002, which is the priority application 
for PCT/FR03/02197. The present application is a national stage application based 
on this PCT application. The diligence evidenced attached shows that we had been 
diligent in pursuing the present invention from prior to Dec. 6, 200 1 (the priority date 
of Shankar) to the date of filing of the PCT application, which is July 1 8, 2002. 

All statements made of my own knowledge are true and that all statements made on 
information and belief are believed to be true, and further that these statements were 
made with the knowledge that willful false statements and the like so made are 
punishable by fine or imprisonment, or both, under Section 1001 of Title 18 of the 
United States Code and thai such willful false statements may jeopardize the validity 
of the application or any patent issued thereon. 



Respectfully SubmiUed, 



Date:, ^! /b ; 2 &>3 




Attachment. 



3 



I§1 

0)0 w 

2 « g. 
„, 3 w 
2 ST M 



o I 8 

c = 

■% "§ 

o S cd 

c . » 

o co £ 

■° "co - 0 

o) E 2 



w cu 5, 

TO- | £) | 

| | § | 

CD E CO £ 

g 2 JH 2 



E <S 



* i 

CD -2 
O) C 
O) — 

S E 



o o ^ 0 

? w co 
£ E 5 £ 




J 



x: cn 



C © 

o o 
cn Q. 



CO 



£2 i 

CD 

1 = 

5 -a 

Si 




1? 



E <d 
2 =5 

Is 



E a) 

CD £ 



if 



25 f 

; CD 5 



= 5 ; 

A 



8> ^ -£ 



.E cd 

'5 I 8 

CD x 

"9 =i= CD 

=5 o £ 

o)*- *: 

CD 03 2 



Ego*; 

5 S s * 
- ® 8 2 

>. JQ 2 o 

© m -° O 

i c/> o x: 



CD | E 



E a) 
5 o) 



E © 



"3 C3) Q. s-I =^ o ! 



n E c 

.9 £ CO 

pi 

CO o 8 



0 



O O) 

9- 03 



8 s 




us 



JI0LU8LU eqi in epnpui qi 



0 



o. w o ; 

1 1 1 L ^ 

i; gait 

■8 « < 



S in 



§ 52 
£ ° 



o "ffl 



o v c 



E | 
o o 



X5 y O 

US 



c O i_ 

E o ° 



E 



5 E 



<D O T3 N 3 

,P 5. w w w 



5 
1 



E "§ ' 

CO o 



SZ CD 
.2 ^ 



> => Si > 



Si 

o E 



^ o s - 

0 *" £ - 

~ CD 
O) o 



© © 8 © 

£ c © £ 

1 1 "3 o 

*- ~ W -n 



* I- 



4 -+S 




\ 



as 



o r; 
£ o 



O (1) 

If 



11 



(1) 3 



il 



g cn 

I f 



9- ° 



3 ° C 

o c 

«5 « .S3 

l> :> o 




sz c 



ii o 

CD O 

E £ 

CO ° 

£ -£= 

t! ^ 

CD "(D 

5 C 



0) 5 



CD .£ 

£ E 

© CD 



7 



> sz 
o cp 
E 'CD 



f -1 





the 


: iber 




N 
CO 


en 


'c 




O) 

c 




'E 




suiti 


N 


CO 






w 


■o 












& 


.9 "a 

Vv 'to 








-a 




c 
o 




a 


o 
cd 






her 


.c 


be 




qet 


00 


am 


too 


o 


o. 




CD 


CD 
CO 


T3 






ZS 


O 
-C 


1 


CD 
.Q 


-fiberj 


(0 
CD 


C 


O 




USIC 


CD 


'c 


CD 




ha 


E 


















o 


CO 






E 


CD 


an 



4- 0 

E ■- -o 

o C 0 

■*= o t3 

m "O CD 

§ 









CD 


c 




CO 


> 






o 




ms 


em 


c 

CD 


_u 






unct 


-o 




an 


rp 


CD 


o 
c 


ibe 










o 


CD 


ST of 


r dete 




_i 


o 





E 



1 § 
is 



2 .2 c 
-P 

3 ffl O 

3 H5 cd 
5 >> ^ 



£(1)3 
<D CD 



■- E 5 
CD W _ 



~ 0) ? 3 

b £ 2 o 

CD --9 = ^ 

= 8 



1 ial-i 



o 




If 

CD 'a} 
* c 

1 i£ 

Hi 



81 



° 8 



:> 

11 



Diligence Document 1 



Modification de l'interface graphique du 
module d'acquisition cell-vlzlo 

Auteur du memo : Florence Meusburger 
Date : 11 03 2002 

Diffusion : interne 

Section : developpement informatique 

Version : Revision : 1.5 

Reference du document : $Id: FMQDIF_interface_cellViZio.tex,v 1.5 2002/03/18 12:51:19 florencem Exp $ 

1 Introduction 

Ce document est le cahier des charges des modifications a apporter a l'interface graphique actuelle 
du prototype Cell-viZio, dans le cadre de la mise en conformite du prototype avec la reglementation 
applicable (APAVE). 

2 Modifications : 

Lancemcnt et sauvcgarde de l'acquisition d'images : 

Seuls deux boutons permettent de commander l'acquisition des images : le bouton Start/Stop et le 
bouton Save. On supprime le bouton Pause/Continue qui n'a pas grand interet pour 1 'utilisateur 
standard. 

Langue : II faudrait que l'interface soit. disponible dans di verses langues : Putilisateur choisira.it 
celle souhaitee. Pour l'instant, Panglais est impose. 

Parametres electroniques : 

Les parametres electroniques une fois regies n'ont plus besoin d'etre modifies. C'est pourquoi, la 
zone de saisie des parametres electroniques disparait de la fenetre principale du module d'a.cqui- 
sition. Mais, il faut pouvoir modifier ces parametres en cas de dereglage. L'entree "Parameters" 
du menu "Calibration" entraine 1'ouverture d'une fenetre de modification des parametres electro- 
niques. L'aeces a cette fenetre doit etre limite aux personnes autorisees. Pour cela, deux possibilites : 
- On peut mettre en place un systeme de login au lancement de l'application permettant de dif- 
ferencier les utilisateurs. On distingue deux "types" d'utilisateurs de l'application : le "simple 
utilisateur" et "radministrateur". Seuls les "administrateurs" ont acces a la fenetre des para- 
metres. Ceci implique qu'un utilisateur se charge de creer un "compte utilisateur" pour chaque 
utilisateur, et que chaque utilisateur se logue pour pouvoir lancer l'application, ce qui est tout de 
meme contraignant. Cependant, cela. permettra.it de garder la trace dans un fichier des differentes 
actions des utilisateurs, mais est-ce vraiment. utile? 
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- Une autre solution consiste a demander un mot de passe a l'ouverture de la fenetre des para- 
metres. Seuls les detenteurs du mot de passe ont acces a la fenetre des parametres. Ce mot de 
passe est entre une seule fois, a la premiere ouverture de la. fenetre. 



La fenetre des parametres ressemble fortement a la zone des parametres actuelle. Le bouton send 
est supprime (pour eviter d'oublier...) : tout changement enframe la modification immediate de 
l'electronique. L'utilisateur a la possibility d'enregistrer une configuration de parametres avec le 
bouton save, configuration qui sera automatiquement chargee au lancement de 1'application. Suite 
a des modifications des parametres, il peut recharger la configuration sauvegardee avec le bouton 
load. Cette fenetre peut rester ouverte pendant, l'acquisition. 

Certains parametres electroniques sont cependant directement accessibles dans la fenetre princi- 
pale du module d'acquisition et modifiables par tout utilisateur : la resolution de l'image (640x640, 
512x512, ou 384x384), et le "horizontal time shift". Leur modification est prise en compte instan- 
tanement. 

Traitements d'images : 

La fenetre de traitement d'images est conservee : l'utilisateur a le choix du nombre d'images a 
"bufferizer", il peut lancer l'acquisition de fond avec le bouton Background, soustraire le fond avec 
le bouton Subtraction, et lancer le traitement d'une image avec ImageCell. En revanche, revolution 
des boutons FiberDetect et TempMea.n va dependre de la faisabilite de traitements en parallele de 
l'acquisition. 

II faudrait mettre en place un systeme qui, au cours de l'acquisition, determine si la detection des 
fibres estimee auparavant est correcte, ou si elle doit etre refaite. Ainsi, on detecte les fibres au 
lancement de rapplication, puis quand le systeme le juge necessaire, suite a l'accord de l'utilisateur. 
Le bouton FiberDetect sera alors supprime. 

Pour la moyenne temporalis, l'ideal serai t. de la calculer regulierement. et de la soustraire aux images 
acquises ; il reste a determiner si on a le temps de le faire en temps reel... Si tel est le cas, le bouton 
TempMean deviendra un "toggle button" : la moyenne sera soustraite si et seulement si le bouton 
est enfonce. 

Ainsi Pinterface ne contiendra plus que quatre boutons de traitement : Background, Subtraction, 
TempMean et ImageCell. 



Securite laser : 

Lorsqu'on demarre le soft, le laser est automatiquement eteint. L'allumage du laser s'effectue deux 
secondes apres le lancement de l'acquisition avec start, pendant lesquelles un message clignotant 
"laser on" apparaTt a l'ecran. Pendant ces deux secondes, la procedure d'activation de 1'emission la- 
ser peut etre interrompue par un stop. Puis un temoin bleu bien visible a l'ecran indique que le laser 
est allume, pendant toute la duree de remission. L'appui sur le bouton stop ent.ra.5ne 1 'arret du laser. 

Sauvegarde du fonctionnement du logiciel : 

II sera.it interessant d'enregistrer dans un fichier quand et pendant combien de temps le logiciel a 
ete utilise, ainsi que les exceptions generees au cours de 1 'utilisation du logiciel ; ceci afin d'estimer 
le bon fonctionnement du logiciel et de deceler les bugs. 

Ameliorations a apporter au viewer : 

Voici une liste non exhaustive des ameliorations envisageables pour le viewer : 

- trouver un moyen de diminuer le scintillement engendre par la visualisation d'un film. 

- dans la ViewPalette, les seuils de la palette de couleur ne sont pas manipulates de maniere 
intuitive. Pour l'instant, on fixe les seuils a l'aide de clics droit ou gauche de souris. II serait 
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peut-etre plus facile d'ajouter des "barres" sur la ViewPalette, materialisant les seuils. 

- On n'analyse pas 1'image ou le film affiche afin de trouver sa palette de couleur. Quand on 
affiche la palette associee a 1'image ou au film ouvert, on suppose toujours que c'est une palette 
en niveaux de gris... 

- permettre a l'utilisateur de choisir la zone a zoomer a l'aide de la souris. Pour Pinstant, la zone 
zoomee est la region en haut, a gauche de 1'image. On peut par la suite se deplacer dans 1'image 
zoomee a l'aide de la souris ou des fleches. 
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Diligence Document 2 



Acquisitions d'Images Novembre 2001 



Auteur du memo : Georges LE GOUALHER 
Date : 27 Novembre 2001 

Diffusion : interne, NA 
Section : optique-image 

Sujet : constitution d'une base de donnees images 
Version : Revision : 1.3 

Reference du document : $Id: MEM0_acquisitionN0V2001.tex, v 1.3 2001/11/27 16:22:38 georges Exp $ 

1 Objectif 

Faire une serie d'acquisitions afin de constituer une base de donnees d'images signifkatives mon- 
trant le savoir faire de MKT en acquisition d'images de retrodiffusion en novembre 2001. 

Ces images feront l'objet d'un memo qui presentera les images acquises en mode non-fibre, en mode 
fibre et des traitements d'images qui leur seront appliques pour montrer l'apport du traitement 
de l'image. On veut done acquerir des images comparables en mode fibre et non-fibre (amplitude 
spatiale, grandissements similaires). 



2 Protocole 

On desire faire des images en mode fibre et non-fibre dans les meilleures conditions possibles. 
Pour le mode fibre, on utilise done la tete depliee et differents grandissements (i.e. differents ob- 
jectifs : 

- x25 (Gtete=l ; Gsysteme=6) ; 

- x40 (Gtete=1.6; Gsysteme=9.6) ; 

- xlOO (Gtete=4; Gsysteme=24) ; 

Pour le mode non fibre, on utilise un objectif : 

- x25 (Gsysteme=6) ; 

- x40 (Gsysteme=10) ; 

- xlOO (Gsysteme=25) ; 

Pour chacun des elements, on enregistre des images a differentes profondeurs. Le pas en Z peut 
etre adapte en fonction de la resolution axiale au grandissement. courant. 
Les elements a imager sont : 

- la mire USAF ; 

- la mire de distorsion ; 
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- le circuit integre ; 

- l'oignon ; 

- muscle de poulet. ; 

- des tissus (a definir) ; 

Points importants lors de 1 'acquisition : 

- Pour chacunes des acquisitions il faut connaltre l'amplitude du champ balaye (afin de determiner 
la resolution de l'image acquise) ; 

- II faut etre capable de connaitre la resolution laterale et. axiale des images ; 

- La profondeur doit etre connue ; 

- Une acquisition du fond (reflexion parasite) est enregistree (20 images) ; 

- Les images ne doivent pas etre trap saturees ; 



3 Organisation des donnees 

/NovOl /fibre/US AFGsysl, USAFGsysl_6,..., Distorsion, Oignon, Tissus(a preciser) 
/NovOl/nonfibre/USAFGsysl, USAFGsysl_6,..., Distorsion, Oignon, Tissus(a preciser) 
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Diligence Document 3 



Developpements Logiciels 



Auteur du memo : GLG 
Date : 02 01 2002 

Diffusion : 
Section : 
Sujet : 

Version : Revision : 1.2 

Reference du document : $Id: HEMD_devlog.tex, v 1.2 2002/01/02 18:15:01 georges Exp $ 

1 Introduction 

Le but de ce document est. de definir les developpements logiciels a effectuer sur 1'application 
principale du tomoscope mktApp. Le texte en italique represents des elements non prioritaires. 
Les references aux Ml, M2 etc... font reference au CR. de Reunion Image avec le E. Borotto. 

2 Application Tomoscope : mktApp 

2.1 Creation du Module de Visualisation : ModVisu 

Le module de visualisation devra.it er.re un composant Qt. utilisable en tant que tel, utilisee dans le 
module d'acquisition modAcq (zone d'affichage) et dans le module de post-traitement modPro- 
cessing. 

Le module de visualisation permet. de visualiser une image fixe (8 ou 10 bits ou plus) de type 
BMP, PNG ou un film INR. ou MNG (8 ou 10 bits ou plus) sous forme de sequence animee. 

Operations pour la visualisation (') : 

- le module de visualisation permet de visualiser une image fixe (8 ou 10 bits ou plus) de type 
BMP, PNG ou un film INR. ou MNG (8 ou 10 bits ou plus) sous forme de sequence animee. 

- effectuer un zoom (avec les boutons i et -) 

- translator 1'image (avec les fleches) 

- appliquer une LUT (niveau de gris ou fausses couleurs). L'ensemble des LUTs pourra etre evolutif 
(ajout facile de nouvelles LUTs). L'utilisateur pent agir sur les niveaux min et. max pour faire 
un fenetrage (lineaire ou logarithmique) . On pent egalement inverser les niveaux de gris. 

- appliquer une transformation gamma sur le fenetrage courant. 

- afficher l'information associee a 1'image (M5) (commentaire contenu dans le fichier cf g si ce- 
lui ci est disponible). L'utilisateur doit pouvoir modifier le commentaire associee (est-ce que 
tout les utilisateurs peuvent le faire?). Cela permet a l'utilisateur de completer et/ou rectifier 
1 'information associee a une image. 

'cf. display.exe 
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- Faire un outil interactif permettant de tracer un trait, entre deux points et d'afficher la longueur 
de ce trait. 

- possibility de selectionner une zone pour y faire une operation (mesure de surface, histogrammed 

- afficher des regies sur les cotes (indices des pixels ou coordonne.es micrometriqu.es) 

- afficher les axes indignant la position du curseur 

- inversion des niveaux de gris 



Operations pour V export : 

- exporter (enregistrer) I'image courante (avec application de traitements eventuels) dans un for- 
mat adequat avec ma.rqu.age MKT optionnel en tant que tel. ou dans un traitement de texte en 
vv.e de I 'edition d'un rapport. 

3 Modifications pour le module d'acquisition : modAcq 

3.1 Modifications identifiers 

- Cf. Sourceforge 

- RealTomoscopePara : singleton 

- Ml 

- M3 

- M10 

3.2 Modifications a approfondir 

- Mil? (Remarque : influence sur le module de viusalisation) 

- Integration du module de visualisation : La zone d'affichage doit etre remplace par le 
module de visualisation. 

- Optimisation de la procedure d'acquisition et du traitement associe : (methode doItQ 
de AcquisitionThread) : il faudrait optimiser les acces aux donnees et calculs effectues. Un 
flowchart (parcours des donnees, calculs) precis devrait etre represents afin d 'optimise la boucle 
d'acquisition. Pour effectuer cette modification, il faudra prendre en compte les besoins en TDI 
(cf. section 6). 

- Choix de l'interface d'acquisition : On doit pouvoir passer rapidemment d'un type d'ac- 
quisition vers un autre : mode sirnule, mode AMCC, mode NI1422, mode FireWire. Cela peut 
etre fait a partir d'une classe generique (interface) et de differentes implementations derivee. 
La classe generique est a definir (elle peut par exemple contenir les methodes open(), close(), 
reset(), getData(),...) -voir dans le cas de la carte NI et du FireWire si cela est compatible-. 
On n'utilisera que les fonctions minimales de la carte d'acquisition (getDataQ), les traitements 
(application de LUTs, inversion des lignes) seraient fait en soft. Le point dur peut etre de faire 
en sorte que le passage 8 vers 10 bits ou plus soit quasi transparent. 

- Fonctionnement en mode 10 bits : Ce probleme est lie a la tache precedents. En bref, il 
s'agit d'etudier les points suivants : Acquisition des donnees en mode 10 bits, Prise en main de 
la carte NI; Format d'enregistrements des donnees: Affichage des donnees (choix d'une LUT, 
d'un gamma) ; Distinction des donnees brutes avec les donnees affichees et/ou traitees. 

4 Modifications pour le module browser : modBrowser 

Les modifications de ce module sont liees a la tache de definition d'une base de donnees d'images 
(cf. ...) . 

Parmis les points deja evoques : 
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- Sauvegarde des donnees : il faut un utilitaire permettant de faire une sauvegarde aisee des 
images sur CDs. II faut pouvoir extraire des paquets de 500 Mo correspondants a des images 
acquises a partir d'une date donnee. 

- Droits d'acces aux donnees. On doit pouvoir visualiser les images a partir de differents postes 
de travail (mktApp sans le module d ''acquisition modAcq) , cependant il faut distinguer la 
personne qui fait l'acquisition des images (operateur) qui doit pouvoir effacer des images, les 
renommer,... des personnes qui consultant les images qui n'ont pas de droit. 

- Previsualisation : aujourd'hui dans le module browser toutes images ont la meme vignette. 
On pourrait associer une vignette de previsualisation a chacune des images (est-ce pertinent ?) 

- Recherche rapide d'images On voudrait retrouver de facon rapide des images a partir de 
mots-cles. Ces mot-cles peuvent se raporter au nom du fichier image et/ou aux informations 
associees aux images (presentent dans le fihier .cfg) 



5 Modifications pour l'application : mktApp 

- Probleme de l'affichage module d'acquisition/browser/visualisation. Lorsque Ton cree un nou- 
veau module d'acquisition depuis le module browser, celui-ci n'apparait pas en premier plan. 

- Doit-on mettre un "setup" permettant de choisir la carte d'acquisition, l'affichage de certains 
elements (mode expert ou mode utilisateur par exemple). 

- Assistance au protocole d'acquisition : peut-on integrer un controle actif d'un protocole dans 
l'application ? 

- Utilisation de XML pour gerer l'ensemble des parametres? 

- Veille complementaire : il faudrait regarder des applications similaires a celle que Ton doit de- 
velopper pour s'en inspirer eventuellement. Qu'en est-il du web? De l'insertion dans un milieu 
medical ? 



6 Traitement de l'image : TDI 

Les traitements d'images identifies sont : 



6.1 TTRQTR 

- soustraction temps reel du fond lie a la reflexion parasite. Le module d'acquisition 
permet d'enregistrer le fond lie a la reflexion parasite en d'en calculer une moyenne tempore.]. Ce 
fond est soustrait a l'image courante (soustraction point a point de l'image brute avec la moyenne 
temporelle du fond). Cela permet d'ameliorer le contrasts de l'image courante. L'utilisateur doit 
etre forme pour acquerir le fond. A l'heure actuelle la soustraction n'est valide que lorsque les 
parametres d'acquisition du fond et de l'objet sont identiques (probleme d 'offset electronique, 
meconnaissance de la courbe de gain). On reparle du talon ? : influence par le module de 
visualisation ?. 

- Correction du decalage inter-Iigne : II s'agit d'un mode utilise par defaut (groupe parametres 
electroniques). Evolution avec la nouvelle electronique et le nouveau boitier? 

- Moyennage temporel : Aujourd'hui seul ce type de moyennage peut etre effectue, est-il judi- 
cieux de laisser la possibilite a. l'utilisateur d'appliquer d'autres types de moyennages? 

- Correction automatique du gain : II est. encore peut utiliser aujourd'hui du fait des variations 
de signal lie a 1 'offset electronique. II peut faire perdre une image ou deux lors de l'acquisition. 
Evolution avec le nouveau boitier? 

- Detection des fibres, calibration, reconstruction... : scenario, interface et implementation ? 
Choix avec ou sans detection des fibres? 
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6.2 TREAI 



- Detection des fibres, calibration, reconstruction... : scenario, interface et implementation ? 

- Correction des distorsions geometriques : On peut envisager un module de traitement dans 
lequel on rentre les informations necessaires (cf. programme prototype) et. une image de la mire 
de distorsion. Un champ de deformation est alors calcule (en X, en Y). Ce champ calcule ainsi 
que les informations necessaires sont. alors utilisees pour supprimer les distorsions sur l'image 
d'interet. On peut alors appliquer les outils specifiques de mesures (interactifs) a l'image corrigee. 



7 Question ouvertes 

Module ActiveX ou JavaBeans pour Tintegration de nos modules dans des applications tiers? 
Export/ Acces par le net ? 
Valeur ajoutee ? 
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Manuel de reference du module ImageCell 

Auteur : Aymeric Perchant 
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Diffusion : interne 

Section : informatique, image 

Sujet : Manuel de reference de la partie libMKTProcessing pour imageCell. Ce manuel contient. 

les informations de traitement d'images, de programmation C++, et d'utilisation. 
Version : Revision : 1.1 

Reference du document : $Id: ReferencelmageCell.tex.v 1.1 2002-02-07 16:30:08 aymeric Exp 

Introduction 

Ce document rassemble les different.es informations de reference pour le module principal intitule 
Image-Cell. Ce module contient les parties suivantes de traitement d'images : 

- detection des fibres, 

estimation des (lux revenarit des fibres, 

- corrections des defauts de l'appareil (biais, fond, ...). 

- calibration de l'appareil. 

- reconstruction d'images. 

Ces different.es parties seront abordees en detail dans chaque partie pour expliquer les algorithmes 
et les parametres (chapitre 1), Fimplantation de ceux-ci (2) et leur utilisation pratique (3). 
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Chapitre 1 



Traitement d'images dans Image-Cell 



Introduction 

Cette partie decrit les algorithmes en commencant par la description generale, jusqu'aux boites 
elementaires. 

1.1 Schema global 

La figure 1.1 represente le schema global du traitement. On peut distinguer quatre blocs, dont 
deux sont identiques : a une parametrisation pres. II existe done trois groupes de traitements qui 
sont. les suivants : 

Detection des fibres : ces traitements permettent de detecter et d'isoler chaque fibre sur une 
image, ainsi que d'analyse la structure d'agencement des fibres du guide d'images ; 

Estimation des flux : une image brute de taille 640 x 640 x 20 (largeur, hauteur, nombre 
d'images temporellement )contient 8 megapixels qui represented l'information vue par 10000 
ou 30000 fibres. Ce bloc permet d'isoler l'information effectivement vue par chaque fibre ; 

Calibration et reconstruction : La connaissance de l'information vue par chaque fibre est 
ensuite traitee pour etre reconstruite sous la forme d'une image debarrassee des defauts de 
l'appareil. 

Chacun de ces trois blocs sont maintenant. detailles. Le prototypage des algorithmes decrits ici est 
detaille dans le rapport de stage de Sandra Marti [?] ; nous renvoyons le lecteur a cette reference 
pour plus de renseignements sur la demarche de developpement des algorithmes, notamment celui 
de detection des fibres. 

Le bloc de gauche d'estimation des flux permet de mesurer le flux sur une image qui represente un 
objet aux proprietes constantes dans Fespace (un milieu diffusant homogene ou un objet presque 
homogene en mouvement aleatoire, un miroir). La sortie de ce bloc est done une image des taux 
d'injection dans chaque fibre, et qui en prend en compte l'integralite de la chaine. 
Le bloc de droite d'estimation des flux permet de l'estimer pour un objet a observer, et a travers 
l'appareil. Le dernier bloc de calibration est une division de l'image des taux d'injection par l'image 
de l'objet observe. Cette operation permet de compenser les mauvaises injections dans certaines 
fibres. Puis l'image est reconstruite. 

1.2 Bloc de detection des fibres 

La detection des fibres s'organise autour de quatre traitements successifs : 
- pretraitements, 

5 




Fig. 1.1 - Schema global du traitement. Les cadres gras representent les entrees. Les cadres poin- 
tilles representent les donnees ou les traitements facultatifs (biais. soustraction du fond...). 

- ligne de partage des eaux (LPE) : la premiere segmentation par region, 

- separation des fibres trop grosses (split), 

- fusion des fibres trop petites (merge). 

Les deux dernieres etapes peuvent Stre bouclees; neanmoins une analyse de performance et de 
convergence doit alors etre menee. Nous n'aborderons pas ce sujet pour le moment. 

1.2.1 Pretraitements 

Les pretraitements realises sont les suivants. 

Diffusion anisotrope : le but est de lisser l'image dans les zones plates, c'est a dire les zones 
inter fibres. Le prototype fonctionnait. mais n''utilisait pas la librairie inrimage, mais une autre 
libriairie (celle utilisee par Sandra, issue de l'ENST). Le reglage de la diffusion est laissee a 
plus tard, lors de la prochaine release de la liblnrimage qui contiendra de la documentation 
a ce sujet ! Parametre [_doDif fusion.] . 

Interpolation x2 au plus proche voisin : on cherche a simuler des elements structurant de 
morphologie mathematique avec un rayon inferieur a un. L'image est doublee pour que l'ou- 
verture qui suit ne touche pas aux maxima isoles, mais seulement ceux qui sont 8-connexe, 
mais non 4 connexe (voisins par une diagonale). L'interet est de faire une selection des maxima 
elimine par l'ouverture. Parametre [_doZoom_] . 
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Ouverture numerique : on cherche a eliminer les maxima parasites si rues sur les fibres. C'est 
un pretraitement classique de la LPE, mais eventuellement modifie par le comportement du 
doublage de l'image explique precedemment. 

Inversion de l'image : la LPE fonctionne a partir de minima, on inverse done l'image pour 
transformer les minima en maxima. 

L'operation d 'interpolation et d'ouverture qui suit devrait pouvoir etre remplace par un nouvel 

operateur de morphologie mathematique qu'il faudrait definir. Cela permettrait un gain important. 

puisque toutes les operations suivantes se feraient sur une image de taille inferieure. 



Le watershed de la liblnrimage necessite la definition de marqueurs dans l'images. Nous obtenons les 
marqueurs comme minima de l'image pretraitee avec de h-domes [?]. Les h-domes pour les minima 
som^efinfsnpl)ur~ura 



avec E%°(b) la reconstruction geodesique par erosion de b dans a. g designe ici Pensemble des h- 
domes, que Ton seuil a 1 afin de tous les recuperer. On obtient ainsi les minima de profondeur au 
moins egale a h. Ce parametre h est donne par [_hDomeWatershed_] . 

Le resultat de cette LPE est une image de composantes connexes representant chaque fibre detectee. 
Afin de resoudre bon nombre de problemes de bord on colle les composantes touchant le bord au 
bord lui-meme. Ces composantes sont considerees comme faisant partie du bord, et du fond et sont 
mise a la composante connexe d'indice 0, qui est l'indice reserve au fond, et au bord. 

Enfin, les composantes connexes sont reorganises pour ne pas presenter de trous, et sont tries par 

ordre decroissant de taille. 

1.2.3 Split 

Les pretraitements, et notamment 1'ouverture numerique servent. a limiter le defaut principal de la 
LPE qui est de sursegmenter. En faisant cela, on a tendance a augmenter le nombre de segments 
qui vont englober plusieurs objets (ici, des fibres). Le premier defaut que Ton corrige est la sous- 
segmentation inevitable au traitement, mais surtout rajoutee par ces pretraitements. Pour cela on 
effectue les traitements suivants. 

Selection des segments suspectes de sous-segmentation : on trouve les fibres dont la taille 
est superieure a un seuil fixe en fonction de la moyenne des tailles normalisees (parametre 
[_oversizedFactor_] ),et en fonction du nombre de voisins (parametre [_overSizedNeighLimitJ ). 

Re-segmentation de ceux-ci : Les segments selectionnes sont isolees et re-segmenter soit sur 
l'image initiale sans pretraitement, soit sur l'image de carte de distance a l'interieur de ces 
regions (parametre LoverSizedResegmentMetliodJ). Dans le premier cas, on suppose que 
les defauts majoritaires proviennent des pretraitement, dans le second cas on suppose qu'ils 
etaient la avant, et que Ton souhaite separer les segment aux endroits de retrecissement 
(forme de la cacahuete). La LPE est effectuee avec les meme parametres que la premiere. 

1.2.4 Merge 

Puis on cherche a corriger le defaut le plus classique (pour une LPE) : la sursegmentation. On 
cherche done a fusionner les segments. Pour cela on effectue les operations suivantes : 

- preselection des candidats a la fusion, 

- parmi ceux-ci trier ceux qui seront obligatoirement fusionnes. 

- pour les autres, eliminer les fusions impossibles 
- repertorier les fusions possibles 



1.2.2 Premier watershed 



g = Ef{f + h)-f, 



(1.1) 
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- eliminer les fusions dormant de mauvais resultats 

si il ne reste pas de fusion apres filtrage, retirer la fibres des fibres a fusionner 
- pour toutes les fibres rest.antes, prendre la meilleur fusion au sens de la compacite des segments 
fusionnes. 

Les points importants sont detailles ci-apres. 
Remarques sur la compacite 

La compacite utilisee est le critere le plus simple existant. Pour un objet de perimetre P et de 

surface S, le critere est Le perimetre est pris comme etant la frontiere de l'objet 4-connexe, et 

done est defini comme tous les pixels de l'objet ayant un voisin 8-connexe avec le fond (defini ici 

comme etant tout le reste de l'image. La surface est simplement. l'objet lui-meme. 

Cependant, ce critere simple a calculer n'est pas invariant, avec Pechelle de Pobjet. On peut montrer 

quepbllr^e^fornWsirnp^ 

Cette croissante n'est cependant qu'en moyenne en ^, avec n la diametre de l'objet (ici, pour une 
boule en 4-connexite). Ce critere reste done tout de meme assez raisonnable, mais plus delicat a 
manier dans le cas de la comparaison de la compacite d'objets de faille differentes. 
II faut remarquer ici que les premiers essais sont realises sur une image deformee par les distortions, 
et qu'il faudrait, pour etre plus precis soit corriger l'image avant, soit tenir compte de l'anisotropie 
non stationnaire spatialement : on a tendance a surestimer a la fois les frontieres et la surface sur 
les bords de l'image. 



Selection des fibres detectees trop petites 

La selection s'effectue en deux etapes. On commence par trouver les segments candidats a la 
fusion avec un seuil sur la taille. On utilise le parametre [_underSizedFactor_] qui est un facteur 
multiplicatif de la moyenne, et est done inferieur a 1 en general. La seconde etape est de trouver 
parmi ces segments lesquels doivent obligatoirement etre fusionnes. et lesquels doivent etre filtres. 



Les segments obligatoirement fusionnes 

La selection des segments (ou fibres) obligatoirement fusionnes s'effectue dans Pespace (nombre de 
voisins) x (taille normalisee). La taille est normalisee par rapport a la moyenne (on divise toutes 
les tallies par la moyenne des tailles). Le nombre de voisins est compte avec la 8-connexite. 
Les figures 1.2 et 1.3 represented Phistogramme conjoint de la taille normalisee et du nombre de 
voisins pour chaque fibre. On remarque qu'il existe un axe principal qui est centre sur la droite 
passant par des points d'interets particuliers (6, 1).. (8,2), et (5, .5). Ces points peuvent s'explique 
geometriquement, comme indique sur la figure 1.4 

Cet axe est done un axe naturel de la structure hexagonale des fibres dans cet espace. En s'orien- 
tant sur cet axe, les petites fibres se situent vers un nombre de voisins faibles et une taille 
faible. On determine une frontiere lineaire pour selectionner les fibres detectees reellement trop 
petites. Apres quelques tests sur des images reelles de la base de donnees «Novembre 2001» et 
en fonction des resultats preliminaries de [?], on a choisi une droite passant par les points (6; |) 
(3; 1). Les parametres de cette droites sont donnees sous la forme y = mx + p, avec m la pente 
[_underSizedMendatoryMergeSlope_] et p l'ordonnee a Porigine [_underSizedMendatoryMergeOf f set_] . 
Tout point en dessous de la droite sera fusionne obligatoirement, et ne va done pas passer par les 
filtres decrits dans la section suivante. Les autres fibres sont par contre filtree comme suit. 

Filtrage des fusions possibles 

On entend par filtrage ici un procede qui retire certaines fusions. II existe trois filtres possibles. 
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Effectif 




Ftg. 1.2 - Vue 3d de l'histogramme 2d conjoint de la taille normalise et du nombre de voisins 
pour chaque fibre. Des lignes de niveaux sont projetees sur la fond du graphique : les lignes sont 
centrees autour de 6 voisins et de la moyenne de la taille des fibres. 



Filtrage sur la taille de la fusion : (parametre [_underSizedFilterSize_]) on elimine les 

fusions dont la taille totale depasse un seuil fixe par le parametre [_underSizedFilterSizeMax_] , 
qui est pris comme facteur multiplicatif de l'ecart type ajoute a la taille moyenne. 

Filtrage sur la compacite totale : (parametre [_underSizedFilterCompact_]) on elimine les 

fusions dont la compacite totale depasse un seuil fixe par le parametre LunderSizedFilterCompactMaxJ , 
qui est pris comme facteur multiplicatif de l'ecart type ajoute a la taille moyenne. 

Filtrage sur le changement de compacite : (parametre [_underSizedFilterMoreCompact_] ) 
on elimine les fusions qui augmentent la compacite de l'objet avec lequel on veut fusionner. 

Dans les parametres par defaut, on n'utilise pour le moment que le premier filtre qui semble suffisant 

dans la plupart des cas. Les autres filtres semblent trop restrictifs et doivent etre encore etudies. 

1.3 Bloc d'estimation des flux 

Le bloc d'estimation des flux comprend plusieurs sous-parties : 

estimation du flux vu par chaque fibre : c'est ce que Ton peu recuperer d'une image a l'aide 

de la detections es fibres ; 
estimation du fond de l'image : ici le fond designe les reflexions parasites ou (ou inclusif !) 

Foffset du a l'electronique et au detecteur; 
soustraction du fond : on retire le fond a l'image ; 
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Norabre de voisins 
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Fig. 1.3 - Lignes de niveau de Thistogramme 2d conjoint de la taille normalises et du nombre de 
voisins pour chaque fibre. La forme 3d de cet histogramme est donnee sur la figure 1.2 




Fig. 1.4 - Variation du nombre de voisins en fonction de la taille du segment detecte comme etant 
une fibre 

correction du biais : les images sont en general biaises, c'est a dire qu'il existe un fond lente- 
ment variable au lieu d'etre constant. 

Le bloc de soustr action du fond est optionnel, ainsi que celui de la correction du biais. 
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1.3.1 Estimation du flux vu par chaque fibre 

Le rapport. [?] decrit une methode pour rattraper les problemes de saturation de l'image. Ces 
probleme vont disparaitre avec l'arrivee de l'electronique mi mai 2002. Nous avons done decide de 
ne pas continuer dans cette voie pour le moment. Ici nous avons simplement choisi d'estimer le 
flux avec la moyenne sur la fibre. 
Plusieurs ameliorations sont envisageables : 

- tenir compte du bruit poissonien dans l'estimateur du max avec la moyenne, 

- tenir compte de la disparite de forme de fibres du guide d'image ; on peut constater que les toutes 
les fibres ne sont pas identiques et que l'estimateur du flux devrait s'adapter a la forme exacte 
de la fibre, 

- verifier les problemes possibles dues a la zone inter-fibre, en fonction de la tache focale sur l'entree 
du guide, 

- et plus geiu'i aU'irieni. ivcuperer le modele complet d'injection et do ret our du spot de crit d ans 
Ces' ameliorations peuvent etre etudiees, mais ne sont, pas prioritaires pour le moment. 



1.3.2 Estimation ou utilisation du fond, puis soustraction 

Le fond peut avoir plusieurs sources; celles-ci sont decrite dans les memos [?, ?, ?]. En resume, 
nous appelons fond ici soit les reflexions parasites sur les optiques, et done y compris sur la sortie 
du guide d'images, mais egalement 1 'offset due a la chame de numerisation. 

Si 1'ofFset est dominant sur 1'image, on ne peut pas obtenir le fond simplement. en retirant l'image. 
car l'offset depend du contenu, et n'est done plus le meme. Dans ce cas, on utilise un quantile de 
Fhistogramme ou j~) pour Pestimer. 

Dans le cas contraire, il faut utiliser le fond acquis lorsqu'on retire l'objet a regarder, et. le soustraire. 
Dans tous les cas il faut bien penser a saturer la soustraction pour ne pas etre gene par des outliers 
negatifs. L'utilisation d'un talon (cf. [?]) n'est pas retenue car le fond diffus n'est pas du tout une 
composante majoritaire du signal. 



1.3.3 Correction du biais 

Le biais doit a priori etre corrige sur l'acquisition des taux d'injection (branche de gauche) et sur 
l'acquisition de l'objet (branche de droite). Pour le premier cas, cela vient du fait que la calibration 
se fait sur un miroir plan, et que la courbure de champ va reduire la qualite d'injection au retour 
sur les bords (qui sont defocalises) . Ce ne serait probablement pas le cas 

Sur l'objet, ou dans un milieu diffusant homogene, l'injection reste moins bonne sur les bords, et 
cela se traduit par un biais tres similaire quant a sa forme au premier. 
Dans tous les cas le biais a une symetrie quasi circulaire. 

L'estimation du biais se fait en divisant l'image en N x N blocs de taille fixe, puis en estimant le 
biais sur chaque bloc. Pour cela il faut considerer la nature de l'objet observe. Dans le cas d'un 
objet homogene, le biais peut etre acquis en prenant la valeur moyenne ou mediane sur le bloc. 
Quand il y a un objet, il faut savoir si cet objet, est plus sombre ou plus claire que le biais. Dans 
notre cas, le biais est. multiplicative, et on prend done plutot un operateur de moyenne ou de 
mediane (par rapport a un max ou min pour un biais additif). 

On obtient alors une image de taille N x N qui est utilisee avec une interpolation pour trouver la 
valeur du biais vue par chaque fibre. On utilise pour le moment une interpolation bilineaire, faute 
de mieux. La librairie inrimage pourrait etre complet.ee pour integrer une interpolation en spline 
cubique point a point. 

On peut. remarquer que le biais peut etre estimee plus finement avec d'autres methodes plus subfiles 
qu'il faudrait. etudier. Notre methode reste cependant rapide et semble suffisante. 
Une fois le biais estime, il faut diviser l'image par son biais. II faut. penser a verifier que le biais 
est. bien superieur a un pour ne pas multiplier en fait ! Dans le cas de l'image de calibration, on va 
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retablir une dynamique Active sur 1000 niveaux afin de ne pas avoir de problemes avec la seconde 
division pour la calibration du taux d'injection (cf. partie 1.4.1). 



1.4 Bloc de calibration et reconstruction 

1.4.1 Calibration des taux d'injection 

Le bloc d'estimation du flux de gauche permet d'obtenir une image des taux de calibration fibre 
a fibre. On peut dont Putiliser pour retablir une injection de 100% sur toutes les fibres. Pour cela 
on divise Fimage de l'objet par Pimage du taux d'injection, apres avoir pris quelques precautions 
sur les valeurs de la calibration (superieures a 1, avec suffisamment de dynamique). 

1.4.2 Reconstruction mosaique 

La reconstruction mosai'que s'effectue en repartissant sur toute la surface de chaque fibre la valeur 
estimee du flux prise apres Pimage de calibration. 

1.4.3 Reconstruction par RBF 

A faire... 

1.5 Resultats et discussion 

Les resultats sont visibles dans la base de donnees d'images... 

1.5.1 Exemple 

Les figures 1.5 et 1.6 illustrent toutes les etapes du traitement de l'image. Dans ce cas de figure, le 
fond de l'image ^ecalHwafiondes fibres est"^trmg;"<re pourqudi" irri'appafait pas"surles"ifnages": 
Le fond de l'objet est par contre fourni. Les images de calibration et l'objet sont toutes les deux 
de-biaisees. Ce resultat a ete obtenu avec les parametres par defaut. Sur cet exemple ; on remarque 
que les images des biais n'ont pas la symetrie circulaire habituelle. La raison est que la majorite 
du biais venait en fait de la saturation du detecteur qui se traduisait par un nouveau fond additif 
non prevu au depart, mais qui est tout de m§me attenue par les traitements. 

1.5.2 Discussion 

De nombreuses idees d'ameliorations ont ete proposee dans ce chapitre. Certaines semblent plus 
importantes que d'autres. Voici la proposition d'un tri de ces idees par ordre decroissant d'impor- 



Reconstruction par base de fonctions radiales : cela permet. d'obtenir une image lisse et de 

faire d'eventuels traitements ulterieurs. 
Supprimer le zoom sur l'image : le zoom est utile pour diminuer 1'effet de Touverture nu- 

merique avant la LPE. II faut trouver l'operateur equivalent sur l'image non zoomee. Une 

hypothese a creuser est que ce serait equivalent a une ouverture avec un element structurant. 

en croix ( ?) . 

Auto-contr6le de la detection des fibres : savoir si la detection est correcte, ou si au cours 

des acquisition on arrive a un moment ou il faut la refaire. 
Auto-contr6le de valeur du fond : savoir sir le fond propose est valable ou si il faut l'estimer. 

Dans le cas ou le fond n'est. pas fournit, en faut-il un? 
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tance. 



(e) Image de-biaisee issue de (d) 



(f) Resultat final 



Fig. 1.6 - Etapes intermediaires de la procedure Image-Cell 
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Amelioration continue de la detection des fibres : on devrait pouvoir utiliser les images 

des objets pour renforcer ou corriger la detection des fibres. 
Estimer le biais avec des splines : le biais est estime par reechantillonnage bilineaire qui pro- 

duit un effet de blocs qui peut etre genant 
On peut aussi prevoir plusieurs axes de recherche pour des fonctionnalites futures : ce sera le sujet 
d'un prochain memo... 
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Chapitre 2 

Programmation d'Image-Cell 



Introduction 

Cette partie rassemble les principaux elements de la librairie. La description de ces element est 
globale, cette documentation venant en complement de la documentation automatique 1 decrivant 
toutes les classes et toutes les fonctions. et des tests unitaires decrivant Futilisation des classes et 
fonctions principales. 

2.1 Une librairie sur les graphes 

La programmation des algorithmes a necessite Tutilisation de structure de graphes. Nous avons 
done recherche des librairies C++ permettant de les manipuler. Nous avons teste les librairies 
suivantes. 

GTL : pour Grap h Template Li brary (http:// www.infosun.f mi.uni-passau.de /GTL/). C ette 
librairie est la plus adaptee a nos besoins et a ete celle utilisee pour le prot6typage~c[ans 
[?]. Le probleme principal de cette librairie est son prix dissuasif pour les entreprises : 5 000 
Euros (licence de developpement pour un seul poste). Elle etait gratuite pour le stage car 
utilisee par une etudiante... Nous n'avons done pas choisi de garder cette librairie pour cette 



LEDA : pour Library of Efficient Data Types and Algorithms (http : / /www. mpi-sb.mpg.de/LEDA/). 
Cette librairie traite de nombreux probleme de type de donnees, dont celui des graphes. Cette 
librairie a cependant beaucoup evoluee vers une sorte de librairie couteau Suisse qui fait beau- 
coup de choses, dont les interfaces graphiques. et qui redefinit de nombreux outils maintenant 
standard, notamment les types STL. II existe done de nombreux cas de conflits potentiels 
avec nos softs. Le prix de cette librairie est de 1500 US$ pour un developpeur, et 8600$ pour 
developper sur tout un site, et vendre des produits derriere. J'ai personnellement utilise cette 
librairie pour ma these (gratuitement car pour une these...). 

BOOST : une sorte de groupware qui travail sur des proposition de futurs standard C++ 
(http://www.boost.org/index.htm). II y a en fait toute un librairie avec de nombreux 
themes. Le point central est que e'est de la programmation avancee, generique, reutilisable, 
versatile. C'est beau mais e'est tres complique a utiliser malgre une documentation tres 
propre. Le probleme vient des pre-requis indispensables en programmation objet de haute 
voltige qui sont un peu trop complexes pour nos besoins. et qui necessiterai un travail sup- 
plemental important. Cette librairie est cependant gratuite, et tres bien faite. Un attrait 
interessant est qu'elle n'est constitute que de fichiers .h, sans partie a compiler. Tout repose 

'Utilisant doc++ 
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sur une gestion complexe de templates qui ne peuvent pas etre precompiler. II n'y a done pas 

de problemes de compilation sur differentes plateformes. 
Toutes ces librairies fonctionnent sous toutes les plateformes courant.es. et avec tous les compilateurs 
courant. De plus elles sont maintenues et stables. II existe de nombreuses autres petite libriaires a 
droite a gauche mais qui sont loin d'approcher la qualite de celles-ci. Nous n'avons cependant pas 
pu en retenir une celle, soit a cause du prix et. des risques lies a Tutilisation d'une grosse librairies 
supplemental, soit a cause de la complexity d'utilisation trap grande. 

Nous avons done choisi de prendre nos propres structures de donnees. adaptes aux traitement 
pour plus de rapidite a la fois de developpement et d'execution. Ces structures de donnees sont 
cependant assez souples car elle fonctionnent sur des types et des algorithmes STL qui peuvent 
etre interchangeables assez rapidement... 



2.2 — Gestion generate des parametres d'algorithmes 
2.2.1 L'objet parametre 

L'objet parametre n'existe pas encore en tant que tel. mais va bientot, voir le jour... II est pour le 
moment a l'etat de concept dans le meme sens que la terminologie STL. 
Une classe parametre est de type suivant : 

class AlgoParam { 

public : 



Ad j GraphPar am ( ){ 
def aultsO ; 

}; 

void defaults () { // Cette fonction fixe les parametres par defaut 
} 

/// Load parameters from a file. If any error occurs , default values are provided. 
void - ioadPar am ( c onst~cbar + fname)-; 

/// Save parameters in a file, 
void saveParam( const char * fname) ; 

/// Output stream 

friend std: : ostream & operator « (std: :ostream &, AlgoParam &) ; 
/// Input stream 

friend std: : istream & operator » (std: : istream &, AlgoParam &) ; 
public: // parameters list here, with comments please ! 
[...] 
}; 

Les fonctions loadParam et saveParam vont. juste creer un flot. dans le fichier dont le nom est 
specifie. lis deleguent done le travaille aux flots ostream et istream. 
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Pour le moments les flot.s sont simples, et. ne font que ecrire le nom du parametre suivie de la valeur 
numerique, le tout sur une ligne. Les parametres sont encadres par 

ALGDPARAM_CONFIGURATIQN 
PARAMETERS 

et 

END_ALGOPARAM_CONFIGURATION 

Ce mecanisme (surement perfectible...) permet de deleguer les entrees et les sorties aux flots qui 
peuvent done etre imbriques. Ainsi un objet parametre peut contenir un autre objet parametre 
sans avoir a recrire les entrees et les sorties de cet objet. 

De telles classes de parametres sont/utilises pour- ImageGell FiberDetector Adj Graph et FiberFlow. 



2.2.2 Le futur des parametres 

La premiere chose a faire est de creer une classe virtuelle (e'est-a-dire une interface) pour la classe 
de parametres. 

Enfin, dans un second temps est a l'etude Futilisation d'XML pour interfacer ces parametres. 

2.3 La detection des fibres : FiberDetector et Adj Graph 

Voir le test unitaire test_f iberdetect pour les details d'utilisation. 

L'objet FiberDetector est un objet-algorithme modulaire pour detecter les fibres. II contient un 
objet specialise AdjGraph qui permet de faire Fetape de fusion des fibres qui est la plus compliquee. 

2.3.1 Fonctionnement detaille 

On retrouve les etapes suivantes dans la classe : 

void init(Inr image *inrln) ; 
void preprocessingO ; 

_vpid water she d () ; 

int splitBiggest () ; 
int mergeSmallest () ; 

L'etape d'initialisation init permet de fixer l'image servant de base a la detection des fibres. Puis 
l'etape preprocessing fait es premiers traitements : zoom, ouverture numerique, inversion de 
l'image, et quelques calculs preliminaire, comme une structure de voisinage des fibres. Ces deux 
etapes permettent de preparer l'image pour les traitements de segmentation. 
La segmentation est. realisee par une ligne de partage des eaux (LPE) dans la fonction watershed. 
L'algorithme pourrait egalement interfacer un autre algorithme... 

Les deux etapes suivantes sont deux etapes de post-processing de split and merge. splitBiggest 
va reperer les fibres trop grosses, puis les resegmenter. mergeSmallest permet de fusionner les 
fibres les plus petites. Cette derniere etape plus complexe et plus critique est elle-meme encapsulee 
dans une classe algorithme AdjGraph expliquee dans la prochaine partie. 

Ces deux etapes contiennent une dizaine de parametres qui permettent d'ajuster la precision et la 
sensibilite de l'algorithme. La structure de l'algorithme permet d'eventuellement. boucler ces etapes 
afin de faire converger l'algorithme vers une segmentation stable au sens de l'idempotence du split 
and merge, pour un jeu de parametre donne 2 . 

2 Nous n'avons pas de preuve de l'existence d : une telle convergence et il semblc difficile a priori d'en obtenir. 
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2.3.2 L'etape de fusion avec AdjGraph 



Cette partie decrit l'etape de fusion. Cette etape est importante car les petites fibres (sur-segmentees) 
sont plus ggnant.es que les grosses (sous-segmentees) sur un resultat final. La complexity de cette 
partie nos a oblige a en faire une partie modulaire. au meme titre que le FiberDetector. Les 
differentes «modules» ont ete les suivants : 

- type de calcul de la compacite d'une region, et de la fusion de deux regions 

- construction du graphe d'adjacence dedie 

- structures de donnees dediee aux calculs des attributs utiles. 

Cette etape etait la plus longue sur le prototype de [?], et a done fait l'objet d'optimisation. L'etape 
de fusion n'est realise que sur des criteres de formes, sans lien directe necessaire avec l'image de 
depart 3 . Le probleme principal est le calcul des attributs utilises : la compacite et la surface. La 
principale difficulty est liee a ce que le graphe va evoluer par fusion de regions, et que Ton cherche 
a mettre a jduirliies afinbuts ;l^^f~^alemeh£'fest€ar" l«"fuaohs 'pMables"et"'db"hc pre calculer les 
attributs possibles. 

Pour cela, on distingue les frontieres exterieures et interieures d'un objet comme indique sur la 
figure 2.1. On divise la frontiere d'un objet en fonction de l'adjacence. La frontiere de A est egale a 
la somme des frontiere interieures de A avec tous ses voisins. La frontiere interieure de A commune 
avec F est egale a la frontiere exterieure de F commune avec A. Cette derniere propriete permet 
de ne parler que de frontiere interieure. 



Fig. 2.1 - Les frontieres exterieurs et interieures d'un objet. La frontiere interieure de A commune 
avec F est egale a la frontiere exterieure de F commune avec A. Ces informations sont portes par 
des arcs du graphe d'adjacence. 

Ces informations sont portes par des arcs du graphe d'adjacence. Un arc allant de A vers B va 
porter l'information de frontiere commue que A detient avec B. 

On peut formaliser ces relations ainsi. Soit G = (N, E) un graphe avec N l'ensemble de ses noeuds 
et E C N x N l'ensemble de ses arcs. G est un digraphe 4 ou deux nceuds n'ont soit pas d'arcs 
en commun, soit deux arcs antisymetriques. E est done un ensemble de paires ordonnees. Soit s : 
N -* R + une fonction permettant d'obtenir la surface. Soit f : E —> R + une fonction permettant 
d'obtenir la frontiere interieure entre les deux nee designe par Fare. Enfin soit c : N — > R + la 
fonction permettant d'obtenir la compacite d'une region. On note V(n), n € N l'ensemble des 
voisins de n. La compacite peut alors s'ecrire : 





(2.1) 




graphe ou 



fait 




insion (lourde !) possible 
sont diriges 



20 



Document, confidentiel - Mauna Kea Technologies - 2002-02-07 



Pour tester la compacite de la fusion eventuelle de n et m, on note la nouvelle compacite c(nUm) : 
_ ( S oe /((n, o)) + E o6 y (m) /(K o)) - /((», ™)) - f((m, n))) ' 



Et si la fusion d 
an): 



s(n) + s(m 

i avec n est effectuee. on modifie les fonctions (on suppose que m est rajoutee 



s(n)^s(n)+s(m) 
f(n, o) <- f(n, o) pour o € V(n) et o £ V(m) 
f{n, o) <- f(n, o) + /(to, o) pour o 6 V(n) n V(m). 



On peut ainsi ne travailler que sur les graphes, sans 



"Tecalculer les attributs 1orsnd'une~iusion— i 
voire eventuellement de tester toutes les fusions possibles. 



(2-3) 
(2.4) 
(2.5) 



2.4 



Structures de donnees associees aux fibres : FiberStruc- 
ture et Fiberlnfo 



Voir le test unitaire test.f iberlnf o pour les details d'utilisation. 

2.4.1 Lien entre la structure des fibres et It 



fibre 

La structure de donnee est specifique au probleme qui nous interesse. Les informations qui nous 
interessent sont des informations attachees aux fibres seulement, les fibres etant. en fait liees a la 
structure detectee. Les deux types d'informations sont done etroitement liees, les informations sur 
les fibres ne pouvant exister independamment de la structure des fibres. 

Modele sujet-observateur 

avons done defini ce mecanisme dans le fichier Observer. h. Le principe est le suivant : 

«Definit une interdependance un a plusieurs, de fagon telle que, quand un objet change 
d'etat, tous ceux qui en dependent en soient notifies et automatiquement mis a jour.» 

Ainsi on distingue le sujet, qui pour nous est la structure des fibres detenue par la classe FiberStructure ; 

et les observateurs qui sont les informations que Ton peut attacher aux fibres (classe Fiber Inf o<T>). 

On a ajoute un comportement supplementaire : un objet Fiber Inf o<T> ne peut exister sans un 

FiberStructure, ce pourquoi il n'existe qu'un constructeur qui prend en argument la FiberStructure. 

L'objet Fiber Inf o<T> va d'ailleurs s'initialiser tout seul en regardant la structure a laquelle il est 

attache. 

Quand la structure vient a changer, tous les Fiber Inf o<T> en sont informes si besoin, et ceux-ci 
vont alors se reinitialiser automatiquement. Quand la structure est detruite, les objets d'information 
sont informes egalement qu'ils ne sont plus valides. 



Structure codee comme un graphe : FiberStructure 

La classe fiberstructure se construit a partir de Timage des composantes connexes. L'objet va alors 
garder en interne une copie de cette image de reference. Elle contient en interne la structure de 
graphe associee alors aux composantes connexes. Elle contient egalement d'autres informations, 
comme les barycentres des regions, en coordonnees pixel, les surface des regions ... 
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Toutes ces informations sont accessibles depuis l'interface de la classe, mais ne sont pas modifiables, 
car elles dependent de la structure meme. 

La classe d'information : Fiberlnfo 

La classe Fiber Inf o<T> est une classe template car les informations peuvent etre de nature tres 
differentes. Elle contient un vecteur STL d'information dont la taille est geree automatiquement 
en fonction de la structure associee. Ce vecteur est un membre public afin d'avoir des acces rapides 
aux informations. La classe sert juste de decorateur 5 pour pouvoir gerer les dependances avec la 
structure des fibres. 

Quand le type du template est numerique. il existe des fonctions pour additionner, soustraire. 
diviser et multiplier deux Fiberlnf o<T>, ou bien un Fiber Info<T> avec un element de T : addFI. 
subFI ; divFI, mulFI. ; ~ ~~ 



2.4.2 Transformation d'images en Fiberlnfo et reciproquement 

Voir le test unitaire test_f iberProcessing pour les details d'utilisation. 

Une image peut etre transformee en Fiberlnf o<float> avec la fonction d'estimation simple des 
flux : 

int f lowEstimation (Inrimage* image, Fiberlnf o<float> & flow, bool bias = false). 
II faut au prealable initialiser de Fiberlnf o<float> flow avec la structure, ainsi Falgorithme 
connait via flow l'image des composantes connexes a utiliser. 

La transformation reciproque est la reconstruction d'image. pour laquelle il n'existe pour le moment 
que la methode de mosaiquage : 

int imageMosaic (Inrimage* image, Fiberlnf o<float> & flow) 

Le principe est similaire a f lowEstimation. La structure est passee via le flow. 

2.5 Soustraction du fond et biais : FiberFlow 

- Voir le test unitaire test_f iberFlow pour les details, d'utilisation. 



2.5.1 Fonctionnement general 

La classe FiberFlow implemente le bloc d'estimation des flux de la partie 1.1, sauf la partie 
de transformation de l'image en Fiberlnfo. La classe comporte deux methodes principales : 
setBackground et set Image, et son comportement est parametre avec la classe de parametre 
FiberFlowParam. La classe permet. de soustraire le fond qui est donnee, ou bien il peut Festimer. 
La classe sert aussi a calculer le biais dans l'image, apres soustraction du fond. 

2.5.2 Le calcul du fond 

Voir le test unitaire test_f iberHisto pour les details d'utilisation. 

Pour calculer le fond, on calcul un quantile de 1'histogramme des fibres. Pour cela il existe une 
classe FiberHisto<T> qui prend un FiberInfo<T> en constructeur, et qui permet. de faire des 
calculs sur 1'histogramme calcule a partir de 1 'information des fibres. Ces calculs sont plus rapides 
que sur l'image elle-meme. 



5 Attention, pas au sens du Design Pattern decorateur [?] ! 
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2.5.3 Le calcul du biais 



Voir le test unitaire test_f iberProcessing pour les details d'utilisation. 

Le calcul du biais n'est en fait pas calcule sur l'image, mais encore une fois sur l'information des 
fibres directement. C'est encore un cas ou on utilise la structure et reformation des fibres. Ce 
calcul est fait dans une fonction : 

int biasEstimate ( Fiberlnf o<f loat> & flowln, 
Fiber Info<float> & bias, 
int nbBlocks, 

biasEstimatorType bet=BET_MEAN) 

L'estirnateur du biais est pour le moment a choisir entre la moyenne (BET_MEAN) et la mediane 
(BET.MEDIAN). : : 

On divise l'espace total des fibres en nbBlocks blocs de taille egale. On parcotirt ensuite l'ensemble 
des blocs et on trie les fibres en fonctions du bloc dans lequel elles tombent. Puis, pour chaqtie 
bloc, on applique l'operateur choisi- d 'estimation du biais. 

On calcule le biais ensuite en interpolant pour chaque fibre uniquement (et non pour toute 1 'image) 
les valeurs des blocs les plus proches. On utilise pour le moment une interpolation bilineaire. car 
c'est la seule disponible avec Inrimage. On fait une proposition pour avoir une interpolation spline 
cubique. 

2.6_. Ameliorations 



- La classe de parametres est k optimiser et a programmer plus proprement. 

- Le biais devrait 6tre interpole avec une spline cubique 

- Une librairie de graphe serait tout de meme la bienvenue. Cela permettrait d'avoir une structure 
commune entre FiberDetector, AdjGrapb., et FiberStructure. ce qui n'est actuellement pas 
exactement le cas. 
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Chapitre 3 



Utilisation d'Image-Cell 





C'est une partie delicate car elle doit contenir le sa\ 
developpement, et les nombreux tests. 


'oir-faire acquis pendant le prototypage. le 






Les petites fibres (sur-segmentees) sont plus genante 

result at final. 


s que les grosses (sous^segmentees) sur un 



3.1 Utilisation de la parametrisation 

II faut ouvrir une fenfitre d'un editeur de texte (notepad; emacs, wordpad...) avec le fichier ima- 
g ecell.txt. Ce fic hier e st relu a chaque utilisation. La liste de s p aramet re. ai nsi que des resumes 
d'indications sont donnees en annexe A. 



3.2 Parametrage de FiberDetector 

3.2.1 Pretraitements 

Dans presque tous les cas, il faut faire le zoom sur 1'image pour que Palgorithme fonctionne cor- 
rectement. 

3.2.2 Parametrage de AdjGraph 

Les petites fibres (sur-segmentees) sont plus genantes que les grosses (sous-segmentees) sur un 
resultat final. 

3.3 Parametrage du FiberFlow pour la calibration 

3.4 Parametrage du FiberFlow pour 1'image a reconstruire 

3.5 Parametrage d'ImageCell : les reste des parametres 
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Annexe A 



Parametres de Image- Cell 



Les parametres internes d'Image-Cell sont modifiables via un fichier de parametres present, dans 
le repertoire source du lancement de Fapplicat.ion, au mgme endroit que tomoscope.cfg. Le fichier 
des parametres imagecell . cf g a la forme suivante : 

27 



| Ligne du fichier | Defaut. Possibles 


Remarques 






IMAGECELLPARAM CONFIGURATION 


Parametres de 1 algorithme ent.ier 


PARAMETERS 




FIBERDETECTORPARA RECONFIGURATION 


Parametres de Palgorithme entier 


PARAMETERS 




doZoom 


1 


{0,1} 


Travail sur l'image doublee (plus precis) 


_doDiffusion_ 


0 


{0,1} 


A regler. 


_hDomeWatershed_ 


1 


{1....10} 


Hauteur des maxima selectionnes pour initia- 
liser les fibres. 


_oversizedFactor_ 


1.8 


[i;2] 


iaille normalisee (par rapport a la moyennej 
au dela de laquelle une fibre va etre divisee. 
-■ Ge critere t j s t utilise avvc le critere- suivant- sur— 






















-la-taiHe-(ET-logique) 




_overSizedNeighLimit_ 


6 


{5,6,7} 


Nombre de voisins au dela duquel une fibre va 
etre divisee. Ce critere est utilise avec le critere 
















_overSizedResegmentMethod_ 


1 


{0,1} 


-piece, en . su (_ g q^ ) ^ 

rescgrnenter dans p S image°orighiale S sanrpre 
traitement 




ADJGRAPHPARAM_CONFIGURAI 


ION 


Parametres pour la fusion des fibres 


PARAMETERS 




_underSizedFactor_ 


0.6 


[0.1; 1] 


Taille normalisee (par rapport a la moyenne) 
en dega de laquelle une fibre va etre une candi- 
date pour la fusion. Ce critere est utilise avec 
le critere suivant sur la taille (ET logique) 


_underSizedMendatoryMergeSlope_ 


-0.22 


[?;?] 


Ce parametre et le suivant determinent 
une droite dans Pespace (nombre de voi- 
sins)x(taille normalisee) delimitant les fibres 
qui doivent obligatoirement §tre fusionnees. 
Mieux vaut ne pas y toucher sans faire de cal- 
cul precis ! La droite par defaut passe par les 
points (6; 1/3) et (3,1). Augmenter la pente 
revient a s'occuper plus des grosses fibres avec 
peu de voisins. 


_underSizedMendatoryMergeOffset_ 


1.667 


[?;?] 


Voir remarque precedente. Ce parametre 
est l'ordonnee a l origine (axe des ordon- 
nees=tailles normalisees). L : augmenter re- 
vient a obligatoirement fusionner plus de 










fibres. 




_underSizedFilterSize_ 


1 


{0,1} 


de la taille totale de la fusion ? 


_underSizedFilterSizeMax_ 


1.17741 


[i;3] 


de fibre normalisee superieure a 


_underSizedFilterCompact_ 


0 


{0,1} 


Filtrage des fibres fusionnees en fonction de la 
compacite 


underSizedFilterCompactMax_ 


2 


[1;10] 


Facteur multiplicatif de l'ecart-type 
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Ligne du fichier 


Defaut 


Possibles 


Remarques 






0 


{0,1} 


Ne garde que les fusions qui rendent les fibres 
fusionnees plus compacte que la fibre voisine 
seule. 




END_ADJGRAPHPARAM_CON 


FIGURATION 


END FIBERDETECTORPARAM CONFIGURATION 


FIBERFLOWPARAM _ CONFIGURATION 


de calibration 


PARAMETERS 




_backgroundMode_ 


1 


{0,1,2} 


0 : pas de soustraction du fond ; 1 : soustrac- 
tion avec le fond fourni. sinon 2 ; 2 : estimation 
du fond sur l'histogramme 




- -_backgroundEstin 1a i ioi iC 1 1 1. 




W 


-] J est-imation-du-fond 






_biasCorrection_ 


1 


{0,1} 


Correction du biais, apres soustarction du 
fond 




' 








et en Y pour son estimation 




_biasCorrectionOperator_ 


1 


{0,1} 


Operateur pour la correction du biais : 0 : 


END FIBERFLOWPARAM CONFIGURATION 


FIBERFLOWPARAM _ CONFIGURATION 


Parametres d'estimation du flux pour l'image 
a reconstruire 


PARAMETERS 




_backgroundMode_ 


1 


{0,1,2} 


0 : pas de soustraction du fond ; 1 : soustrac- 
tion avec le fond fourni, sinon 2 ; 2 : estimation 
du fond sur l'histogramme 


_backgroundEstimationCut_ 


0.01 


[0;0.1] 


Quantile de l'histogramme de l'image pour 
l'estimation du fond 


biasCorrection _ 


1 


{0,1} 


Correction du biais, apres soustarction du 
fond 


_biasCorrectionNbBlocks_ 


4 


{2,. .,64} 


Si correction du biais, nombre de blocks en X 
et en Y pour son estimation 


_biasCorrectionOperator_ 


1 


{0,1} 


Operateur pour la correction du biais : 0 : 
moyenne, 1 : mediane 


END FIBERFLOWPARAM CONFIGURATION 


_tempoFilterTypeFib_ 


3 


{0,..,4} 


Type du filtre temporel pour l'estimation des 
fibres :0 : moyenne ;1 : moyenne coupee; 2 : 
medianne ; 3 : max ; 4 : min 




calibration 


1 


{0,1} 


Effectuer la calibration 




END IMAGECELLPARAM COF 


^FIGURATION 
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Manuel de reference du module ImageCell 



Auteur : Aymeric Perchant 
Date : 2002-03-12 



Diffusion : interne 

Section : informatique, image 

Sujet : Manuel de reference de la partie libMKTProcessing pour imageCell. Ce manuel contient 

les informations de traitement d'images, de programmation C++, et. d'utilisation. 
Version : Revision : 1.2 

Reference du document : $Id: ReferencelmageCell.tex.v 1.2 2002-03-12 10:05:01 aymeric Exp $ 

Introduction 

Ce document rassemble les differentes informations de reference pour le module principal intitule 
Ima ge-Cell . Ce module contient les parties suivant.es de traitement d'images : _ 

.^-.detection, des Jibrea, - ■ ; 

- estimation des flux revenant des fibres. 

- corrections des defauts de Pappareil (biais, fond, ...), 
— calibration de Pappareil. — 

- reconstruction d'images. 

Ces differentes parties seront abordees en detail dans chaque partie pour expliquer les algorithmes 
et les parametres (chapitre 1), Pimplantation de ceux-ci (2) et leur utilisation pratique (3). 

Corrections 

- rajout d'un filtrage passe bas, mars 2002. 
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Chapitre 1 

Traitement d'images dans Image-Cell 



Introduction 



Cette partie decrit les algorithmes en commencant par la description generate, jusqu'aux boites 
elementaires. 



1.1 Schema global 

La figure 1.1 represente le schema global du traitement. On peut. distinguer quatre blocs, dont 
deux sont identiques, a une parametrisation pres. II existe done trois groupes de traitements qui 
sont les suivants : 

Detection des fibres : ces traitements permettent de detecter et d'isoler chaque fibre sur une 
image, ainsi que d 'analyse la structure d'agencement des fibres du guide d'images ; 

Estimation des flux : une image brute de taille 640 x 640 x 20 (largeur. hauteur, nombre 
d'images temporellement )contient 8 mega pixels qui represented l'information vue par 10000 
ou 30000 fibres. Ce bloc permet d'isoler l'information effectivement vue par chaque fibre ; 

Calibration et reconstruction : La connaissance de l'information vue par chaque fibre est 
ensuite traitee pour etre reconstruite sous la forme d'une image debarrassee des defauts de 
l'appareil. 

Chacun de ces trois blocs sont maint.enant detailles. Le prototypage des algorithmes decrits ici est 
detaille dans le rapport de stage de Sandra Marti [?] ; nous renvoyons le lecteur a cette reference 
pour plus de renseignements sur la demarche de developpement des algorithmes, notamment celui 
de detection des fibres. 

Le bloc de gauche d'estimation des flux permet de mesurer le flux sur une image qui represente un 
objet aux proprietes constantes dans l'espace (un milieu diffusant homogene ou un objet presque 
homogene en mouvement aleatoire, un miroir). La sortie de ce bloc est done une image des taux 
d'injection dans chaque fibre, et qui en prend en compte l'integralite de la chaine. 
Le bloc de droite d'estimation des flux permet. de Festimer pour un objet a observer, et a travers 
l'appareil. Le dernier bloc de calibration est une division de 1'image des taux d'injection par l'image 
de Fobjet observe. Cette operation permet de compenser les mauvaises injections dans certaines 
fibres. Puis l'image est reconstruite. 

1.2 Bloc de detection des fibres 

La detection des fibres s'organise autour de quatre traitements successifs : 
- pretraitements, 
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FIG. 1.1 - Schema global du traitement. Les cadres gras representent les entrees. Les cadres poin- 
tilles representent les donnees ou les traitements facultatifs (biais, soustraction du fond...). 

- ligne de partage des eaux (LPE) : la premiere segmentation par region, 

- separation des fibres trop grosses (split), 

- fusion des fibres trop petites (merge). 

Les deux dernieres etapes peuvent etre bouclees: neanmoins une analyse de performance et de 
convergence doit alors etre menee. Nous n'aborderons pas ce sujet pour le moment. 

1.2.1 Pretraitements 

Les pretraitements realises sont les suivants. 

Diffusion anisotrope : le but est de lisser l'image dans les zones plates, c'est. a dire les zones 
inter fibres. Le prototype fonctionnait mais n'utilisait pas la librairie inrimage, mais une autre 
libriairie (celle utilisee par Sandra, issue de l'ENST). Le reglage de la diffusion est laissee a 
plus tard, lors de la prochaine release de la liblnrimage qui contiendra de la documentation 
a ce sujet ! Parametre [_doDif fusion.] . 

Interpolation x2 au plus proche voisin : on cherche a simuler des elements structurant de 
morphologie mathematique avec un rayon inferieur a un. L'image est. doublee pour que Pou- 
verture qui suit ne touche pas aux maxima isoles, mais seulement. ceux qui sont 8-connexe, 
mais non 4 connexe (voisins par une diagonale). L'interet est de faire une selection des maxima 
elimine par Pouverture. Parametre [_doZoom_] . 
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Ouverture numerique : on cherche a eliminer les maxima parasites situes sur les fibres. C'est 
un pretraitement classique de la LPE, mais eventuellement modifie par le comportement. du 
doublage de l'image explique precedemment. 

Inversion de l'image : la LPE fonctionne a partir de minima, on inverse done l'image pour 
transformer les minima en maxima. 

L'operation d'interpolation et d'ouverture qui suit devrait pouvoir etre remplace par un nouvel 

operateur de morphologie mathematique qu'il faudrait definir. Cela permettrait un gain important. 

puisque toutes les operations suivantes se feraient sur une image de taille inferieure. 

1.2.2 Premier watershed 

Le watershed de la liblnrimage necessite la definition de marqueurs dans l'images. Nous obtenons les 
marqueur s comme minima de l'image pretrai tee aver, de h-d 6mes | ?]: Les h-domes pourles -mfarima- 
sont definis pour un entier h, une fonction / : 2Z x 2Z, -> 2Z representant~Timage concerriee : 

g=mf+v-f> M 

avec E%°(b) la reconstruction geodesique par erosion de b dans a. g designe ici l'ensemble des fa- 
domes, que Ton seuil a 1 afin de tous les recuperer. On obtient ainsi les minima de profondeur au 
moins egale a h. Ce parametre h est donne par [_M>omeWatershed_] . 

"Le resultat de cette LPE est une image de composantes connexes repres^mn:faaiiu-eiTbTeid'eteirtree: 
Afin de resoudre bon nombre de problemes de bord on colle les composantes touchant le bord au 
bord lui-meme. Ces composantes sont considerees comme faisant partie du bord, et du fond et sont 
mise a la composante connexe d'indice 0, qui est l'indice reserve au fond, et au bord. 
Enfin, les composantes connexes sont reorganises pour ne pas presenter de trous, et sont tries par 
ordre decroissant de taille. 

1.2.3 Split 

Les pretraitements, et notamment l'ouverture numerique servent a limiter le defaut principal de la 
LPE qui est de sursegmenter. En faisant cela. on a tendance a augmenter le nombre de segments 
qui vont englober plusieurs objets (ici, des fibres). Le premier defaut que 1'on corrige est la sous- 
segmentation inevitable au traitement. mais surtout rajoutee par ces pretraitements. Pour cela on 
effectue les traitements suivants. 

Selection des segments suspectes de sous-segmentation : on trouve les fibres dont la taille 
est superieure a un seuil fixe en fonction de la moyenne des tailles normalisees (parametre 
Lovers izedFactorJ ), et en fonction du nombre de voisins (parametre [_overSizedNeighLimit_] ). 

Re-segmentation de ceux-ci : Les segments selectionnes sont. isolees et re-segmenter soit sur 
l'image initiale sans pretraitement, soit sur l'image de carte de distance a Pinterieur de ces 
regions (parametre [_overSizedResegmentMethod_] ). Dans le premier cas, on suppose que 
les defauts majoritaires proviennent des pretraitement, dans le second cas on suppose qu'ils 
etaient la avant. et que Ton souhaite separer les segment aux endroits de retrecissement 
(forme de la cacahuete). La LPE est effectuee avec les meme parametres que la premiere. 

1.2.4 Merge 

Puis on cherche a corriger le defaut le plus classique (pour une LPE) : la sursegmentation. On 
cherche done a fusionner les segments. Pour cela on effectue les operations suivantes : 

- preselection des candidate a la fusion, 

- parmi ceux-ci trier ceux qui seront obligatoirement fusionnes, 

- pour les autres, eliminer les fusions impossibles 
- repertorier les fusions possibles 
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- eliminer les fusions dormant de mauvais resultats 

- si il ne reste pas de fusion apres filtrage, retirer la fibres des fibres a fusionner 
- pour toutes les fibres restantes, prendre la meilleur fusion au sens de la compacite des segments 

fusionnes. 

Les points importants sont. detailles ci-apres. 
Remarques sur la compacite 

La compacite utilisee est le critere le plus simple existant. Pour un objet de perimetre P et de 
surface S, le critere est Le perimetre est pris comme etant la frontiere de l'objet 4-connexe, et 
done est defini comme tous les pixels de l'objet ayant un voisin 8-connexe avec le fond (defini ici 
comme etant tout le reste de l'image. La surface est simplement l'objet lui-meme. 
Cependant. cc ci here sim ple a calculer n'est. pas invariant, avec 1'echelle de l'obje t. On p ent montrer 
que pour des formes simple, le critere est croissant avec la taille de l'objet au lieu d'etre invariant. 
Cette croissante n'est cependant qu'en moyenne en ^. avec n la diametre de l'objet (ici, pour une 
..boule„en.4-coimexite)^Ce„mt.ere reste done tout de meme assez raisonnable, mais plus d elicat a 
rnanier dans le cas de la comparaison do la compacite d'objets de taille different.es: 
II faut remarquer ici que les premiers essais sont realises sur une image deformee par les distortions, 
et qu'il faudrait, pour etre plus precis soit corriger l'image avant, soit tenir compte de l'anisotropie 

les bords de l'image. 

Selection des fibres detectees trop petites 

La selection s'effectue en deux etapes. On commence par trouver les segments candidats a la 
fusion avec un seuil sur la taille. On utilise le parametre [_underSizedFactor_] qui est un facteur 
multiplicatif de la moyenne, et est done inferieur a 1 en general. La seconde etape est de trouver 
parmi ces segments lesquels doivent obligatoirement etre fusionnes. et lesquels doivent etre filtres. 

Les segments obligatoirement fusionnes 

La selection des segments (ou fibres) obligatoirement fusionnes s'effectue dans l'espace (nombre de 
voisins) x (taille normalises). La taille est normalisee par rapport a la moyenne (on divise toutes 
les tailles par la moyenne des tallies). Le nombre de voisins est compte avec la 8-connexite. 
Les figures 1.2 et 1.3 represented l'histogramme conjoint de la taille normalisee et du nombre de 
voisins pour chaque fibre. On remarque qu'il existe un axe principal qui est centre sur la droite 
passant par des points d'interets particuliers (6, 1), (8,2), et. (5, .5). Ces points peuvent s'explique 
geometriquement, comme indique sur la figure 1.4 

Cet axe est done un axe naturel de la structure hexagonale des fibres dans cet espace. En s'orien- 
tant sur cet axe, les petites fibres se situent vers un nombre de voisins faibles et une taille 
faible. On determine une frontiere lineaire pour selectionner les fibres detectees reellement trop 
petites. Apres quelques tests sur des images reelles de la base de donnees «Novembre 2001» et 
en fonction des resultats preliminaries de [?], on a choisi une droite passant par les points (6; §) 
(3; 1). Les parametres de cette droites sont. donnees sous la forme y = mx + p, avec m la pente 
LunderSizedMendatoryMergeSlopeJ et. p l'ordonnee a Forigine [_underSizedMendatoryMergeOf f set_] 
Tout point en dessous de la droite sera fusionne obligatoirement, et ne va done pas passer par les 
filtres decrits dans la section suivante. Les autres fibres sont par contre filtree comme suit. 

Filtrage des fusions possibles 

On entend par filtrage ici un procede qui retire certaines fusions. II existe trois filtres possibles. 
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Effectif 




Fig. 1.2 - Vue 3d de l'histogramme 2d conjoint de la taille normalisee et du nombre de voisins 
pour chaque fibre. Des lignes de niveaux sont projetees sur la fond du graphique : les lignes sont 
centrees autour de 6 voisins et de la moyenne de la taille des fibres. 



Filtrage sur la taille de la fusion : (parametre [_underSizedFilterSize_]) on elimine les 

fusions dont la taille totale depasse un seuil fixe par le parametre LunderSizedFilterSizeMaxJ , 
qui est pris comme facteur multiplicatif de l'ecart. type ajoute a la taille moyenne. 

Filtrage sur la compacite totale : (parametre [_underSizedFilterCompact_] ) on elimine les 

fusions dont la compacite totale depasse un seuil fixe par le parametre [_underSizedFilterCompactMax_] , 
qui est pris comme facteur multiplicatif de l'ecart type ajoute a la taille moyenne. 

Filtrage sur le changement de compacite : (parametre [_underSizedFilterMoreCompact_] ) 
on elimine les fusions qui augmentent la compacite de Pobjet avec lequel on veut fusionner. 

Dans les parametres par defaut, on n'utilise pour le moment que le premier filtre qui semble suffisant 

dans la plupart des cas. Les autres filtres semblent trop restrictifs et doivent etre encore etudies. 

1.3 Bloc d'estimation des flux 

Le bloc d'estimation des flux comprend plusieurs sous-parties : 

estimation du flux vu par chaque fibre : c'est ce que Ton peu recuperer d'une image a l'aide 

de la detections es fibres ; 
estimation du fond de l'image : ici le fond designe les reflexions parasites ou (ou inclusif !) 

Poffset du a l'electronique et au detecteur ; 
soustraction du fond : on retire le fond a Pimage ; 
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FIG. 1.3 - Lignes de niveau de l'histogramme 2d conjoint de la taille normalisee et du nombre de 
voisins pour chaque fibre. La forme 3d de cet histogramme est donnee sur la figure 1.2 





Unc double fibre a 8 voisins 



Fig. 1.4 - Variation du nombre de voisins en fonction de la taille du segment detecte comme etant 
une fibre 

correction du biais : les images sont en general biaises, c'est a dire qu'il existe un fond lente- 
ment variable au lieu d'etre constant. 

Le bloc de soustraction du fond est optionnel, ainsi que celui de la correction du biais. 
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1.3.1 Estimation du flux vu par chaque fibre 

Le rapport [?] decrit une methode pour rattraper les problemes de saturation de l'image. Ces 
probleme vont disparaitre avec Parrivee de Pelectronique mi mai 2002. Nous avons done decide de 
ne pas continuer dans cette voie pour le moment. Ici nous avons simplement choisi d'estimer le 
flux avec la moyenne sur la fibre. 
Plusieurs ameliorations sont envisageables : 

- tenir compte du bruit poissonien dans Pestimateur du max avec la moyenne, 

- tenir compte de la disparite de forme de fibres du guide d'image ; on peut constater que les toutes 
les fibres ne sont pas identiques et que Pestimateur du flux devrait s'adapter a la forme exacte 
de la fibre, 

- verifier les problemes possibles dues a la zone inter-fibre, en fonction de la tache focale sur Pentree 
du guide, 

et. plus gwuiralcinerit re cuperer le rhodele coinplet tf 'injection et de ret our du spot decrit dans 



Ces ameliorations peuvent etre et.udiees, mais ne sont pas prioritaires pour le moment. 



1.3.2 Estimation ou utilisat ion du fond, puis soustraction 

Le fond peut avoir plusieurs sources; celles-ci sont decrite dans les memos [?, ?, ?]. En resume, 
nous appelons fond ici soit les reflexions parasites sur les optiques, et done y compris sur la sortie 
du guide d'images, mais egalement l'offset due a la chaine de numerisation. 

Si Poffset est dominant sur l'image, on ne peut pas obtenir le fond simplement en retirant Pimage, 
car Poffset depend du contenu, et n'est done plus le meme. Dans ce cas, on utilise un quantile de 
Phistogramme ou pour Pestimer. 

Dans le cas contraire, il faut utiliser le fond acquis lorsqu'on retire Pobjet a regarder, et le soustraire. 
Dans tous les cas il faut bien penser a saturer la soustraction pour ne pas etre gene par des outliers 
negatifs. L'utilisation d'un talon (cf. [?]) n'est pas retenue car le fond diffus n'est pas du tout une 
composante majoritaire du signal. 

1.3.3 Correction du biais 

Le biais doit a priori etre corrige sur l'acquisition des taux d'injection (branche de gauche) et sur 
Pacquisition de Pobjet (branche de droite). Pour le premier cas, cela vient du fait que la calibration 
se fait sur un miroir plan, et que la courbure de champ va reduire la qualite d'injection au retour 
sur les bords (qui sont defocalises) . Ce ne serait probablement. pas le cas 

Sur Pobjet, ou dans un milieu diffusant homogene, Pinjection reste moins bonne sur les bords, et 
cela se traduit par un biais tres similaire quant a sa forme au premier. 
Dans tous les cas le biais a une symetrie quasi circulaire. 

L'estimation du biais se fait en divisant Pimage en N x N blocs de taille fixe, puis en estimant le 
biais sur chaque bloc. Pour cela il faut considerer la nature de Pobjet observe. Dans le cas d'un 
objet homogene, le biais peut etre acquis en prenant. la valeur moyenne ou mediane sur le bloc. 
Quand il y a un objet, il faut savoir si cet objet est plus sombre ou plus claire que le biais. Dans 
notre cas, le biais est multiplicative, et on prend done plutot un operateur de moyenne ou de 
mediane (par rapport a un max ou min pour un biais additif). 

On obtient alors une image de taille N x N qui est. utilisee avec une interpolation pour trouver la 
valeur du biais vue par chaque fibre. On utilise pour le moment, une interpolation bilineaire, faute 
de mieux. La librairie inrimage pourrait etre completee pour integrer une interpolation en spline 
cubique point a point. 

On peut remarquer que le biais peut etre estimee plus finement avec d'autres methodes plus subfiles 
qu'il faudrait etudier. Notre methode reste cependant rapide et semble suffisante. 
Une fois le biais estime, il faut. diviser Pimage par son biais. II faut penser a verifier que le biais 
est bien superieur a un pour ne pas multiplier en fait ! Dans le cas de l'image de calibration, on va 
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retablir une dynamique Active sur 1000 niveaux afin de ne pas avoir de problemes avec la seconde 
division pour la calibration du taux d'injection (cf. partie 1.4.1). 

1.4 Bloc de calibration et reconstruction 
1.4.1 Calibration des taux d'injection 

Le bloc d'estimation du flux de gauche permet d'obtenir une image des taux de calibration fibre 
a fibre. On peut dont l'utiliser pour retablir une injection de 100% sur toutes les fibres. Pour cela 
on divise l'image de l'objet par l'image du taux d'injection, apres avoir pris quelques precautions 
sur les valeurs de la calibration (superieures a 1, avec suffisamment de dynamique). 



1.4.2 Reconstruction mosa'i'que 

La reconstruction mosai'que s'effectue en repartissant sur toute la surface de chaque fibre la valeur 
estimee du flux prise apres l'image de calibration. 

"Afin de donner un aspect plus "lisse,. on- a. rajoute .un....filtrage_recursif...passe baa de type Deriche. 
Le parametre [_outputSigmaFilter_] permet de regler l'ecart-type. Noter que la complexite de 
l'algorithme est independante de l'ecart type de la gaussiene. 



1.4.3 Reconstruction par RBF 

A faire... 

1.5 Resultats et discussion 

Les resultats sont visibles dans la base de donnees d'images... 

1.5.1 Exemple 

Les figures 1.5 et 1.6 illustrent toutes les etapes du traitement de l'image. Dans ce cas de figure, le 
fond de l'image de calibration des fibres est estime, ce pourquoi il n'apparatt pas sur les images. 
Le fond de l'objet est par contre fourni. Les images de calibration et l'objet sont toutes les deux 
de-biaisees. Ce resultat a ete obtenu avec les parametres par defaut. Sur cet exemple ; on remarque 
que les images des biais n'ont pas la symetrie circulaire habituelle. La raison est que la majorite 
du biais venait en fait de la saturation du detecteur qui se traduisait par un nouveau fond additif 
non prevu au depart, mais qui est tout de meme attenue par les traitements. 

1.5.2 Discussion 

De nombreuses idees d'ameliorations ont ete proposee dans ce chapitre. Certaines semblent plus 
importantes que d'autres. Void la proposition d'un tri de ces idees par ordre decroissant d'impor- 
tance. 

Reconstruction par base de fonctions radiales : cela permet d'obtenir une image lisse et de 

faire d'eventuels traitements ulterieurs. 
Supprimer le zoom sur l'image : le zoom est utile pour diminuer Peffet de l'ouverture nu- 

merique avant la LPE. II faut trouver l'operateur equivalent sur l'image non zoomee. Une 

hypothese a creuser est que ce serait equivalent a une ouverture avec un element structurant 

en croix ( ?) . 

Auto-contr6le de la detection des fibres : savoir si la detection est correcte, ou si au cours 
des acquisition on arrive a un moment ou il faut la refaire. 
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(e) Image de-biaisee issue de (d) 



(f) Resultat. final 



Fig. 1.6 - Etapes intermediaires de la procedure Image-Cell 
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Fig. 1.7 - Resultat final apres reconstruction (cetteimage a juste, ete lissee, et n'est pas une re- 
construction par RBF). 



Auto-contr6le de valeur du fond : savoir sir le fond propose est valable ou si il faut l'estimer. 

Dans le cas ou le fond n'est pas fournit, en faut-il un ? 
Amelioration continue de la detection des fibres : on devrait pouvoir utiliser les images 

des objets pour renforcer ou corriger la detection des fibres. 
Estimer le biais avec des splines : le biais est estime par reechantillonnage bilineaire qui pro- 

duit un effet de blocs qui peut etre genant 
On peut aussi prevoir plusieurs axes de recherche pour des fonctionnalites futures : ce sera le sujet 
d'un prochain memo... 
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Chapitre 2 

Programmation d'Image-Cell 



Introduction 



Cette partie rassemble les principaux elements de la librairie. La description de ces element est 
globale, cette documentation venant en complement de la documentation automatique 1 decrivant 
toutes les classes et toutes les fonctions, et des tests unitaires decrivant l'utilisation des classes et 
TFon^iras^principalesr 

2.1 Une librairie sur les graphes 

La programmation des algorithmes a necessite l'utilisation de structure de graphes. Nous avons 
done recherche des librairies C++ permettant de les manipuler. Nous avons teste les librairies 
suivantes. 

GTL : pour Graph Template Library (http: //www. infosun.f mi .uni-passau.de/GTL/). Cette 
librairie est la plus adaptee a nos besoins et a ete celle utilisee pour le prototypage dans 
[?]. Le probleme principal de cette librairie est son prix dissuasif pour les entreprises : 5 000 
Euros (licence de developpement pour un seul poste). Elle etait gratuite pour le stage car 
utilisee par une etudiante... Nous n'avons done pas choisi de garder cette librairie pour cette 
raison. 

LEDA : pour Library of Efficient Data Types and Algorithms (http://www.mpi-sb.mpg.de/LEDA/). 
Cette librairie traite de nombreux probleme de type de donnees, dont celui des graphes. Cette 
librairie a cependant beaucoup evoluee vers une sorte de librairie couteau Suisse qui fait beau- 
coup de choses, dont les interfaces graphiques ; et qui redefmit de nombreux outils maintenant 
standard, notamment les types STL. II existe done de nombreux cas de conflits potentiels 
avec nos softs. Le prix de cette librairie est de 1500 US$ pour un developpeur. et 8600$ pour 
developper sur tout un site, et vendre des produits derriere. J'ai personnellement utilise cette 
librairie pour ma these (gratuitement car pour une these...). 

BOOST : une sorte de groupware qui travail sur des proposition de futurs standard C++ 
(http://www.boost.org/index.htm). II y a en fait toute un librairie avec de nombreux 
themes. Le point central est que e'est de la programmation avancee, generique, reutilisable, 
versatile. C'est beau mais e'est tres complique a utiliser malgre une documentation tres 
propre. Le probleme vient des pre-requis indispensables en programmation objet de haute 
voltige qui sont un peu trop complexes pour nos besoins. et qui necessiterai un travail sup- 
plemental important. Cette librairie est cependant gratuite, et tres bien faite. Un attrait 
interessant est qu'elle n : est constitute que de fichiers -h, sans partie a compiler. Tout repose 

^tilisant doc++ 
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sur une gestion complexe de templates qui ne peuvent pas etre precompilee. II n'y a done pas 



de problemes de compilation sur differentes plateformes. 
Toutes ces librairies fonctionnent sous toutes les plateformes courantes. et avec tous les compilateurs 
courant. De plus elles sont maintenues et. stables. II existe de nombreuses autres petite libriaires a 
droite a gauche mais qui sont loin d'approcher la qualite de celles-ci. Nous n'avons cependant pas 
pu en retenir une celle, soit a cause du prix et des risques lies a Tutilisation d'une grosse librairies 
supplemental, soit a cause de la complexity d'utilisation trap grande. 

Nous avons done choisi de prendre nos propres structures de donnees, adaptes aux traitement 
pour plus de rapidite a la fois de developpement et d'execution. Ces structures de donnees sont 
cependant assez souples car elle fonctionnent sur des types et des algorithmes STL qui peuvent 
etre interchangeables assez rapidement... 



2.2 — Gestion generate des parametres d'algorithmes 

2.2.1 L'objet parametre 

L'objct parametre n 'existe" pffi 

moment a l'etat de concept dans le meme sens que la terminologie STL. 
Une classe parametre est de type suivant : 

class AlgoParam { 

public : 

AdjGraphParamO-C 
defaults () ; 

}; 

void defaults () { I / Cette fonction fixe les parametres par defaut 



/// Load parameters from a file. If any error occurs, default values are provided, 
void loadParamC const char * fname) ; 

/// Save parameters in a file, 
void saveParam( const char * fname); 

/// Output stream 

friend std::ostream & operator « (std: :ostream ft, AlgoParam &) ; 
/// Input stream 

friend std: : istream & operator » (std: :istream &, AlgoParam &) ; 
public: // parameters list here, with comments please ! 



Les fonctions loadParam et saveParam vont juste creer un not dans le fichier dont le nom est 
specifie. lis deleguent done le travaille aux flots ostream et istream. 



> 



[...] 



>; 
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Pour le moments les flots sont simples, et ne font, que ecrire le nom du parametre suivie de la valeur 
numerique, le tout sur une ligne. Les parametres sont encadres par 

ALGOPARAM_CONFIGURATION 
PARAMETERS 

et 

END_ALG0PARAM_C0NFIGURATION 

Ce mecanisme (surement perfectible...) permet de deleguer les entrees et les sorties aux flots qui 
peuvent done etre imbriques. Ainsi un objet parametre peut contenir un autre objet parametre 
sans avoir a recrire les entrees et les sorties de cet objet. 

-De-telles-classes de-parametres sont utilises pour ImageCell Fiber-DeteGtor-Ad^Gr-aph-et-EiberElow. 



2.2.2 Le futur des parametres 

I , a premiere chose a fains est de creer une.cTasse.virt.ueirelCc.!est-a-dire line interface) pour la classe 
de parametres. 

Enfin, dans un second temps est a Petude Putilisation d'XML pour interfacer ces parametres. 

2.3 La detection des fibres : FiberDetector et AdjGraph 

Voir le test unitaire test.f iberdetect pour les details d'utilisation. 

L'objet FiberDetector est un objet-algorithme modulaire pour detecter les fibres. II contient un 
objet specialise AdjGraph qui permet de faire l'etape de fusion des fibres qui est la plus compliquee. 

2.3.1 Fonctionnement detaille 

On retrouve les etapes suivantes dans la classe : 

void init(Inrimage *inrln) ; 
void preprocessingO ; 
void watershedO ; 
int splitBiggest (); 
int mergeSmallest () ; 

L'etape d'initialisation init permet de fixer l'image servant de base k la detection des fibres. Puis 
l'etape preprocessing fait es premiers traitements : zoom, ouverture numerique, inversion de 
l'image, et quelques calculs preliminaire, comme une structure de voisinage des fibres. Ces deux 
etapes permettent de preparer l'image pour les traitements de segmentation. 
La segmentation est realisee par une ligne de partage des eaux (LPE) dans la fonction watershed. 
L'algorithme pourrait egalement interfacer un autre algorithme... 

Les deux etapes suivantes sont deux etapes de post-processing de split and merge. splitBiggest 
va reperer les fibres trop grosses, puis les resegmenter. mergeSmallest permet de fusionner les 
fibres les plus petites. Cette derniere etape plus complexe et. plus critique est elle-meme encapsulee 
dans une classe algorithme AdjGraph expliquee dans la prochaine partie. 

Ces deux etapes contiennent une dizaine de parametres qui permettent d'ajuster la precision et la 
sensibilite de l'algorithme. La structure de l'algorithme permet d'eventuellement boucler ces etapes 
afin de faire converger l'algorithme vers une segmentation stable au sens de Pidempotence du split 
and merge, pour un jeu de parametre donne 2 . 

2 Nous n'avons pas de preuve de l'existence d'une telle convergence et il semble difficile a priori d'en obtenir. 
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2.3.2 L'etape de fusion avec AdjGraph 

Cette partie decrit l'etape de fusion. Cette etape est importante car les petites fibres (sur-segmentees) 
sont plus genantes que les grosses (sous-segmentees) sur un resultat final. La complexity de cette 
partie nos a oblige a en faire une partie modulaire, au meme titre que le FiberDetector. Les 
differentes «modules» ont ete les suivants : 

- type de calcul de la compacite d'une region, et de la fusion de deux regions 

- construction du graphe d'adjacence dedie 

- structures de donnees dediee aux calculs des attributs utiles. 

Cette etape etait la plus longue sur le prototype de [?], et a done fait l'objet d'optimisation. L'etape 
de fusion n'est realise que sur des criteres de formes, sans lien directe necessaire avec l'image de 
depart 3 . Le probleme principal est le calcul des attributs utilises : la compacite et la surface. La 
principale difficulte est liee a ce que le graphe va evoluer par fusion de regions, et que Ton cherche 
a mfettrfe &. jour c^ areritaflsrciP" vtrar^atenrent testCTies-fasioimTJOssibles^t-donc^r^ealetdb^e^ 
attributs possibles. 

Pour cela. on distingue les frontieres exterieures et interieures d'un objet comme indique sur la 
figure 2.1. On divise la front \C-rv d'un objei en foncrion de l'adjacence. La frontiere de A est e gale a 
la somme d(>s frontiere interieuroH A avec i.ous ses voisitus. La front itire interienro de A commune 
avec F est egale a la frontiere exterieure de F commune avec A. Cette derniere propriete permet 
de ne parler que de frontiere interieure. 




FlG. 2.1 - Les frontieres exterieurs et interieures d'un objet. La frontiere interieure de A commune 
avec F est egale a la frontiere exterieure de F commune avec A. Ces informations sont portes par 
des arcs du graphe d'adjacence. 

Ces informations sont portes par des arcs du graphe d'adjacence. Un arc allant de A vers B va 
porter l'information de frontiere commue que A detient avec B. 

On peut formaliser ces relations ainsi. Soit G = (AT, E) un graphe avec N l'ensemble de ses noeuds 
et E C N x N l'ensemble de ses arcs. G est un digraphe 4 ou deux nceuds n'ont soit pas d'arcs 
en commun, soit deux arcs antisymetriques. E est done un ensemble de paires ordonnees. Soit s : 
N — » R + une fonction permettant d'obtenir la surface. Soit / : E -> R + une fonction permettant 
d'obtenir la frontiere interieure entre les deux nee designe par Tare. Enfin soit c : N -» R + la 
fonction permettant d'obtenir la compacite d'une region. On note V(n), n e N l'ensemble des 
voisins de n. La compacite peut alors s'ecrire : 



3 Ce qui est en fait une extension (lourde !) possible 
4 un graphe ou tous les arcs sont diriges 
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Pour tester la compacite de la fusion eventuelle de n et m, on note la nouvelle compacite c(n U m) : 
, ^ (E oe y ( n) HM) + £ oe v(^) f(M ) - f((n,m)) - /((m,n))) 2 

Et si la fusion de m avec n est. effectuee, on modifie les fonctions (on suppose que m est rajoutee 
an): 

s(n)^s(n) + s(m) (2.3) 
f(n, o) «- /(n, o) pour o e V(n) et o ^ V(m) (2.4) 
/(n, o) «- /(n, o) + /(m, o) pour o e V(n) n V(m). (2.5) 

On peut ainwi ne travailkr quo >sur les graphcs. sans avoir a rcyonir a la segmentation initiate pour 

ster beaucoup i>1hk <le fusion. 



e eventuellement de tester toutes les fusions possibles. 



2.4 Structures de donnees associees aux fibres : FiberStruc- 
ture et Fiberlnfo 

Voir le test unitaire test_f iberlnf o pour les details d'utilisation. 

2.4.1 Lien entre la structure des fibres et les informations sur chaque 
fibre 

La structure de donnee est specifique au probleme qui nous interesse. Les informations qui nous 
interessent sont des informations attachees aux fibres seulement. les fibres etant en fait liees a la 
structure detectee. Les deux types d 'informations sont done etroitement liees, les informations sur 
les fibres ne pouvant exister independamment de la structure des fibres. 

Modele sujet-observateur 

Nous avons done choisi de les relier sous le modele de l'observateur des design patterns [?]. Nous 

avons done defini ee mecanisme dans le fichier Observer. h. Le principe est le suivant : 

«Definit une interdependance un a plusieurs, de facon telle que, quand un objet change 
d'etat, tous ceux qui en dependent, en soient notifies et automatiquement mis a jour.» 

Ainsi on distingue le sujet. qui pour nous est la structure des fibres detenue par la classe FiberStructure ; 

et les observateurs qui sont les informations que Ton peut attacher aux fibres (classe Fiber Inf o<T>). 

On a ajoute un comportement supplementaire : un objet Fiber Inf o<T> ne peut exister sans un 

FiberStructure, ce pourquoi il n'existe qu'un constructeur qui prend en argument la FiberStructure. 

L'objet Fiberlnf o<T> va d'ailleurs s'initialiser tout seul en regardant la structure a laquelle il est 

attache. 

Quand la structure vient a changer, tous les Fiberlnf o<T> en sont informes si besoin, et ceux-ci 
vont alors se reinitialiser automatiquement. Quand la structure est detruite, les objets d'information 
sont informes egalement qu'ils ne sont. plus valides. 



Structure codee comme un graphe : FiberStructure 

La classe fiberstructure se construit a partir de Fimage des composantes connexes. L'objet va alors 
garder en interne une copie de cette image de reference. Elle contient en interne la structure de 
graphe associee alors aux composantes connexes. Elle contient egalement d'autres informations, 
comme les baryeentres des regions, en coordonnees pixel, les surface des regions ... 
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Tout.es ces informations sont. accessibles depuis l'interface de la classe, mais ne sont pas modifiables. 
car dies dependent de la structure meme. 

La classe d'information : Fiberlnfo 

La classe Fiber Inf o<T> est une classe template car les informations peuvent etre de nature tres 
differentes. Elle contient un vecteur STL d'information dont la taille est geree automatiquement 
en fonction de la structure associee. Ce vecteur est un membre public afin d'avoir des acces rapides 
aux informations. La classe sert juste de decorateur 5 pour pouvoir gerer les dependances avec la 
structure des fibres. 

Quand le type du template est numerique. il existe des fonctions pour additionner, soustraire ; 
diviser et multiplier deux Fiber Inf o<T>, ou bien un Fiber Inf o<T> avec un element de T : addFI. 
subFI/divFL mulTT 



2.4.2 Transformation d'images en Fiberlnfo et reciproquement 

Voir le test unitaire test_f iberProcessing pour les details d 'utilisation. 

Une image peut etre transformee en Fiberlnf o<f loat> avec la fonction d'estimation simple des 
flux : 

int f lowEstimation (Inrimage* image, Fiberlnf o<float> & flow, bool bias = false). 
II faut au prealable initialiser de Fiberlnf o<float> flow avec la structure, ainsi l'algorithme 
connalt via flow l'image des composantes connexes a utiliser. 

La transformation reciproque est la reconstruction d'image. pour laquelle il n'existe pour le moment 
que la methode de mosa'fquage : 

int imageMosaic (Inrimage* image, Fiberlnf o<float> & flow) 

Le principe est similaire a f lowEstimation. La structure est passee via le flow. 

2.5 Soustraction du fond et biais : FiberFlow 

Voir le test unitaire test.f iberFlow pour les details d'utilisation. 

2.5.1 Fonctionnement general 

La classe FiberFlow implemente le bloc d'estimation des flux de la partie 1.1, sauf la partie 
de transformation de l'image en Fiberlnfo. La classe comporte deux methodes principales : 
setBackground et set Image, et son comportement est parametre avec la classe de parametre 
FiberFlowParam. La classe permet de soustraire le fond qui est donnee, ou bien il peut l'estimer. 
La classe sert aussi a calculer le biais dans l'image, apres soustraction du fond. 

2.5.2 Le calcul du fond 

Voir le test unitaire test_f iberHisto pour les details d'utilisation. 

Pour calculer le fond, on calcul un quantile de Phistogramme des fibres. Pour cela il existe une 
classe FiberHisto<T> qui prend un FiberInfo<T> en constructeur, et qui permet de faire des 
calculs sur Phistogramme calcule a partir de Pinformation des fibres. Ces calculs sont plus rapides 
que sur l'image elle-meme. 

5 At.tention, pas au sens du Design Pattern decorateur [?] ! 
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2.5.3 Le calcul du biais 



Voir le test unitaire test_f iberProcessing pour les details d'utilisation. 

Le calcul du biais n'est en fait pas calcule sur Timage, mais encore une fois sur l'information des 
fibres directement. C'est encore un cas ou on utilise la structure et rinformation des fibres. Ce 
calcul est fait dans une fonction : 

int biasEstimate ( Fiber Inf o<f loat> & flowln, 
Fiber Inf o<float> & bias, 
int nbBlocks, 

biasEstimatorType bet=BET_MEAN ) 

L'estimateur du biais est pour le moment a choisir em re la tuoyentie (BET.KEAN) et. la mediane 

(BET.KEDIAN). : " _ 

On divise l'espace total des fibres en nbBlocks blocs de taille egale. On parcourt ensuite Pensemble 
des blocs et on trie les fibres en fonctions du bloc dans lequel elles tombent. Puis, pour chaque 

-bloc,-on applique l'operateiu-. choisi. d!estimation-du..biais 

On calcule le biais ensuite en ihterpolanTpouT chaque fibre Vihiqtiement (et non pour toute l'image) 
les valeurs des blocs les plus proches. On utilise pour le moment une interpolation bilineaire. car 
c'est la seule disponible avec Inrimage. On fait une proposition pour avoir une interpolation spline 
cubique. 

2.6 Ameliorations 

- La classe de parametres est a optimiser et a programmer plus proprement. 

- Le biais devrait etre interpole avec une spline cubique 

- Une librairie de graphe serait tout de meme la bienvenue. Cela permettrait d'avoir une structure 
commune entre FiberDetector, Adj Graph, et FiberStructure. ce qui n : est actuellement pas 
exactement le cas. 
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Chapitre 3 

Utilisation d'Image-Cell 



C'est une partie delicate car elle doit contenir le savoir-faire acquis pendant le prototypage, le 

developpement et les nombreux tests. 

t petit! hi sur-s gmentee t> plu intef |u< le 51 3S m se mentees m 111 
resultat final. 

3.1 Utilisation de la parametrisation 

II faut ouvrir une fengtre d'un editeur de texte (notepad, emacs, wordpad...) avec le fichier ima- 
gecell.txt. Ce fichier est relu a chaque utilisation. La liste des parametre, ainsi que des resumes 
d 'indications sont donnees en annexe A. 

3.2 Parametrage de FiberDetector 

3.2.1 Pr^traitements 

Dans presque tous les cas, il faut faire le zoom sur l'image pour que l'algorithme fonctionne cor- 
rect ement. 

3.2.2 Parametrage de AdjGraph 

Les petites fibres (sur-segmentees) sont plus genantes qUe les grosses (sous-segmentees) sur un 
resultat final. 

3.3 Parametrage du FiberFlow pour la calibration 

3.4 Parametrage du FiberFlow pour l'image a reconstruire 

3.5 Parametrage d'ImageCell : les reste des parametres 



25 



Document confidentiel - Mauna Kea Technologies - 2002-03-12 



Annexe A 



Parametres de Image-Cell 



Les parametres internes d'lmage-Cell sont modifiables via un fichier de parametres present dans 
le repertoire source du lancement de Papplication, au meme endroit. que tomoscope.cfg. Le fichier 
des parametres imagecell . cf g a la forme suivante : 

27 





Ligne du fichier | Defaut | Possibles 


Remarques 




IMAGECELLPARAM CONFIGURATION 


Parametres de l'algorithme entier 




PARAMETERS 




FIBERDETECTORPARAM_CONFIGURATION 


Parametres de Palgorithme entier 


PARAMETERS 




doZoom 


1 


{0,1} 


Travail sur l'image doublee (plus precis) 


_doDiffusion_ 


0 


{0,1} 


A regler. 


_hDomeWatershed_ 


1 


{1,-,10} 


Hauteur des maxima selectionnes pour initia- 
liser les fibres. 


_oversizedFactor_ 


1.8 


[i;2] 


Taille normalisee (par rapport a la moyenne) 
au dela de laquelle une fibre va etre divisee. 
Co. crkero est utilise av(!(: le critere suivant sur 










la taille (ET-logkrae) 






_ over SizedNeighLimit _ 


6 


{5,6,7} 


Nombre de voisins au dela duquel une fibre va 
etre divisee. Ce critere est utilise .avee le critere 


















_overSizedResegmentMethod_ 


1 


{0,1} 


resegmenteTda^ 
traitement 




ADJGRAPHPARAM_CONFIGURAT 


ION 


Parametres pour la fusion des fibres 


PARAMETERS 




_underSizedFactor 


0.6 


[0.1; 1] 


en dega de laquelle une fibre va etre une candi- 
date pour la fusion. Ce critere est utilise avec 
le critere suivant sur la taille (ET logique) 


_underSizedMendatoryMergeSlope_ 


-0.22 


[?;?] 


Ce parametre et le suivant determinent 
une droite dans l'espace (nombre de voi- 
sins)x(taille normalisee) delimitant les fibres 
qui doivent obligatoirement etre fusionnees. 
Mieux vaut ne pas y toucher sans faire de cal- 

points (6; 1/3) et (3, 1). Augmenter la pente 
revient a s'occuper plus des grosses fibres avec 
peu de voisins. 


_underSizedMendatoryMergeOffset_ 


1.667 


[?;?] 


Voir remarque precedente. Ce parametre 
nees=tailles normalisees). L'augmenter re- 
fibres. ' 


_ underSizedFilterSize _ 




Ri} 


de la taille totale de la fusion ? 


_underSizedFilterSizeMax_ 


1.17741 


[i;3] 


de fibre normalisee superieure a 


_underSizedFilterCompact_ 


0 


{0,1} 


Filtrage des fibres fusionnees en fonction de la 
compacite 


_underSizedFilterCompactMax_ 


2 


[1;10] 


Facteur multiplicatif de l'ecart-type 
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Ligne du fichier | Defaut 


Possibles 


Remarques 


_underSizedFiIterMoreCompact_ 


0 


Ri} 


Ne garde que les fusions qui rendent les fibres 
fusionnees plus compacte que la fibre voisine 
seule. 


END ADJGRAPHPARAM CON 


FIGURA r 


noN 



END FIBERDETECTORPARAM CONFIGURATION 



FIBERFLOWPARAM_CONFIGURATION 



PARAMETERS 



_backgroundMode_ 



_biasCorrectionOperatoi 



{0,1,2} 



{0,1} 



{0,1} 



END_FIBERFLOWPARAM_CONFIGURATION 



Parametres d'estimation du flux pour l'image 
de calibration 



0 : pas de soustraction du fond ; 1 : soustrac- 
tion avec le fond fourni, sinon 2 ; 2 : estimation 
du fond sur rhistogramme 



Quantilo de V 
* 



Correction du biais. 
fond 



apres soustarction du 



a Y pour son estimation 



Operateur pour la correction du 1 
moyenne ; 1 : mediane 



FIBERFLOWPARAM_CONFIGURATION 



PARAMETERS 



_backgroundMode_ 



_backgroundEstimationCut _ 



biasCorrectionNbBlocks 



_biasCorrectionOperator_ 



{0,1,2} 



[0;0.1] 



{0,1} 



{2, ..,64} 



{0,1} 



END_FIBERFLOWPARAM_CONFIGURATION 



Parametres d'estimation du flux pour l'image 
a reconstruire 



0 : pas de soustraction du fond ; 1 : soustrac- 
tion avec le fond fourni, sinon 2 ; 2 : estimation 
du fond sur rhistogramme 



Quantile de Phistogramme de l'image pour 
l'estimation du fond 



Correction du biais, apres soustarction du 
fond 



Si correction du biais, nombre dp blocks 
et en Y pour son estimation 



Operateur pour la correction du 1 
moyenne, 1 : mediane 



_tempoFilterTypeFib_ 



_outputSigmaFilter_ 



{0,..,4} 



{0,1} 



[0;10]U{-1} 



END_IMAGECELLPARAM_CONFIGURATION 



Type du filtre temporel pour l'estimation des 
fibres :0 : moyenne ;1 : moyenne coupee; 2 : 
medianne : 3 : max ; 4 : min 



Effectuer la calibration 



Ecart type de la gaussiene pour un filtrage 
passe bas en sortie de image cell. Si la valeur 
est negative (-1), il n'y a pas de filtrage. 
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Sujet : Manuel de reference de la partie libMKTProcessing pour imageCell. Ce manuel contient 
les informations de traitement d'images, de programmation C++, et d'utilisation. 

Version : Revision : 1.3 Attention, cette version est la version du fichier principal. Chaque 
chapitre possede sa propre version car les chapitres sont independants les uns des autres. 

Reference du document : $Id: Ref erenceImageCell.tex, v 1.3 2002-04-24 07:14:48 aymeric Exp $ 

Introduction 

Ce document, rassemble les differentes informations de reference pour le module p rincipal intitule 
Jmage- Cell. Ce module contient les parties suivantes de traitement d'imag es : 

- detection des fibres. 

- estimation des flux revenant des fibres. 

-^-corrections des defauts de -l.!appateil-(biais,..fond,— ..)., 

- calibration de rappareil, 

- reconstruction d'images. 

Ces differentes parties seront abordees en detail dans chaque partie pour expliquer les algorithmes 
et les parametres (chapitre 1). Pimplantation de ceux-ci (2) et leur utilisation pratique (3). 

Corrections 

- rajout d'un filtrage passe bas, mars 2002. 
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Chapitre 1 

Traitement d'images dans Image-Cell 



Introduction 

Auteur CVS (du chapitre) : Author : aymeric 
Version (du chapitre) : Revision : 1.4 

Reference du document : $Id: chapTdi . tex , v 1.4 2002-04-24 07:14:48 aymeric Exp $ 

Cette partie decrit les algorithmes en commencant par la description generate, jusqu'aux boites 
elementaires. 

1.1 Schema global 

La figure 1.1 represente le schema global du traitement. On peut distinguer quatre blocs, dont 
deux sont identiques, a une parametrisation pres. II existe done trois groupes de traitements qui 
sont les suivants : 

Detection des fibres : ces traitements permettent de detecter et d'isoler chaque fibre sur une 
image, ainsi que d'analyse la structure d'agencement des fibres du guide d'images ; 

Estimation des flux : une image brute de taille 640 x 640 x 20 (largeur, hauteur, nombre 
d'images temporellement )contient 8 megapixels qui represented information vue par 10000 
ou 30000 fibres. Ce bloc permet d'isoler Finformation effectivement vue par chaque fibre ; 

Calibration et reconstruction : La connaissance de l'information vue par chaque fibre est 
ensuite traitee pour etre reconstruite sous la forme d'une image debarrassee des defauts de 
l'appareil. 

Chacun de ces trois blocs sont maintenant detailles. Le prototypage des algorithmes decrits ici est 
detaille dans le rapport de stage de Sandra Marti [?] ; nous renvoyons le lecteur a cette reference 
pour plus de renseignements sur la demarche de developpement des algorithmes, notamment celui 
de detection des fibres. 

Le bloc de gauche d'estimation des flux permet de mesurer le flux sur une image qui represente un 
objet aux proprietes constantes dans l'espace (un milieu diffusant. homogene ou un objet presque 
homogene en mouvement aleatoire, un miroir). La sortie de ce bloc est. done une image des taux 
d'injection dans chaque fibre, et qui en prend en compte Pintegralite de la chaine. 
Le bloc de droite d'estimation des flux permet de Pestimer pour un objet a observer, et a travers 
l'appareil. Le dernier bloc de calibration est une division de l'image des taux d'injection par Pimage 
de l'objet observe. Cette operation permet de compenser les mauvaises injections dans certaines 
fibres. Puis Pimage est reconstruite. 
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Fig. 1.1 - Schema global du traitement. Les cadres gras representent les entrees. Les cadres poin- 
tilles representent les donnees ou les traitements facultatifs (biais, soustraction du fond...). 

1.2 Bloc de detection des fibres 

La detection des fibres s'organise autour de quatre traitements successifs : 

- pretraitements, 

- ligne de partage des eaux (LPE) : la premiere segmentation par region, 

- separation des fibres trop grosses (split), 

- fusion des fibres trop petites (merge). 

Les deux dernieres etapes peuvent etre bouclees: neanmoins une analyse de performance et de 
convergence doit alors etre menee. Nous n'aborderons pas ce sujet pour le moment. 

1.2.1 Pretraitements 

Les pretraitements realises sont les suivants. 

Diffusion anisotrope : le but est. de lisser Fimage dans les zones plates, c'est a dire les zones 
inter fibres. Le prototype fonctionnait mais n'utilisait. pas la librairie inrimage. mais une autre 
libriairie (celle utilisee par Sandra, issue de l'ENST). Le reglage de la diffusion est laissee a 
plus tard, lors de la prochaine release de la liblnrimage qui contiendra de la documentation 
a ce sujet ! Parametre [_doDif f usion_] . 

Interpolation x2 au plus proche voisin : on cherche a simuler des elements structurant de 
morphologie mathematique avec un rayon inferieur a un. L'image est doublee pour que l'ou- 
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verture qui suit ne touche pas aux maxima isoles, mais seulement ceux qui sont 8-connexe, 

mais non 4 connexe (voisins par une diagonale). L'interet est. de faire une selection des maxima 

elimine par Pouverture. Parametre LdoZoomJ. 
Ouverture numerique : on cherche a eliminer les maxima parasites situes sur les fibres. C'est 

un pretraitement classique de la LPE. mais eventuellement modifie par le comportement du 

doublage de 1 'image explique precedemment. 
Inversion de Pimage : la LPE fonctionne a partir de minima, on inverse done Pimage pour 

transformer les minima en maxima. 
L'operation d 'interpolation et d'ouverture qui suit devrait pouvoir etre remplace par un nouvel 
operateur de morphologie mathematique qu'il faudrait definir. Cela permettrait un gain important, 
puisque toutes les operations suivantes se feraient sur une image de taille inferieure. 



1.2 .2 — P remier w atershed 

Le watershed de la liblnrimage necessite la definition de marqueurs dans Pimages. Nous obtenons les 
-marqueursxcmmejmruma-deJJmag 

sont dermis pour un entier h, une fonction / : 2Z x 2Z -» ZZ representant Pimage concernee : 

g = EfU + h)-f, (1-1) 

avec E£°(b) la reconstruction geodesique par erosion de b dans a. g designe ici Pensemble des h- 
domes, que Ton seuil a 1 afin de tous les recuperer. On obtient ainsi les minima de profondeur au 
moins egale a h. Ce parametre h est donne par [_hDomeWatershed_] . 

Le resultat de cette LPE est une image de composantes connexes representant chaque fibre detectee. 
Afin de resoudre bon nombre de problemes de bord on colle les composantes touchant le bord au 
bord lui-meme. Ces composantes sont considerees comme faisant partie du bord, et du fond et sont 
mise a la composante connexe d : indice 0, qui est Pindice reserve au fond, et au bord. 
Enfin, les composantes connexes sont reorganises pour ne pas presenter de trous. et sont tries par 
ordre decroissant de taille. 

1.2.3 Split 

Les pretraitements, et notamment 1'ouverture numerique servent a limiter le defaut principal de la 
LPE qui est de sursegmenter. En faisant, cela, on a tendance a augmenter le nombre de segments 
qui vont englober plusieurs objets (ici, des fibres). Le premier defaut que Pon corrige est la sous- 
segmentation inevitable au traitement, mais surtout rajoutee par ces pretraitements. Pour cela on 
effectue les traitements suivants. 

Selection des segments suspectes de sous-segmentation : on trouve les fibres dont la taille 
est superieure a un seuil fixe en fonction de la moyenne des tailles normalisees (parametre 
[_oversizedFactor_] ), et en fonction du nombre de voisins (parametre LoverSizedNeighLimitJ ). 

Re-segmentation de ceux-ci : Les segments selectionnes sont isolees et re-segmenter soit sur 
Pimage initiale sans pretraitement, soit sur Pimage de carte de distance a Pinterieur de ces 
regions (parametre [_overSizedResegmentMethod_] ). Dans le premier cas, on suppose que 
les defauts majoritaires proviennent des pretraitement, dans le second cas on suppose qu'ils 
etaient la avant, et que Pon souhaite separer les segment aux endroits de retrecissement 
(forme de la cacahuete). La LPE est effectuee avec les meme parametres que la premiere. 

1.2.4 Merge 

Puis on cherche a corriger le defaut le plus classique (pour une LPE) : la sursegmentation. On 
cherche done a fusionner les segments. Pour cela on effectue les operations suivantes : 
- preselection des candidats a la fusion, 
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- parmi ceux-ci trier ceux qui seront. obligatoirement fusionnes, 

- pour les autres, eliminer les fusions impossibles 

- repertorier les fusions possibles 

- eliminer les fusions donnant de mauvais resultats 

- si il ne reste pas de fusion apres filtrage, retirer la fibres des fibres a fusionner 

- pour toutes les fibres restant.es, prendre la meilleur fusion au sens de la compacite des segments 
fusionnes. 

Les points importants sont detailles ci-apres. 
Remarques sur la compacite 

La compacite utilisee est le crite re le plus simple existant. Pour un objel de perimetre P et de 

surfaxie^SUfi-oiiare est ^, L e p er imetre est p ri s comme etant la frontiere de l'objet 4-connexe^ et 

done est defini comme tous les pixels de l'objet ayant un voisin 8-cormexe avec le fond (defini ici 
comme etant tout le reste de l'image. La surface est simplement l'objet lui-meme. 
Gependant^ce-critere simple a ealculer-n^st-pas-invariaHt-avee-Peehelle de l'objet; On peut-montrer — - 



que pour des formes simple, le critere est croissant avec la faille de l'objet au lieu d'etre invariant. 
Cette croissante n'est cependant qu'en moyenne en avec n la diametre de l'objet (ici, pour une 
boule en 4-connexite). Ce critere reste done tout de meme assez raisonnable, mais plus delicat a 
manier dans le cas de la comparaison de la compacite d'objets de taille differentes. 
II faut remarquer ici que les premiers essais sont realises sur une image deformee par les distortions, 
et qu'il faudrait, pour etre plus precis soit corriger l'image avant, soit tenir compte de l'anisotropie 
non stationnaire spatialement : on a tendance a surestimer a la fois les frontieres et la surface sur 
les bords de l'image. 

Selection des fibres detectees trop petites 

La selection s'effectue en deux etapes. On commence par trouver les segments candidats a la 
fusion avec un seuil sur la taille. On utilise le parametre LimderSizedFactorJ qui est un facteur 
multiplicatif de la moyenne, et est done inferieur k 1 en general. La seconde etape est de trouver 
parmi ces segments lesquels doivent obligatoirement etre fusionnes. et lesquels doivent etre filtres. 

Les segments obligatoirement fusionnes 

La selection des segments (ou fibres) obligatoirement fusionnes s'effectue dans l'espace (nombre de 
voisins) x (taille normalisee). La taille est. normalisee par rapport a la moyenne (on divise toutes 
les tailles par la moyenne des tailles). Le nombre de voisins est compte avec la 8-connexite. 
Les figures 1.2 et. 1.3 represented l'histogramme conjoint de la taille normalisee et. du nombre de 
voisins pour chaque fibre. On remarque qu'il existe un axe principal qui est centre sur la droite 
passant par des points d'interets particuliers (6,1), (8,2), et (5, .5). Ces points peuvent, s'explique 
geometriquement, comme indique sur la figure 1.4 

Cet axe est done un axe naturel de la structure hexagonale des fibres dans cet. espace. En s'orien- 
tant. sur cet axe, les petites fibres se situent vers un nombre de voisins faibles et une taille 
faible. On determine une frontiere lineaire pour selectionner les fibres detectees reellement trop 
petites. Apres quelques tests sur des images reelles de la base de donnees «Novembre 2001» et. 
en fonction des resultats preliminaires de [?], on a choisi une droite passant par les points (6: |) 
(3; 1). Les parametres de cette droites sont donnees sous la forme y = mx + p, avec m la pente 
[_underSizedMendatoryMergeSlope_] et p l'ordonnee a l'origine LunderSizedMendatoryMergeOff set_] . 
Tout point en dessous de la droite sera fusionne obligatoirement, et ne va done pas passer par les 
filtres decrits dans la section suivante. Les autres fibres sont par centre filtree comme suit. 
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Effectif 




Fig. 1.2 - Vue 3d de Phistogramme 2d conjoint de la taille normalisee et du nombre de voisins 
pour chaque fibre. Des lignes de niveaux sont projetees sur la fond du graphique : les lignes sont 
centrees autour de 6 voisins et de la moyenne de la taille des fibres. 



Filtrage des fusions possibles 

On entend par filtrage ici un precede qui retire certaines fusions. II existe fcrois filtres possibles. 

Filtrage sur la taille de la fusion : (parametre [_underSizedFilterSize_] ) on elimine les 

fusions dont la taille totale depasse un seuil fixe par le parametre [_underSizedFilterSizeMax_] , 
qui est pris comme facteur multiplicatif de l'ecart type ajoute a la taille moyenne. 

Filtrage sur la compacite totale : (parametre [_underSizedFilterCompact_] ) on elimine les 

fusions dont la compacite totale depasse un seuil fixe par le parametre [_underSizedFilterCompactMax_] , 
qui est pris comme facteur multiplicatif de l'ecart type ajoute a la taille moyenne. 

Filtrage sur le changement de compacite : (parametre LunderSizedFilterMoreCompact_] ) 
on elimine les fusions qui augmentent la compacite de Fobjet avec lequel on veut fusionner. 

Dans les parametres par defaut. on n'utilise pour le moment que le premier filtre qui semble suffisant 

dans la plupart des cas. Les autres filtres semblent trop restrictifs et. doivent etre encore etudies. 

1.3 Bloc d'estimation des flux 

Le bloc d'estimation des flux comprend plusieurs sous-parties : 

estimation du flux vu par chaque fibre : c'est ce que Ton peu recuperer d'une image a l'aide 

de la detections es fibres ; 
estimation du fond de l'image : ici le fond designe les reflexions parasites ou (ou inclusif !) 

1 'offset du a Telectronique et au detecteur; 
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1.5 

Taille normalisee 
















.5 













Nombre de voisins 



1.4e+03 600 80 

1.2e+03 400 40 

800 120 



Fig. 1.3 - Lignes de niveau de l'histogramme 2d conjoint de la taille normalisee et du nombre de 
voisins pour chaque fibre. La forme 3d de cet histogramme est donnee sur la figure 1.2 




Fig. 1.4 - Variation du nombre de voisins en fonction de la taille du segment detecte comme etant 
une fibre 

soustraction du fond : on retire le fond a l'image ; 

correction du biais : les images sont en general biaises, c'est. a dire qu'il existe un fond lente- 
ment variable au lieu d'etre constant. 



10 



Document confidentiel - Mauna Kea Technologies - 2002-04-24 



Le bloc de soustraction du fond est optionnel, ainsi que celui de la correction du biais. 



1.3.1 Estimation du flux vu par chaque fibre 

Le rapport [?] decrit une methode pour rattraper les problemes de saturation de l'image. Ces 
probleme vont disparaitre avec l'arrivee de l'electronique mi mai 2002. Nous avons done decide de 
ne pas continuer dans cette voie pour le moment. Ici nous avons simplement choisi d'estimer le 
flux avec la moyenne sur la fibre. 
Plusieurs ameliorations sont. envisageables : 

- tenir compte du bruit poissonien dans l'estimateur du max avec la moyenne, 

- tenir compte de la disparite de forme de fibres du guide d'image ; on peut constater que les toutes 
les fibres ne sont. pas identiques et que l'estimateur du flux devrait s'adapter a la forme exacte 
de~Ia~ftbre; 

- verifier les problemes possibles dues a la zone inter-nbre, en tonction de la tacKefocale sur l'entree 
du guide, 

- et p lus generalement recuperer le modele c omp let. d'injection et de retour du spot dec rit dans 



Ces ameliorations peuvent etre etudiees, mais ne sont pas prioritaires pour le moment. 

1.3.2 Estimation ou utilisation du fond, puis soustraction 

Le fond peut avoir plusieurs sources; celles-ci sont decrite dans les memos [?, ?, ?]. En resume, 
nous appelons fond ici soit les reflexions parasites sur les optiques, et done y compris sur la sortie 
du guide d'images, mais egalement l'offset due a la chalne de numerisation. 

Si l'offset est dominant sur l'image, on ne peut pas obtenir le fond simplement en retirant l'image, 
car l'offset depend du contenu. et n'est done plus le meme. Dans ce cas, on utilise un quantile de 
Phistogramme (ygg ou pour l'estimer. 

Dans le cas contraire, il faut utiliser le fond acquis lorsqu'on retire l'objet a. regarder, et le soustraire. 
Dans tous les cas il faut bien penser a saturer la soustraction pour ne pas 6tre gene par des outliers 
negatifs. L'utilisation d'un talon (cf. [?]) n'est pas retenue car le fond diffus n'est pas du tout une 
composante majoritaire du signal. 

1.3.3 Correction du biais 

Le biais doit a priori gtre corrige sur l'acquisition des taux d'injection (branche de gauche) et sur 
l'acquisition de l'objet (branche de droite). Pour le premier cas, cela vient du fait que la calibration 
se fait sur un miroir plan, et que la courbure de champ va reduire la qualite d'injection au retour 
sur les bords (qui sont defocalises) . Ce ne serait probablement pas le cas 

Sur l'objet, ou dans un milieu diffusant homogene, 1'injection reste moins bonne sur les bords, et 
cela se traduit par un biais tres similaire quant a sa forme au premier. 
Dans tous les cas le biais a une symetrie quasi circulaire. 

L'estimation du biais se fait en divisant l'image en N x N blocs de taille fixe, puis en estimant. le 
biais sur chaque bloc. Pour cela il faut considerer la nature de l'objet observe. Dans le cas d'un 
objet homogene, le biais peut etre acquis en prenant la valeur moyenne ou mediane sur le bloc. 
Quand il y a un objet, il faut savoir si cet objet est plus sombre ou plus claire que le biais. Dans 
notre cas, le biais est multiplicative, et on prend done plutot un operateur de moyenne ou de 
mediane (par rapport a un max ou min pour un biais additif). 

On obtient alors une image de taille N x N qui est utilises avec une interpolation pour trouver la 
valeur du biais vue par chaque fibre. On utilise pour le moment une interpolation bilineaire, faute 
de mieux. La librairie inrimage pourrait etre completee pour integrer une interpolation en spline 
cubique point, a point. 
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On peut remarquer que le biais peut etre estimee plus finement. avec d'autres methodes plus subfiles 
qu'il faudrait etudier. Notre methode reste cependant rapide et semble suffisante. 
Une fois le biais estime, il faut. diviser l'image par son biais. II faut penser a verifier que le biais 
est bien superieur a un pour ne pas multiplier en fait ! Dans le cas de Pimage de calibration, on va 
retablir une dynamique Active sur 1000 niveaux afm de ne pas avoir de problemes avec la seconde 
division pour la calibration du taux d'injection (cf. partie 1.4.1). 

1.4 Bloc de calibration et reconstruction 

1.4.1 Calibration des taux d'injection 

LeJ3lo.c_dlesiimatioii_di^ fibre 

a-fibre . On p e w t don t l' u t ilise r p o u r re t ablir u ne-irjjection^e-U)Q%-Sui^^ Pour cela 

on divise I'image de l'objet par l'image du taux d'injection, apres avoir pris quelques precautions 
sur Ies valeurs de la calibration (superieures a 1, avec suffisamment de dynamique). 



1.4.2 Reconstruction mosaique 

La reconstruction mosaique s'effectue en repartissant sur toute la surface de chaque fibre la valeur 
estimee du flux prise apres l'image de calibration. 

Afm de donner un aspect plus lisse, on a rajoute un filtrage recursif passe bas de type Deriche. 
Le parametre [_outputSigmaFilter_] permet de regler l'ecart-type. Noter que la complexite de 
l'algorithme est independante de l'ecart type de la gaussiene. 

1.4.3 Reconstruction par RBF 

A faire... 



1.5 Resultats et discussion 

Les resultats sont visibles dans la base de donnees d'images... 

1.5.1 Exemple 

Les figures 1.5 et 1.6 illustrent toutes les etapes du traitement de l'image. Dans ce cas de figure, le 
fond de l'image de calibration des fibres est estime, ce pourquoi il n'apparait pas sur les images. 
Le fond de l'objet est par contre fourni. Les images de calibration et l'objet sont toutes les deux 
de-biaisees. Ce resultat a ete obtenu avec les parametres par defaut. Sur cet exemple ; on remarque 
que les images des biais n'ont pas la symetrie circulaire habituelle. La raison est que la majorite 
du biais venait. en fait de la saturation du detecteur qui se traduisait par un nouveau fond additif 
non prevu au depart, mais qui est tout de meme attenue par les traitements. 

1.5.2 Discussion 

De nombreuses idees d 'ameliorations ont ete proposee dans ce chapitre. Certaines semblent plus 
importantes que d'autres. Void la proposition d'un tri de ces idees par ordre decroissant d'impor- 
tance. 

Reconstruction par base de fonctions radiales : cela permet d'obtenir une image lisse et de 
faire d'eventuels traitements ulterieurs. 
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Fig. 1.7 - Resultat final apres reconstruction (cetteimage a juste ete lissee, et n'est pas une re- 
construction par RBF). 



Supprimer le zoom sur l'image : le zoom est utile pour diminuer Peffet de l'ouverture nu- 
merique avant la LPE. II faut trouver l'operateur equivalent sur l'image non zoomee. Une 
hypothese a creuser est que ce serait equivalent a une ouverture avec un element structurant 
en croix ( ?). 

Auto-contr6le de la detection des fibres : savoir si la detection est correcte, ou si au cours 

des acquisition on arrive a un moment ou il faut la refaire. 
Auto-controle de valeur du fond : savoir sir le fond propose est valable ou si il faut l'estimer. 

Dans le cas ou le fond n'est pas fournit. en faut-il un ? 
Amelioration continue de la detection des fibres : on devrait pouvoir utiliser les images 

des objets pour renforcer ou corriger la detection des fibres. 
Estimer le biais avec des splines : le biais est estime par reechantillonnage bilineaire qui pro- 

duit un effet de blocs qui peut etre genant 
On peut aussi prevoir plusieurs axes de recherche pour des fonctionnalites futures : ce sera le sujet 
d'un prochain memo... 
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Chapitre 2 

Programmation d'Image-Cell 



Introduction - _. _ 

Auteur CVS (du chapitre) : Author : aymeric 
Version (du chapitre) : Revision : 1.2 

Reference du document : $Id: chaplnfo.tex.v 1.2 2002-04-24 07:14:48 aymeric Exp $ 

Cette paxtie rassemble les principaux elements de la librairie. La description de ces element est 
globale, cette documentation venant en complement de la documentation automatique 1 decrivant 
toutes les classes et toutes les fonctions, et des tests unitaires decrivant Putilisation des classes et 
fonctions principales. 

2.1 Une librairie sur les graphes 

La programmation des algorithmes a necessite Putilisation de structure de graphes. Nous avons 
done recherche des librairies C++ permettant de les manipuler. Nous avons teste les librairies 
suivantes. 

GTL : pour Graph Template Library (http://www.infosun.fmi.uni-passau.de/GTL/). Cette 
librairie est la plus adaptee a nos besoins et a ete celle utilisee pour le prototypage dans 
[?]. Le probleme principal de cette librairie est son prix dissuasif pour les entreprises : 5 000 
Euros (licence de developpement pour un seul poste). Elle etait gratuite pour le stage car 
utilisee par une etudiante... Nous n 'avons done pas choisi de garder cette librairie pour cette 
raison. 

LEDA : pour Library of Efficient Data Types and Algorithms (http://www.mpi-sb.mpg.de/LEDA/). 
Cette librairie traite de nombreux probleme de type de donnees, dont celui des graphes. Cette 
librairie a cependant beaucoup evoluee vers une sorte de librairie couteau Suisse qui fait beau- 
coup de choses, dont les interfaces graphiques.. et qui redefinit de nombreux outils maintenant 
standard, notamment les types STL. II existe done de nombreux cas de conflits potentiels 
avec nos softs. Le prix de cette librairie est de 1500 US$ pour un developpeur. et 8600$ pour 
developper sur tout un site, et vendre des produits derriere. J'ai personnellement utilise cette 
librairie pour ma these (gratuitement. car pour une these...). 

BOOST : une sorte de groupware qui travail sur des proposition de futurs standard C--I- 
(http://www.boost.org/index.htm). II y a en fait toute un librairie avec de nombreux 
themes. Le point central est que e'est de la programmation avancee, generique, reutilisable, 
versatile. C'est beau mais e'est tres complique a utiliser malgre une documentation tres 

1 Utilisant doc++ 
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propre. Le probleme vient des pre-requis indispensables en programmation objet de haute 
voltige qui sont un peu trop complexes pour nos besoins, et qui necessiterai un travail sup- 
plemental important. Cette librairie est cependant gratuite, et tres bien faite. Un attrait 
interessant est qu'elle n'est constitute que de fichiers .h, sans partie a compiler. Tout repose 
sur une gestion complexe de templates qui ne peuvent pas etre precompilee. II n'y a done pas 
de problemes de compilation sur differentes plateformes. 
Toutes ces librairies fonctionnent sous toutes les plateformes courantes, et avec tous les compilateurs 
courant. De plus elles sont maintenues et stables. II existe de nombreuses autres petite libriaires a 
droite a gauche mais qui sont loin d'approcher la qualite de celles-ci. Nous n'avons cependant pas 
pu en retenir une celle, soit a cause du prix et des risques lies a 1 'utilisation d'une grosse librairies 
supplemental, soit a cause de la complexity d'utilisation trop grande. 

Nous avons done choisi de prendre nos propres structures de donnees, adaptes aux traitement 
pour plus de rapidite a la""fois"d"e"developpem^ 

cependant assez souples car elle fonctionnent sur des types et des algonthmes STL qui peuvent 
etre interchangeables assez rapidement... 



2.2 Gestion generate des parametres d'algorithmes 
2.2.1 L'objet parametre 

L'objet parametre n'existe pas encore en tant que tel. mais va bientot voir le jour... II est pour le 
moment a Fetat de concept dans le meme sens que la terminologie STL. 
Une classe parametre est de type suivant : 

class AlgoParam { 
public : 

Ad j GraphPar am ( ) { 
defaults () ; 

}; 

void defaults () { // Cette fonction fixe les parametres par defaut 
> 

/// Load parameters from a file. If any error occurs, default values are provided, 
void loadParam( const char * f name) ; 

/// Save parameters in a file, 
void saveParam(const char * f name) ; 

/// Output stream 

friend std::ostream & operator « (std: :ostream &, AlgoParam &) ; 
/// Input stream 

friend std::istream & operator » (std: : istream &, AlgoParam &) ; 
public: // parameters list here, with comments please ! 
[..-] 
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}; 

Les fonctions loadParam et saveParam vont juste creer un flot. dans le fichier dont le nom est 
specifie. lis deleguent done le travaille aux flots ostream et istream. 

Pour le moments les flots sont simples, et ne font que ecrire le nom du parametre suivie de la valeur 
numerique, le tout sur une ligne. Les parametres sont encadres par 

ALGOPARAM_CONFIGURATION 
PARAMETERS 

et 

END_ALGOPARAM_C0NFIGURATION 
Ce mecanisme (surement perfectible...) permet de deleguer les entrees et les sorties aux flots qui 
^euventrlontretreTmbriques:^insiTm^ 



crire les entrees et les softies de cet objet. 
De telles classes de parametres sont utilises pour ImageCell FiberDetector Adj Graph et FiberFlow. 



2.2.2 Le futur des parametres 

La premiere chose a faire est de creer une classe virtuelle (e'est-a-dire une interface) pour la classe 
de parametres. 

Enfin, dans un second temps est a l'etude Futilisation d'XML pour interfacer ces parametres. 



2.3 La detection des fibres : FiberDetector et AdjGraph 

Voir le test unitaire test_f iberdetect pour les details d 'utilisation. 

L'objet FiberDetector est un objet-algorithme modulaire pour detecter les fibres. II contient un 
objet specialise AdjGraph qui permet de faire Fetape de fusion des fibres qui est la plus compliquee. 

2.3.1 Fonctionnement detaille 

On retrouve les etapes suivantes dans la classe : 

void init(Inr image *inrln) ; 

void preprocessing () ; 

void water shed () ; 

int splitBiggest () ; 

int mergeSmallest () ; 

L'etape d'initialisation init permet de fixer Fimage servant, de base a la detection des fibres. Puis 
Fetape preprocessing fait es premiers traitements : zoom, ouverture numerique, inversion de 
Fimage, et quelques calculs preliminaire, comme une structure de voisinage des fibres. Ces deux 
etapes permettent de preparer Fimage pour les traitements de segmentation. 
La segmentation est realisee par une ligne de partage des eaux (LPE) dans la fonction watershed. 
L'algorithme pourrait egalement interfacer un autre algorithme... 

Les deux etapes suivantes sont deux etapes de post-processing de split and merge. splitBiggest 
va reperer les fibres trop grosses, puis les resegmenter. mergeSmallest permet de fusionner les 
fibres les plus petites. Cette derniere etape plus complexe et plus critique est elle-meme encapsulee 
dans une classe algorithme AdjGraph expliquee dans la prochaine partie. 

Ces deux etapes contiennent une dizaine de parametres qui permettent d'ajuster la precision et. la 
sensibilite de l'algorithme. La structure de l'algorithme permet d'eventuellement boucler ces etapes 
afin de faire converger l'algorithme vers une segmentation stable au sens de Fidempotence du split 
and merge., pour un jeu de parametre donne 2 . 

2 Nous n'avons pas de preuve de 1'existencc d'unc telle convergence et il semble difficile a priori d'en obtenir. 
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2.3.2 L'etape de fusion avec AdjGraph 



Cette partie decrit. l'etape de fusion. Cette etape est importante car les petites fibres (sur-segmentees) 
sont plus genantes que les grosses (sous-segmentees) sur un resultat final. La complexity de cette 
partie nos a oblige a en faire une partie modulaire. au meme titre que le FiberDetector. Les 
differentes «modules» ont ete les suivants : 

- type de calcul de la compacite d'une region, et de la fusion de deux regions 

- construction du graphe d'adjacence dedie 

- structures de donnees dediee aux calculs des attributs utiles. 

Cette etape etait la plus longue sur le prototype de [?], et a done fait l'objet d'optimisation. L'etape 
de fusion n'est realise que sur des criteres de formes, sans lien directe necessaire avec l'image de 
depart 3 . Le probleme principal est le calcul des attributs utilises : la compacite et la surface. La 
principale difficulty est liee a ce que le graphe va evoluer par fusion de regions, et que 1'on cherche 
a mettre a jour ces attributs ; on veut egalement tester les fusions possibles et done pre-calculer les 
attributs possibles. 

Pour cela, on distingue les frontieres exterieures et interieures d'un objet comme indique sur la 

figureJLl.JD.rLdivise.lato A est .egale a ... 

la somme des frontiere interieures de A avec tous ses voisins. La frontiere interieure de A commune 
avec F est egale a la frontiere exterieure de F commune avec A. Cette derniere propriete permet 
de ne parler que de frontiere interieure. 




FlG. 2.1 - Les frontieres exterieurs et. interieures d'un objet. La frontiere interieure de A commune 
avec F est egale a la frontiere exterieure de F commune avec A. Ces informations sont portes par 
des arcs du graphe d'adjacence. 

Ces informations sont portes par des arcs du graphe d'adjacence. Un arc allant de A vers B va 
porter 1'information de frontiere commue que A detient avec B. 

On peut formaliser ces relations ainsi. Soit G = (N, E) un graphe avec N l'ensemble de ses noeuds 
et E C N x N l'ensemble de ses arcs. G est un digraphe 4 ou deux nceuds n'ont soit pas d'arcs 
en commun, soit deux arcs antisymetriques. E est. done un ensemble de paires ordonnees. Soit s : 
N -* R + une fonction permettant d'obtenir la surface. Soit / : E — > R + une fonction permettant 
d'obtenir la frontiere interieure entre les deux nee designe par l'arc. Enfin soit c : N — > R + la 
fonction permettant d'obtenir la compacite d'une region. On note V(n), n e N l'ensemble des 
voisins de n. La compacite peut alors s'ecrire : 

(EW<(^))) 2 

s(n) 



3 Ce qui eat en fe 
4 un graphe ou t( 
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Pour tester la compacite de la fusion eventuelle de n et m. on note la nouvelle compacite c(n U m) 
(E oeV (n) /((".«)) + E„ 6 y W /(K")) - f((n,m)) - f((m,n))) ? 



c(nUm) s{n) + s(m) 

Et si la fusion de m avec n est effectuee, on modifie les fonctions (on suppose que n 



(2-2) 
rajoutee 



f(n,o) <— /(n,o) pour o G V(n) et o ^ V(m) 
/(n, o) «- /(n, o) + /(m, o) pour o 6 V(n) n V(m). 



(2.3) 
(2.4) 
(2.5) 



— Orrpe ut ainsi ne Lra^aiiler~qire"SUT-les~graprfes7-saTrs a v o ir ' 
recalculer les attributs lors d'une fusion. Cela nous a permis d 



e eventuellement de tester toutes les fusions possibles 



2.4 Structures de donnees associees aux fibres : FiberStruc- 
ture et Fiberlnfo 

Voir le test unitaire test_f iberlnf o pour les details d'utilisation. 

2.4.1 Lien entre la structure des fibres et les informations sur chaque 
fibre 

La structure de donnee est specifique au probleme qui nous interesse. Les informations qui nous 
interessent sont des informations attachees aux fibres seulement. les fibres etant en fait liees a la 
structure detectee. Les deux types d 'informations sont done etroitement liees, les informations sur 
les fibres ne pouvant exister independamment de la structure des fibres. 



Modele sujet-observateur 

Nous avons done choisi de les relier sous le modele de l'observateur des design patterns [?]. Nous 

avons done defini ce mecanisme dans le fichier Observer. b.. Le principe est. le suivant : 

«Definit une interdependance un a plusieurs. de facon telle que, quand un objet change 
d'etat, tous ceux qui en dependent en soient notifies et, automatiquement mis a jour.» 

Ainsi on distingue le sujet, qui pour nous est la structure des fibres detenue par la classe FiberStructure ; 

et les observateurs qui sont les informations que Ton peut attacher aux fibres (classe Fiber Inf o<T>). 

On a ajoute un comportement supplementaire : un objet Fiber Inf o<T> ne peut. exister sans un 

FiberStructure, ce pourquoi il n'existe qu'un constructeur qui prend en argument la FiberStructure. 

L'objet. Fiberlnf o<T> va d'ailleurs s'initialiser tout seul en regardant, la structure a laquelle il est 

attache. 

Quand la structure vient, a changer, tous les Fiberlnf o<T> en sont. informes si besoin, et ceux-ci 
vont alors se reinitialiser automatiquement. Quand la structure est detruite, les objets d'information 
sont informes egalement qu'ils ne sont plus valides. 



Structure codee comme un graphe : FiberStructure 

La classe fiberstructure se construit a partir de 1 'image des composantes connexes. L'objet. va alors 
garder en interne une copie de cette image de reference. Elle contient. en interne la structure de 
graphe associee alors aux composantes connexes. Elle contient egalement d'autres informations, 
comme les bary centres des regions, en coordonnees pixel, les surface des regions ... 
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Toutes ces informations sont. accessibles depuis Pinterface de la classe, mais ne sont pas modifiables. 
car elles dependent de la structure meme. 

La classe d'information : Fiberlnfo 

La classe Fiber Inf o<T> est une classe template car les informations peuvent etre de nature tres 
differentes. Elle contient un vecteur STL d'information dont la taille est geree automatiquement 
en fonction de la structure associee. Ce vecteur est un membre public afin d'avoir des acces rapides 
aux informations. La classe sert juste de decorateur 5 pour pouvoir gerer les dependances avec la 
structure des fibres. 

Quand le type du template est numerique, il existe des fonctions pour additionner, soustraire, 
diviser et multiplier deux Fiberlnf o<T>, ou bien un Fiberlnf o<T> avec un element de T : addFI. 
subFI, divFI, mulFI. 



2.4.2 Transformation d'images en Fiberlnfo et reciproquemen t ■;— 

Voir le test unitaire test_f iberProcessing pour les details d 'utilisation. 

Une image peut 6tre transformee en Fiberlnf o<float> avec la fonction d'estimation simple des 
flux : 

int f lowEstimation (Inrimage* image, Fiberlnf o<float> & flow, bool bias = false). 
II faut au prealable initialiser de Fiberlnf o<f loat> flow avec la structure, ainsi Palgorithme 
connait via flow l'image des composantes connexes a utiliser. 

La transformation reciproque est la reconstruction d'image. pour laquelle il n'existe pour le moment 
que la methode de mosaiquage : 

int imageMosaic (Inrimage* image, Fiberlnf o<float> & flow) 

Le principe est similaire a f lowEstimation. La structure est passee via le flow. 



2.5 Soustraction du fond et biais : FiberFlow 

Voir le test unitaire test_f iberFlow pour les details ({'utilisation. 

2.5.1 Fonctionnement general 

La classe FiberFlow implemente le bloc d'estimation des flux de la partie 1.1. sauf la partie 
de transformation de l'image en Fiberlnfo. La classe comporte deux methodes principales : 
setBackground et setlmage, et son comportement est parametre avec la classe de parametre 
FiberFlowParam. La classe permet de soustraire le fond qui est donnee, ou bien il peut l'estimer. 
La classe sert aussi a calculer le biais dans l'image. apres soustraction du fond. 

2.5.2 Le calcul du fond 

Voir le test unitaire test_f iberHisto pour les details d'utilisation. 

Pour calculer le fond, on calcul un quantile de l'histogramme des fibres. Pour cela il existe une 
classe FiberHisto<T> qui prend un Fiberlnf o<T> en constructeur, et qui permet. de faire des 
calculs sur l'histogramme calcule a partir de l'information des fibres. Ces calculs sont plus rapides 
que sur l'image elle-meme. 



'Attention, pas au sens du Design Pattern decorateur [?] ! 
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2.5.3 Le calcul du biais 



Voir le test unitaire test_f iberProcessing pour les details d'utilisation. 

Le calcul du biais n'est en fait pas calcule sur l'image, mais encore une fois sur l'information des 
fibres directement. C'est encore un cas ou on utilise la structure et Pinformation des fibres. Ce 
calcul est fait dans une fonction : 

int biasEstimate ( Fiberlnf o<f loat> & flowln, 
Fiber Info<float> k bias, 
int nbBlocks, 

biasEstimatorType bet=BET_MEAM) 

L'estima.teur du hiais est^our_le_moirient a chQisir_ejatr£Ja_m.Q,y-enne_CBET_HEAN) et. la mediane 

-(BBEJiE&IAN), 

On divise l'espace total des fibres en nbBlocks blocs de taille egale. On parcourt ensuite Pensemble 
des blocs et, on trie les fibres en fonctions du bloc dans lequel elles tombent. Puis, pour chaque 



On calcule le biais ensuite en interpolant pour chaque fibre uniquement (et non pour toute Pimage) 
les valeurs des blocs les plus proches. On utilise pour le moment une interpolation bilineaire, car 
c'est la seule disponible avec Inrimage. On fait une proposition pour avoir une interpolation spline 
cubique. 

2.6 Ameliorations 

- La classe de parametres est a optimiser et a programmer plus proprement. 

- Le biais devrait etre interpole avec une spline cubique 

- Une librairie de graphe serait tout de mane la bienvenue. Cela permettrait d'avoir une structure 
commune entre FiberDetector. AdjGraph, et FiberStructure, ce qui n'est actuellement pas 
exactement le cas. 
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Chapitre 3 

Utilisation d'lmage-Cell 



Auteur CVS (du chapitre) : Author : aymeric 

Version (du chapi£re)T^eFmon'rT.2' ; ;_ ZZZH 

Reference du document : $Id: chapUtil.tex, v 1.2 2002-04-24 07:14:48 aymeric Exp $ 
C'est une partie delicate car elle doit contenir le savoir-faire acquis pendant le prototypage. le 
developpement et les nombreux tests. 

Les petites fibres (sur-segmentees) sont plus genantes que les grosses (sous-segmentees) sur un 
resultat final. 

3.1 Utilisation de la parametrisation 

II faut ouvrir une fenetre d'un editeur de texte (notepad, emacs. wordpad...) avec le fichier ima- 
gecell.txt. Ce fichier est relu a chaque utilisation. La liste des parametre, ainsi que des resumes 
d'indications sont donnees en annexe A. 

3.2 Parametrage de FiberDetector 

3.2.1 Pretraitements 

Dans presque tous les cas. il faut faire le zoom sur l'image pour que l'algorithme fonctionne cor- 
rectement. 

3.2.2 Parametrage de Adj Graph 

Les petites fibres (sur-segmentees) sont plus genantes que les grosses (sous-segmentees) sur un 
resultat final. 

3.3 Parametrage du FiberFlow pour la calibration 

3.4 Parametrage du FiberFlow pour l'image a reconstruire 

3.5 Parametrage d'ImageCell : les reste des parametres 
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Annexe A 



Parametres de Image- Cell 



Les parametres internes d'Image-Cell sont modifiables via un fichier de parametres present dans 
le repertoire source du lancement de Tapplication, au meme endroit que tomoscope.cfg. Le fichier 
des parametres imagecell. cfg a la forme suivante : 
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Ligne du fichie 



IMAGECELLPARAM_CONFIGURATION 



| Defaut | Possibles 



PARAMETERS 



FIBERDETECTORPARAM_CONFIGURATION 



PARAMETERS 



_hDomeWatershed_ 



oversizedFact.or_ 



{0,1} 



{0,1} 



{WO} 



[i;2] 



Parametres de l'algorithme entier 



Parametres de 1' 



Travail sur l'image doublee (plus precis) 



Option non pris en compte pour le moment. 
A regler. 



Hauteur des 
liser les fibre 



selectionnes pour ii 



Taille normalisee (par rapport a la moyenne) 
au dela de laquelle une fibre va etre divisee. 
^encrit^e"esrutilTse~ave(rle cri lere~srrivarrtrsu:r 



a. taiHe (ETlogique) 



_overSizedNeighLimit_ 



{5,6,7} 



Nombre de voisins au dela duquel une fibre va 
etre divisee. Ce critere est utilise avec le critere 



precedent, su r la taille ( K T logique) 



_overSizedResegmentMethod_ 



{0,1} 



•esegmenter dans une fonction distance ; 1 : 
resegmenter dans l'image originale sans pre- 
traitement 



ADJGRAPHPARAM_CONFIGURATION 



Parametres pour la fusion des fibres 



PARAMETERS 



_underSizedFactor_ 



_underSizedMendatoryMergeSlope_ 



underSizedFilterSiz 



underSizedFilterSizeMax_ 



_underSizedFilterCompact_ 



_underSizedFilt.erCompact.Max_ 



[0.1; 1] 



[?;?] 



Taille normalisee (par rapport a la moyenne) 
en dega de laquelle une fibre va etre une candi- 
date pour la fusion. Ce critere est. utilise avec 
le critere suivant sur la taille (ET logique) 



Ce parametre et le suivant determinent 
une droite dans l'espace (nombre de voi- 
sins)x(taille normalisee) delimitant les fibres 
qui doivent obligatoirement etre fusionnees. 
Mieux vaut ne pas y toucher sans faire de cal- 
cul precis ! La droite par defaut passe par les 
points (6; 1/3) et (3,1). Augmenter la pente 
revient k s'occuper plus des grosses fibres avec 
peu de voisins. 



Voir remarque precedente. Ce parametre 
est l'ordonnee a l'origine (axe des ordon- 
nees=tailles normalisees). L'augmenter re- 
vient a obligatoirement fusionner plus de 
fibres. 



{0,1} 



Doit-on filter les fusions possibles en fonction 
de la taille totale de la fusion ? 



[i;3] 



{0,1} 



Si oui, interdire les fusions donnant une taille 
de fibre normalisee superieure a 



[i;io] 



Filtrage des fibres fusionnees en fonction de 1 
compaeite 



Facteur multiplicatif de Pecart-type 
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Ligne du fichier | Defaut 






_underSizedFilterMoreCompact_ 


0 


{0,1} 


Ne garde que les fusions qui rendent les fibres 
fusionnees plus compact.e que la fibre voisine 
seule. 


END ADJGRAPHPARAM CON 


FIGURATION 



END_FIBERDETECTORPARAM_CONFIGURATION 



FIBERFLOWPARAM CONFIGURATION 



PARAMETERS 



_backgroundMode_ 



"3backgroundEstimationCut_^ 



{0,1,2} 



{0,1} 



Parametres d 'estimation du flux pour l'image 
de calibration 



0 : pas de soustraction du fond ; 1 : sous trac- 
tion avec le fond fourni, sinon 2 ; 2 : estimation 
du fond sur Fhistogramme 



■■Quant He-- de 1 ; histogramme-deH4mage-pour- 
Fos t i i 1 1 a t ion dir fond ; 



Correction du biais, apres soustarction du 
fond 

Si correction du biais, nomine de blocks en X 
n Y pour son estimation 



{2,. ,,6 



_biasCorrectionOperator_ 



{0,1} 



Operateur pour la correction du biais : 0 : 



END_FIBERFLOWPARAM_CONFIGURATION 



FIBERFLOWPARAM_CONFIGURATION 


Parametres d'estimation du flux pour l'image 
a reconstruire 


PARAMETERS 




_backgroundMode_ 


1 


{0,1,2} 


0 : pas de soustraction du fond ; 1 : soustrac- 
tion avec le fond fourni, sinon 2 ; 2 : estimation 
du fond sur Fhistogramme 


backgroundEstimationCut _ 


0.01 


[0;0.1] 


Quantile de Fhistogramme de l'image pour 
Pestimation du fond 


_biasCorrection_ 


1 


{0,1} 


Correction du biais, apres soustarction du 
fond 


_biasCorrectionNbBlocks_ 


4 


{2,. .,64} 


Si correction du biais noml u d blocks en X 
et en Y pour son estimation 


biasCorrectionOperator 


1 


{0,1} 


Operateur pour la correction du biais : 0 : 
moyenne, 1 : mediane 


END FIBERFLOWPARAM CONFIGURATION 


_tempoFilterTypeFib_ 


3 


{0,..,4} 


Type du filtre temporel pour Pestimation des 
fibres :0 : moyenne ;1 : moyenne coupee; 2 : 
medianne : 3 : max ; 4 : min 


calibration 


1 


{0,1} 


Effectuer la calibration 


_outputSigmaFilter_ 


3 


[0;10]U{-1} 


Ecart type de la gaussiene pour un filtrage 
passe bas en sortie de image cell. Si la valeur 
est negative (-1), il n'y a pas de filtrage. 



END_IMAGECELLPARAM_CONFIGURATION 
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1 Presentation generale de la section logicielle 



1.1 Specifications techniques 

La section logicielle du tomoscope doit- consister a : 

- permettre le controle d'une partie des composants opto-filectroniques ; 

- recupfirer des donnees (images 2D) ; 

- assurer la gestion, la visualisation, le trait.ement et l'analyse de ces images. 



1.2 Solutions possibles 

- solutionl : Une solution consiste a construire une application logicielle tournant sur un PC communicant 
avec le boitier opto-electronique. 

- solutio n2 : L es con cep te ur s de s y s t emes endoscopiques offrent des systemes embarques per mettant l'ac- 



■ j 3 — Analyse -'. — : — — : — : — 

La conception d'un systeme embarque avec affichage des images sur un ficran TV suffit dans un contexte 
trfis bien determine ou l'opfirat.eur n'a pas le temps (ou la passibilit.fi) d'interagir avec une application trop 
complexe. Le dfiveloppement d'un tel systeme peut fitre menfi par des ingfinieurs possedant des connaissa.nc.es 
a la fois en hardware et en software. II implique la recherche et la mise en oeuvre de carte electroniques, la 
programmation sur OS specialise gfinfiralement bas niveau (assembleur). 

La mise en oeuvre d'un systeme base sur une composante embarqufie (boitier opto-filectonique) connectfi a 
un PC dedie permet a des ingfinieurs software de developper une application modulaire facilement modi- 
fiable. L'utilisation d'un PC standard, utilisant un OS standard permet de se placer dans une configuration 
fivolutive ou les composants dfiveloppfis par des entreprises tiers peuvent fitre utilises. Le travail de proto- 
typage effectue sur un systeme PC/boitier opto-electronique peut fitre transcrit sous la forme d'un systeme 
embarqufi lorsque les besoins sont bien standardisfis. 



1.4 Choix 

Etant dans une phase de prototypage, nous jugeons que la conception d'un systeme embarque, peu modulaire 
ne peut fitre retenue. Cette solution pourra fitre fitudifie dans une phase futur si le design de l'ensemble du 
produit s'avfire stable. II est fort probable que la realisation d'une integration optimale impliquera un accord 
avec un constructeur. 



Nous optons done pour une solution composfie d'un boitier optaelectronique connectfi a un PC sur lequel 
tourne une application logiciel dediee repondant aux services finonces ci-dessus dans la section specifications 
techniques (section 1.1) . 



2 Scenario de mise en oeuvre 

2.1 Scenario de mise en oeuvre et delais associes 

Le scenario de mise en oeuvre du systeme est constitue des etapes suivantes (lesquelles doivent fitre realisfies 
dans 1'ordre) : 

1. connexion du boitier opto-electronique vers le PC : connexion sfirie sur le port COM2, et branchement 
de la carte d'acquisition CNAM embarquee sur le PC avec une nappe au format LVDS ; 

| Confidentiel/Confidential | Manna Kea Technologies, 9, me d'Enghien 75010 PAR.IS, France | page: 3 | 



2. mise sous tension et le demarrage du boiter opto-electronique ; 

3. mise sous tension et le demarrage du PC ; 

4. login de l'utilisateur et demarrage de Implication logicielle du tomoscope. 

Le delai de mise en oeuvre peut done etre considere (pour la partie software du moins) , comme : 

1. temps de boot du PC (si necessaire) : 3 minutes; 

2. temps de lancement de 1'application : 5 secondes; 

3. choix par 1'operateur de la zone disque ou vont etre stockees les images (cf. module Browser mod- 
BROWSER section 4.2) : 1 minute ; 

4. temps du lancement du module d'acquisition modBROWSER (cd. section 4.3) : 5 secondes. 

-Soit un temps cumiile d'environ 4 minutes si le PC doit etre boote ou d'environ 1 min ute si le systemc 
-pC/rroTttef opto-eleetroniqire- e75rdeja"mTS-eTrT)(Ennvrer; 7— ; ; 



2.2 Configuration materielle 

La configuration du PC est la suivante : 

- boitier rack 19" 4U ; 

- ecran plat 17" ; 

- bras de support ecran mobile (2 axes, reglable en hauteur) ; 

- processeur PHI > 700 Mhz ; 

-- disque dur de > 30 Go 7200 TRs ; 

- 256 MO de memoire SDRAM ; 

- lecteur de disquette ; 

- lecteur CD ; 

- graveur de CD ; 

- carte video MATROX G450 16 MO SDRAM dual head (permettant d'utiliser un second ecran) ; 

- carte reseau ; 

- carte d'acquisition CNAM ; 

- clavier industriel (incluant un trackball) ; 

- OS : Windows NT 4 service pack 6 ; 

II s'agit done d'un PC d'un cout approximate de 12 000 FF (sans la carte d'acquisition CNAM) et d'une 
configuration complete de l'ordre de 25 000F. Les elements couteux etant actuellement le boitier rack 19", le 
bras mobile et 1 'ecran plat. L'OS devra evohier vers Windows2000 et eventuellement Linux. 



2.3 Integration hospitaliere 

Pour la pratique hospitaliere on peut envisager de mettre le PC et le boitier electronique sur une colonne 
dediee. L'ecran doit etre suffisamment proche du praticien. Le point, critique est de positionner le connecteur 
de la sonde de telle facon a ne gener ni l'utilisateur du tomoscope, ni une eventuelle personne assise utilisant 
le PC. 
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3 Conception & developpement de l'application logicielle 



3.1 Specifications techniques 

La partie logicielle fournie avec le tomoscope doit. : 

- offrir un service minimal de gcstion dos donnees (images/commentairos) permettant de stocker, visualiser 
les acquisitions realisees avec le tomoscope ; 

- permottre 1'acquisition et. la visualisation en temps reel des images lors d'une session d'acquisition dans 
de bonnes conditions ; 

- pcrmettre a l'utilisateur de visualiser et d'effectuer des post-traitements sur les images sauvegardees lors 
d'une session. 

3.2 — Solutions possibles 

- solution 1 : Les besoins specifies daw \h~mA\mrpf&sti^^mrilfX besdiHs relativement standards: la 
premiere des solutions consisterait a reut.iliser une application existante offrant ces services. 

- solution 2 : Recherche d'un toolkit permettant de developper par nous memo une telle application. 
3.3 Analyse 

- solution 1 : Trouver une telle application, suffisamment, modulable pour y integrer nos propres specificites 
(comme en particulier la communication entre les cartes electroniques et le PC), implique la recherche 
d'un partenariat (consultant/sous-traitant) avec une entreprise fournissant ce genre d'offre. Cela peut etre 
couteux et, aucun choix evident n'apparait au niveau auquel nous nous trouvons actuellement. II est en 
effet, tres difficile aujourd'hui de pouvoir determiner avec certitude les evolutions du soft. 

- solution 2 : Plusieurs toolkits existants sur le marche. Parmi les plus connus : les toolkits Microsoft bases 
sur les MFCs, les produits Borland, des produits gratuits comme t.cl/t.k. Trolltech www.trolltech.com 
propose Qt un framework jouissant d'une bonne reputation (bien documents, prise en main facile, ecrit. 
en C i I-). Le grand avantage de Qt par rapport, aux autres toolkit, est. son aspect multi-plateforme (en 
particulier le code ecrit tourne sur Linux et. WindowsNT). Qt, possede des fonctions pour la localisation 
(choix de la langue). 



3.4 Choix 

Nous avons opte pour un developpement, base sur Qt (www . trolltech . com). Nous essa.yons de garder le plus 
possible 1 'aspect multi-plateforme cependant, les modules d'acquisitions et, de commande de 1'electronique ne 
fonctionnent, pas sous Linux car nous ne disposons pas de drivers pour la carte d'acquisition du CNAM. Nous 
gardons cependant une compatibilite Linux/WindowsNT pour le reste. Le choix de l'aspect, multi-plateforme 
permet. de reduire les couts de developpement, en utilisant les produits gratuits existants sous Linux. Le fait 
de compiler sous deux plateformes different.es permet. egalemenf. de mettre en evidence des bugs. 



4 Description du logiciel 

Les specifications techniques de la partie precedente nous donne directement, le decoupage modulaire de 
l'application : une Interface Generate et, trois modules principaux : 

- modBROWSER : module de gestion des donnees (browser) ; 

- modACQ : module d'acquisition ; 
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- modPROCESS : un module de visualisation/traitement de l'image; 



4.1 Interface Generale 

L'interface principale contient un browser permettant de gerer le stockage et la previsualisation des images 
ainsi que les parametres et commentaires associees aux images (modBROWSER). Elle contient egalement 
deux icones. La premiere permet d'ouvrir le module d'acquisition (modACQ). La seconde permet d'ouvrir 
le module processing (modPROCESS) . La langue utilisee pour l'interface et les commandes est 1 'anglais. 



4.2 Module de gestion des donnees (modBROWSER) 
4.2.1 Specifications techniques 



Le module'de gestion de donnees a pour objectif defaciliter la fache de stockage des images acquises ainsi 
que de leurs parametres d'acquisitions et les commentaires textes associes. II ne constitue cependant pas une 



4.2.2 Solutions possibles 

- solution 1 : Utilisation d'une base de donnees (type Oracle) ; 

- solution 2 : Recherche d'une application existante dans laquelle nous pouvons nous integrer; 

- solution 3 : utilisation d'une arborescence de repertoires classique. 

4.2.3 Analyse 

Un module de gestion de donnees est. generalement compose d'une interface permettant a l'utilisateur d'ac- 
ceder a un base de donnees. Les images obtenues avec le tomoscope constituent une partie d'une information 
qui devra etre associee a un patient. Ces images vont etre soumises a des contraintes medico-legales. Ces 
images devront etre egalement accessibles eventuellement par plusieurs personnes (eventuellement hors-site 
pour une analyse conjointe par un gastroenterologue et. un anapath). L'ensemble de ces contraintes et be- 
soins montre l'importance du module de gestion des images. Plusieurs entreprises developpent des solutions 
pertinentes pour ce genre de problemes (par exemple Etiam pour les aspects Dicom et, la telemedecine, le 
Pack-PMSI de i-care pour la gestion administrative du dossier patient). 

Les besoins a l'etape actuelle sont relativement simples : il s'agit d'offrir a l'utilisateur un minimum d'orga- 
nisation le guidant dans l'archivage des donnees et la. post-analyse de ceux-ci. 

4.2.4 Choix 

Considerant. que les besoins a l'etape actuelle sont relativement. simples : il s'agit d'offrir a l'utilisateur un 
minimum d'organisation le guidant. dans l'archivage des donnees et. la post-analyse de ceux-ci, nous optons 
pour la constitution d'un browser sur une hierarchie de repertoires predefinie. 

Le browser permet de visualiser une hierarchie composee de 3 niveaux, nominees par exemple 

- Patient : identifiant du patient 

- Examen : identifiant de l'examen 

- Study : identifiant de 1 'etude 

Les images et commentaires associees obtenues lors de l'acquisition k sont done stockees dans le repertoire 
Patientp/Examene/StudyA,. 
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Les noms des differents repertoires peuvent etre changes, la profondeur de la hierarchie peut egalement, etre 
modifiee en fonction des besoins. 



Scenario d 'utilisation 

Avant, de commencer une acquisition, l'utilisateur renseigne les champs mentionnes plus haut. Les images 
obtenues lors de Facquisition seront alors sauvegardees dans le repertoire correspondant a ces champs. Les 
images d'une meme session sont done tout.es sauvegardees dans un meme repertoire. Les noms des images 
seront generes automatiquement ou specifiees par Foperateur. 



Configuration materielle 

Les-oWiees-Sonl-stodceea-&urJ.e„disque dur du_IjC_d-acquiiiition. Le simple partage du repertoire de sto- 

-eka g c p e rm e ttra ^r^.c^s ^omm 4ir^in-s-iin--r-fcsftan- locaL- 

Nous effectuons une veille pour la recherche de solutions plus completes pour ce module (juge pen prioritaire a 
Fheure a~ctrrelie)r^ 

4.3 Module d'acquisition (modACQ) 

4.3.1 Besoins 

Les different.es concertations avec les experts medicaux, ont conduit aux specifications suivantes : 

- la procedure d'acquisition consiste a poser la tete endoscopique sur la zone d'interet, rester immobile a 
peu pres une seconde, declencher la sauvegarde de 1 'image ; 

- le resultat d'une session d'acquisition est un ensemble d'images 2D, chacune des images correspondant a 
un echantillon selectionne par le praticien operateur ; 

- la selection des images se fait par rapport, a la qualite diagnostique de celles-ci ; 

- il n'est pas necessaire d'enregistrer en continu les acquisitions fait.es par le tomoscope : 

Scenario d'utilisation 

1. les images acquises par le tomoscope defilent. sur Fecran (10 images/secondes) ; 

2. le praticien pose la tete endoscopique sur une zone d'interet, ; 

3. si le praticien juge que les images qu'il voit sur Fecran sont interessant.es, il commande un arret de 
Fimage; 

4. le systeme dispose alors d'un petit intervalle de temps (maximum 5 secondes) pour appliquer des 
traitements rapides visant a ameliorer la qualite de Fimage. Les traitements a appliquer sont, decides 
par les experts MKT ; 

5. Fimage de qualite superieure apparait a Fecran ; 

6. le praticien peut, sauvegarder Fimage ou decider de continuer; 

7. si le praticien decide de sauvegarder Fimage, cette image est stockee automatiquement dans le repertoire 
approprie. Le praticien (ou son assistant) a la possibilite d'associer quelques commentaires text.e a cette 
image ; 

4.3.2 Specifications 

A partir des besoins decrits ci-dessus, on peut, deduire un ensemble de fonctionnalites pour le module d'ac- 
quisition : 

1. reception des images depuis la carte d'acquisition du CNAM ; 

2. commande de Felectronique CNAM (carte PM, carte de modulation, carte de synchronisation) : 
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3. affichage en temps reel des images ; 

4. possibility de bufferiser quelques images; 

5. commande de controle de 1' affichage/ sauvegarde ; 

6. traitement, temps reel ou quasi temps reel ; 

Reception des images depuis la carte d'acquisition du CNAM : 

Lorsque Ton est dans le module modACQ le tomoscope effectue en permanence des acquisitions d 'images. 
Les images sont constitutes de mesures de taille 640x640 (TBC), codees sur 8 bits (TBC). Ces mesures ne 
constituent pas des images isotropes (cf. section 5.2.3). On garantie une acquisition au rythme de 10 images 
par seconde. Physiquement le transport de donnees se fait par une nappe au format LVDS connectant le 
boltier opto-electronique a la carte d'acquisition CNAM sur bus PCI du PC. 



Commande de l'electronique-GNAM-: 

Une interface dans le module d'acquisition permet. de changer les parametres de l'electronique (cf. document 

-du^N-AM^-Xe^cont^ 

controle sera reserve au mode expert. Physiquemont. on utilise le port, serie (COM2) et. une nappe au format 
LVDS pour commander respectivement les cartes electroniques et recevoir les images. 

Affichage en temps reel des images : 

On affiche en temps reel les images brutes (c'est a dire sans traitement) 640x640 (TBC), 10 images/seconde 
(TBC). devolution de cet affichage consistera a apporter les traitements temps reel permettant d'ameliorer 
la lisibilite de l'image par le praticien. 

Possibilite de bufferiser quelques images : 

L'expert MKT peut decider d'un nombre d'image que Ton garde en memoire pour des besoins de traitement 
eventuel. La modification de ce parametre est possible lors de la phase de test. A terme, le contrdle de ce 
parametre sera reserve au mode expert, et preregle automatiquement. 

Commande de contrSle de l'affkhage/sauvegarde : 

Cette section permet d'offrir les services necessaires aux parties 6 et 7 du scenario d'utilisation (cf. section 
4,3.1). L'utilisateur dispose d'une commande deportee (type pedale, joystick, -TBD-) lui permettant de : 

- commander l'arret sur l'image courante ("toggle pause") ; 

- commander la sauvegarde de l'image courante. L'image est alors sauvegardee dans le repertoire specifies 
par le module modBROWSER. 

- lors de la sauvegarde de l'image l'utilisateur (ou son assistant) a la possibilite de sauvegarder un commen- 
taire associe. 

Traitement temps reel/quasi temps reel (TTRQTR) 

• Specification technique : Afficher en temps reel une image de qualite superieure a celles des images 
brutes provenant de la carte d'acquisition. 

• Solutions possibles : II peut s'agir d'effectuer des traitements purement software comme des manipu- 
lations de LUT par exemple permettant d'ameliorer la lisibilite des images. L'autre aspect concerne les 
asservissement de l'electronique a. partir de mesures effectuees sur les images acquises visant a optimiser 
les parametres de controle de l'electronique en fonction de la profondeur par exemple. Ces deux solutions 
sont bien entendue complementaire. 

• Analyse : Les images brutes provenant. du tomoscope contiennent un certain nombre de defauts qui 
peuvent gener le praticien dans la lecture et l'interpretation de l'image (manque de contraste par exemple, 
se referer aussi a la section 4.4). De plus l'electronique de balayage peut se deregler et conduire a des 
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acquisitions de mauvaises qualities. 



• Choix : Le module d 'acquisition doit posseder Ies traitements suivants : 

- TTRQTR 1 : correction automatique du gain : permettant d'ajuster automatiquement le controle du 
gain en function de la profondeur. Cette operation s'effectue par analyse de l'histogramme et consiste 
en une boucle de controle entre les images acquises et les cartes electroniques. 

- TTRQTR 2 : calibration en amplitude. C'est une operation qui implique une phase de calibration 
effectuee en debut d'acquisition et permet de s'affranchir du pattern des fibres (cf. section 4.4). 

Traitement rapide effectue lors d'un arret sur image (TREAI) 

• Specification technique : Afficher rapidement. (delai inferieure a 5 secondes) une image de qualite 
superieure a celle de l'image brute provenant de la carte d'acquisition lors d'un arret sur image. 

. Solutions possibles : Les imagesTTute^o^nanTWl^^^ un certain nomhre de 
- defauts qui peuvent gener le praHe»en-dans4a-lee tur c o t 1'intRrprctation de l'image (cf, section 4.4) 

-•-Analyse : cf. section 4.4. ; __i ; — - _ 

• Choix : Le module d'acquisition doit posseder les traitements suivants : 

- TREAI 1 : attenuation du bruit poissonien par moyennage temporal : il s'agit d'une operation qui a 
pour objectif d'attenuer un bruit intrinseque a l'acquisition afin de rendre plus facilement lisible l'image. 
Cette operation implique que plusieurs images soient bufferisees. 



Notes Ces traitements rapides permettent d'ameliorer la qualite de l'image. Us sont effectues uniquement 
sur demande de 1'expert MKT a partir d'une ou de plusieurs images bufferisees. Lorsque les traitements ne 
sont pas temps reel, ils interrompent l'acquisition. 

4.4 Module Processing (modPROCESS) 

Ce module permet de visualiser une image stockee dans la hierarchie et de lui appliquer certains traitements. 
Ces traitements sont effectues en post-processing, sans contrainte de temps reel, mais avec une contrainte de 
rapidite suffisante pour ne pas ralentir l'examen. Un certain nombres des modules de traitements presents 
dans cette section peuvent etre utilises pour les traitements lors d'un arret sur image (les TREAIs). 

5 Traitements des images 

Cette partie decrit les differents traitements realisees sur les acquisitions. Les traitements sont repartis en 
deux groupes : les Traitements Temps Reels / Quasi Temps Reels (TTRQTR) , et les Traitements Effectues 
lors d'un Arret sur Image (TREAI). Ce decoupage ne concerne pas directement l'organisation des modules. 
Les TTRQTR ne sont realises que dans le module d'acquisition. Les TREAI peuvent, suivant les besoins, 
etre realises dans le module d'acquisition, ou dans la module processing. 

5.1 Traitements Temps Reels / Quasi Temps Reels (TTRQTR) 

5.1.1 Correction automatique du gain 

Specification technique 

Suivant la profondeur et. la nature du tissus, la valeur de la puissance regue sur le detecteur peut varier. La 
chaine de detection prevoit plusieurs moyen de controler ce niveau lors de sa numerisation. Afin d'etre dans 
des conditions optimales du point de vue image, ce gain doit etre controle automatiquement en fonction de 
l'image acquise. 
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Solutions possibles 

Le gain peut etre controle via la puissance laser, la carte PM et. le gain propre du detecteur (si il est 
reglable). II y a deux facons de controler le gain : une boucle de retroaction "semi-aveugle", ou bien en 
anticipant le resultat sur l'image. La premiere solution propose juste de diminuer le gain s'il est trop fort, 
et de Paugmenter si il est trop faible. La seconde solution trouve la bonne valeur du gain par un calcul qui 
incorpore un modele numerique de la chaine d'acquisition et de controle. 

Analyse 

La boucle de retroaction est la solution classique quand une partie du systeme n'est. pas bien connue. 
La difficultc du controle depend alors de l'ampleur de la non-connaissance de cette partie. II devient plus 
difficile de prodire le temps de reponse et. la stability du controle. La seconde solution requiert une bonne 
connaissance de la chaine d'acquisition et. des moyens de controle de cette chaine. On peut alors anticiper par 
1 'observation de l'histogrammo de l'image la valeur correcte du gain a appliquer. II suffit. ensuitc do traduire 
cette valeur du gain en termes de^ontrfile dTl^chaTne d'acqui^tT6"n:i^eTt^"derniere met.hode a~T avant,age 
de-converger en une seule iteration de-contrele-Elle sera-done- plus rapide dans tous -les cas que laboucle de 
retroaction. 



Choix 

Le choix se porte sur la seconde solution : par anticipation de la boucle de controle. II faut bien noter que 
ce choix requiert une bonne connaissance de la chaine de controle. 

5.1.2 Calibration 

Specification technique 

Le phenomene d'injection du spot, laser dans la fibre combine a l'echantillonnage des spots fait apparaitre le 
pattern des fibres sur les acquisitions. Une reflexion parasite sur la partie distale du guide va etre enregistre 
par les detecteurs. L'amplitude de chaque point de mesure de l'image est soumis k ces phenomene opto- 
mecaniques qui va creer une texture genante pour la lisibilite de l'image. Une procedure de calibration est 
ici envisagee pour eliminer ces phenomenes. 

Solutions possibles 

Le probleme est complexe car on ne peut pas tout eliminer en meme temps. Ce que Ton peut faire : 

1. soustraire l'image de la reflexion parasite ( obtention de celle-ci quand le guide est hors focus) 

2. soustraire l'image d'un milieu diffusant homogene 

3. soustraire l'image d'un miroir 

Analyse 

La solution 1 est simple a mettre en oeuvre car elle necessite tres peu de manipulations. Elle ne permet 
que d'enlever un fond sur l'image. Comrae cette soustraction se situe apres la quantification (8 bits) elle 
ne permettra en aucun cas de voir des choses que Ton ne pouvait pas voir avant. Cela peut juste aider la 
visualisation La solution 2 necessite un contact optimal avec un milieu liquide. Soit par immersion, soit avec 
une interface. Cette solution permet de retirer le plus de defauts (echantillonnage, guide, systeme optique...). 
Par contre sa mise en ceuvre est tres difficile. La solution 3 suppose un alignement. parfait. entre le miroir 
et 1'axe optique, ainsi qu'une courbure de champ plate au niveau de la tete endoscopique. On ne connait. 
pas encore la courbure, et les problemes d'alignement devraient etre realises au moyen de pieces mecaniques 
de haute precision. Les solutions 2 et 3 fond perdre la sterilite de r'appareil, sauf si la mire ou le milieu 
sont eux-memes steriles, done jetable probablement. La solution 1 est la plus souple. Une combinaison des 
solutions 1 avec 2 ou 3 est, eventuellement possibles. Les solutions 2 et 3 sont plus difficilement compatibles. 
Quelle que soit, la solution, il faut bien voir que le vrai facteur limitant pour nous est, la quantification 8bits. 
II ne faut done pas attendre de miracle de la soustraction... 
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Choix 

La solution 1 est la moins contraignante. Et. en l'absence de miracle on choisit la moins contraignante. 
(TBC). 

5.1.3 Calcul automatique du horizontal time shift 

Specification technique 

Le horizontal time shift (HTS) est un parametre qui regie le debut de 1'echantillonnage des lignes ([?]). 
Ainsi une mauvaise valeur du HTS va creer un entrelacement des lignes. La valeur de ce parametre varie en 
fonction des miroirs et des conditions externes de son mouvement (temperature, pression...). Ce parametre 
doit etre regie automatiquement. pour qu'il n'y ait. pas d'entrelacement. 

Solutions possibles 

— I^enfi^lacemenrpimretr^^ 



Analyse 

La presence des fibre s va c reer une composante hau te frequence (spati alement....) q ui v a permettre d'obt.enir . 

un pic de correlation 6!cve. Le pic de correlation est unique car le guide est semi-regulier et on ne pent pas 
se retrouver dans une configuration periodique sur une ligne horizontale. 

Choix 

Le choix est rapide puisqu'il n'y a qu'une seule solution. 

5.2 Traitements Effectues lors d'un Arret sur Image (TREAI) 

5.2.1 Augmenter le SNR 
Specification technique 

Les acquisitions sont entaches de plusieurs bruits : le bruit poissonien intrinseque au comptage de photons, 
les bruits de la chaine electronique, de la quantification... Le bruit majoritaire est le bruit poissoniens aux 
niveaux de puissance ou on travaille. Ce bruit gene la lisibilite de 1'image et. doit, etre reduit. 

Solutions possibles 

Les solutions envisagees ne prennent pas en compte les modification possible des elements optiques et/ou 
electronique d 'acquisition. La reduction du bruit peut se faire de deux facons : soit par un filtrage temporel, 
soit par un filtrage spatial. 

Analyse 

Le filtrage temporel de type passe-bas permet de rejeter les frequences hautes qui sont majoritairement 
du bruit dans le cas d'une image fixe. En l'absence d'information sur les images (bouges,...), on propose un 
moyennage uniforme. Les traitements spatiaux sont plus complexes et dependent du contenu de 1'image, et 
de reformation que Ton peut obtenir a priori sur ce contenu. Dans notre cas, les hautes frequences spatiales 
sont occupees par le bruit, mais aussi par le pattern des fibres. Un filtrage passe-bas spatial rendrait 1'image 
floue, ce qui n'est pas acceptable. 

Choix 

Le choix se porte done naturellement sur un moyennage temporel uniforme. 

5.2.2 Suppression du pattern des fibres 
Specification technique 

Fournir une image non pixelisee : s'affranchir du phenomene d'injection dans les fibres. La lecture de 
1'image ne doit pas etre genee par rapparition du guide d'image sur celle-ci. 
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Solutions possibles 

Trois solutions sont envisagees. La premiere est. un filtrage spatial de type passe-bas qui utilise la forme 
frequentielle des fibres. La seconde est. une analyse de la forme des fibres qui permet de retablir un niveau 
d'injection constant sur toute la fibre. La troisieme solution est. une calibration sur une mire uniforme connue 
qui permet de connaitre, pour chaque point de mesure, son attenuation globale dans la chaine d'acquisition. 

Analyse 

La premiere solution a l'avantage de pouvoir etre implant.ee rapidement a l'aide de filtres recursifs sepa- 
rables pour un coiit de 15 operations par pixels, c'est, a dire presque en temps reel. De plus les parametres 
du filtres sont determines automatiquement en fonction des specifications techniques du guide d'image. Le 
traitement fait bien disparaitre le pattern, mais fait apparaitre un artefact reparti sous la forme d'un bruit 
moyenne frequence isotrope visible dans les zones homogenes. Ce filtre est deja contenu dans la bibliotheque 

inrimag e. 

La seconde solution propose u n resultat diffe rent. Le principe est d'isoler chaque fibre a l'ai de de techniques 

-de morphologie Tnathematique. Gette-technique utilise -la forme- de la fonction d'injection 2D, ainsLque-la — 
reflexion parasite en bout, de guide. On peut ainsi isoler sans calibration et pour chaque image individuelle, 

^haqu^fibr^u^ide-d^mag^ — 
avec une analyse statistique locale. Le temps de traitement envisage pour chaque image (512x512) est d'au 
plus 0.5 secondes. L'image resultat est une sorte de mosai'que. Le cout. de developperncnt est. d 'environ une 
semaine. 

La troisieme solution permet de compenser globalement l'attenuation de chaque point de mesure. Une hypo- 
these forte doit etre verifiee : le balayage et Pechantillonnage doivent etre tres stable (inferieure a 500nm). 
Dans ce cas, on peut corriger chaque point, de mesure a partir des donnfees calibree. L'image resultat est 
proche d'une image sans guide d'image, car on simule une injection presque parfaite en chaque point de 
mesure. II peut rester un bruit de quantification pour des points ou 1 'injection etait mauvaise. Cette methode 
necessite tres peu d'operation par pixel, et. un seul scan de l'image (done temps reel). Cette solution a deja 
ete developpee dans le simulateur 

Les deux premieres solutions ont 1'avantage de n'etre que logicielles. La troisieme necessite une manipulation 
de calibration avant l'examen. Sous reserve d'une tres bonne stability du balayage/echantillonnage, et d'une 
bonne injection sur l'ensemble du guide, la troisieme methode devrait. fournir de tres bons result.ats. Cette 
methode doit cependant etre sensible aux hypotheses de stabilite. Des criteres imperatifs de robustesse 
demanderait alors plut6t une des deux premieres methodes (TBC). 

Choix 

TBD. De nombreux tests et. l'avis de plusieurs medecins seront. necessaires. Notamment la stabilite du 
balayage doit etre verifiee sur une echelle de temps importante (de l'ordre de 2h00) et en situation reelle 
pour la calibration systematique. La troisieme methode semble la plus robuste et est le choix par defaut. 

5.2.3 Correction des distortions geometriques 
Specification technique 

La forme specifique du balayage combine a l'echantillonnage ainsi qu'a la chaine de traitement optique 
et electronique induit des distortions non lineaires de l'image. Ces distortions doivent etre corrigees afin de 
presenter une image non deformee (espace image isotrope). 

Solutions possibles 

Deux solutions sont envisagees. La premiere est. une modelisat.ion analytique de la chaine optique-electronique. 
La seconde solution est une estimation globale des distortions a partir d'une mire (observation d'une grille de 
taille connue). Ces modelisations permettent de deduire la transformation et. de realiser un re-echantillonnage 
bilineaire de lMmage dans une grille pixel carree connue. 

Analyse 

La premiere solution (analytique) est. simple a mettre en ceuvre car la forme theorique du balayage et celle 
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de l'echantillonnage tempore! sont, bien connues. Cette methode ne demande pas de calibration. Le resultat 
n'est cependant pas garanti car on ne peut pas prendre en compte toutes les distortions de la chaine. 
La seconde solution demande une calibration initiale de l'appareil, et une verification periodique (TBC). 
On image une grille connue, et on effectue le recalage de toutes les intersections de la grille par rapport a 
un modele analytique de cette grille. Cela permet de quantifier le deplacement de chacun de ces points et. 
d'interpoler une deformation entre ces points. La distortion etant a priori spatialement lentement variable, 
son evaluation en differents points d'une grille de calibration doit, donner un tres bon resultat. La calibration 
devrait de plus etre entierement automatique (TBC). Cette distortion est estimee une fois pour toute pour 
un jeu de parametres donne, et non a chaque examen. 

Choix 

TBD. La solution de calibration sur une grille devrait etre la meilleure et devra au moins etre realisee en 
laboratoire pour co nfiimcr le modele thjorique. __ 



5:274 — v ._ gl1 ^_ H __ m , fe . diagnost i ( . — — 

Specification technique .. . _____ 

L'image doit etre presentee de facon a ne pas presenter d'artefacts de traitement "(ou a defaut ils doivent 
etre reperes), et de facon a faciliter l'interpretation. 

Solutions possibles 

Les solutions concernent la visualisation de l'image deja acquise. Certaines solutions classiques sont em- 
pruntees au monde de la radiologie conventionnelle numerique. Les traitements envisagees sont : fausses 
couleurs (proches des colorants classiques utilises par les anathomo-pathologist.es) , fen 6t rage lineaire et lo- 
garithmique des niveaux de gris, zoom numerique. 



Les fausses couleurs permettent. selon l'observateur et. les conditions de luminosite de la piece, de se placer 
dans des conditions connues d 'observation. Le fenetrage lineaire permet de modifier le contraste de l'image 
en fonction de la plage de niveau de gris. Le fenetrage logarithmique permet de privilegier le contraste des 
zones sombres par rapport a celui des zones claires. Sur une image de type gradient (comme les notres), cela 
permet de faire attention aux zones de transition faibles, par rapport, aux zones de transition fortes (qui sont 
des informations evidentes a priori pour le praticien (TBC). Enfin, le zoom numerique permet de grossir 
certaines zones d'interet. 

Ces traitements de visualisation sont faciles a mettre en ceuvre et sont des outils largement repandus. De 
plus ils impliquent des traitements simples qui ne creent pas d'artefact. 

Choix 

Tous ces outils sont implantes. mais leur utilisation est opt.ionnelle. Un mode de visualisation par defaut 
est propose. 
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Reference du document : $Id: imageCellSA_cdc.tex, v 1.1 2002/11/04 13:24:19 aymeric Exp $ 

1 Introduction 

Ce document est un mini cahier des charges pour la construction d'une interface Qt pour pouvoir 
utiliser facilement imagecell. Le but est que plusieurs personnes de MKT puissent trailer des images 
de fagon autonome. Le fonctionnement standard d'imageCell ne demande que peu de changements, 
et qui sont. faciles a maitriser. 

2 Interface texte existante 

II existe deux versions d'imageCell : une version pour reconstruire des images, et une autre pour 
les films. 

2.1 ImageCell, images fixes 

L'application est dans SoftProj/testApp/imageCell/imageCell. cpp. 

Usage : imageCell -i <image> [-b <image>] -f <image> I -g <image> [-c <image>] 
[-o <image>] [-u <image>] [-p <file>] [-d <file>] 

-i <image> image 

-b <image> background image 

-f <image> fiber image 

-g <image> fiber image component 

-c <image> background fiber image 

-o <image> output image (default imageCell_out . inr) 

-u <image> output fiber image (default: no output for fiber image) 

-m <image> output fiber component image (default: no output for fiber image) 

-p <file> parameter file 

-d <file> dump point coordinate and flow into file 

Entrees: (-i, -b, -f, -g, -c, -p) on rentre l'image (accompagnee d'un fond ou non), 1'image 
des fibres (accompagnee d'un fond ou non egalement), et un fichier de parametres. L'image des 
fibres peut prendre deux formes : soit la forme brute (-f , -c), soit l'image des composantes 
connexes des fibres (-g). 
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Sorties : (-o, -u, -m, -d) on a au moins la sortie de l'image courante (-o). On peut egalement 
sortir l'image des taux d'injection (-u), l'image des composantes connexes (-m), l'image des 
points (-d) : on a pour chaque fibre, les coordonnees x,y,z (coord, axiale ! la tranche tem- 
porelle), et le flux mesure, soit un quadruplet, par fibre, a raison d'une fibre par ligne en 
ASCII. 

L'image de sortie est 2D et 8 bits. L'avantage des options -m et -g est. la rapidite de traitement : 
cela evite d'avoir a recalculer la detection des fibres plusieurs fois. 



2.2 ImageCellMovie, films 

L'application est dans SoftProj/testApp/imageCellMovie/imageCellMovie.cpp. 
-Us ag e : i -mageGell-MQV-A-e-— i^-<4mag-e>— [-b--<-image>-]— --f--<lmag.e>--[- ^c < i m aged— 



-i <image> image 

-b <image> background image 

T^f <imageS 

-g <image> fiber image component 

-c <image> background fiber image 

-o <image> output image (default imageCellMovie_out . inr) 

-u <image> output fiber image (default: no output for fiber image) 

-m <image> output fiber component image (default: no output for fiber image) 

-p <file> parameter file 

-d <file> dump point coordinate and flow into file 
C'est la meme interface que precedemment, sauf que la sortie est 3D. 



3 Interface future 

3.1 Generality 

L'interface doit permettre d'utiliser les fonctionnalites de l'interface en ligne de commande, avec 
en plus la possibility de visualiser les images qui rentrent dans le calcul. Les parametres doivent 
egalement etre modifiables et recuperables. Pour connaitre leurs roles, on s'appuie sur derniere 
partie. Le meme document explique egalement ['implementation d'imageCell. 
La seule difference entre imageCell et imageCellMovie est, au niveau de la sortie qui est 2D ou 3D. 
On peut done integrer les deux dans la meme interface, meme si la realisation utilisera deux ob jets 
distincts imageCell et imageCellMovie. 

L'utilisation des objets est mise en ceuvre dans les deux programmes en ligne de commande, et il 
suffit de s'appuyer dessus... 

3.2 Parametres 

Parmi les parametres, certains sont plus importants que d'autres : 

_doDiffusion_ : prefiltrage pour reduire le bruit (pas encore implante, mais a prevoir), 
_biasCorrection_ : correction du biais (en double : pour l'image des fibres, et pour l'image a 
reconstruire) 

biasCorrectionNbBlocks_ : finesse de la correction du biais (si besoin) : entre 2 et 20. 
tempoFilterTypeFib : filtrage tempore] en entrees (pour images fixes) moyenne, moyenne 
coupee, mediane, max, min 
_ca!ibration_ : calibration des taux d'injection 
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_outputSigmaFilter_ : filtrage passe bas en sortie. 

Ces parametres doivent. apparaitre clairement. Les autres doivent, juste etre accessibles, au cas ou 
et par soucis d'exhaustivite. 

References 

[1] Aymeric Perchant. Manuel de reference pour imagecell. Technical report, Mauna Kea Techno- 
logies, Base CVS, SoftTex/UtilisationSoutien/ManuelReference/R.eferencelmageCell.tex, 2002. 
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introduction 

Les images acquises par le tomoscope sont bruitees. La principale composante de ce bruit semble etre le 
bruit poissonnien issu du photodetecteur. D'autres source de bruit sont envisages comme le speckle, ou bien 
le bruit d'autres appareils (electronique, opt.iques, quantification, controle du gain...)- Ce document a pour 
but de trouver la meilleur operation de type moyennage permettant a partir de plusieurs aquisistions fixes 
"(pas de~bouger;rde diminuer la variance du^g^alrtfoTTs^^oTTs~am 

Ides^op.er.ation„essentiellerrient temporelles. . 

Le principe de reconstruction de 1'image a partir du jeu de mesures fait, intervenir plusieurs spots dans un 
pixel. Le prob lemo ici est de savoir co mment calibrer et combine;- ces different s spots pour potivoi r obtenir 
le pixel le moins bruite. Le deuxieme probleme est d'optimiser le moyennage tempore! de plusieurs jeux de - ■ 
mesures. Ces deux problemes sont evidemment. lies. 
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1 Cas Simple 



On veut. calculer une calibration en intensite et, un mapping geometrique pour corriger respectivement, la 
modulation du signal par le guide d'image et la distorsion geometrique induite par la vitesse sinusoidale du 
miroir ligne. 

Le cas le plus simple consiste a considerer que Ton a pas de bruit sur la mesure au niveau du PM, que le 
balayage est reproductible et que l'on peut calculer la position du spot dans le repere de balayage. C'est cette 
methode qui est implementee dans la version actuelle du Simulateur SimuApp (cf. CalibrationMapping et 
Reconstruct ionMapping). 

On considere : 

- cal(k) : le k' Eme element du flot de donnees obtenu en observant un miroir. II s'agit d'un ensemble de 
512X512 elements par exemple. 

- obs{k) : le k a ' eTOe -etemeiil. du flol. de donnees obtenu en obsen^t^^imnfcition-d%rt6r6fc - 
_ • caP(k) (on pbs J (k)):. le k tem ^lemenrc!U^c^de _ t , arquisiti"oii j r Iritulrce" j repT-feErrteTic^^ ; 

porel. 

Initialiser les elements de 3 vecteurs : 

pour tout k 

X[k] = Y[k] = coeff [k] = INVALIDE 

pour tout k 
si cal[k] != 0 

calculer l'indice du pixel (X[k], Y[k]) de 1' image reconstruite dans lequel tombe ce spot, 
si cet indice est valide : 

coeff [k] = (l/cal[k]) si on veut faire une calibration en amplitude 
coeff [k] = 1.0 sinon 
pond(X[k] , Y[k]) +=lj 

pondQ est une image de meme taille que l'image reconstruite qui indique le nombre de spot valide qui tombe dans 
chacun des pixels. 

Pour tout element du flot calculer : 

coeff [k] *= 1.0/pond(X[k] , Y[k]) 
La calibration calcule done 3 vecteurs : coeff[k] , X[k| , Y[k]. 

Reconstruction : 

Pour tout element du flot obs(k) : 
si X[k], Y[k] sont valides: 
si coeff [k] est valide 

image (X[k], Y[k] ) += obs [k] *coef f [k] ; 
Dans le cas ou Ton dispose de N images d'acquisitions et de K images pour calculer la calibration. On 
remplace les flots obs(k) et cal(k) par leur moyenne temporelle : 
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2 Modele 



2.1 Distribution des amplitudes dues au guide d'image 

Dans cette etude, le terme amplitude designe la valeiir mesuree du signal en sortie du tomoscopo. Ce qui 
nous interosse est dans le modele de bruit est le rapport entre l'amplitude regue si on avait ni guide d'image, 
ni trou de filtrage, et l'amplitude theoriquoment. recue en presence do ces deux composant.es du systeme. 
Nous avons d'abord etudie la distribution des amplitudes maximales do l'echantillonnage du guide d'images. 
La trajectoire du laser passe par l'entree du guide d'image. ce qui module l'amplitude do la puissance du 
spot en fonction de 1 'intersection des disques (spot et. fibre optique). Puis, 1'amplitude est oncore modifiee 
au retour par le trou de filtrage on fonction oncore une fois de Intersection de disquos. 
Le balayage laser est suppose a puissance (pulsee) identiquos partout. Lo systeme de balayage va creer une 
certaine distribution d'amplitudos dues aux deux attenuations, dont la distribution do probability d'' amplitude 
-^st^orfflee^su-r-kJi^a-r-p^ 



Probabilite d'amplitude 
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Fig. 1 - Distribution de probabilite des amplitudes due au phenomene d 'acquisition via l'echantillonnage du 
guide d'image, puis le trou de filtrage, avec les parametres par defaut du tomoscope. Cette distribution est 
estimee a partir de l'histogramme des mesures simulees. 

Cette distribution a ete calculee a partir de donnee non bruitees issues du simulateur, pour les parametres 
par defaut de celui-ci, ainsi que pour une normalisation de l'amplitude de sortie a 200. La figure 2 represents 
la fonction de repartition de la distribution de probabilite p„.. C'ost-a-dire pour une amplitude donnee X n , 
on represente P(x < Xq), 

Dans les experiences qui snivent, nous ut.iliserons cette probabilite p n afin de simuler plusieurs amplitudes 
theoriques, c'est-a-dire on bruitees, de plusieurs spots. Le principe de reconstruction de 1 'image a partir du 
jeu de mesures fait intervenir plusieurs spots dans un pixel. Les amplitudes theoriques de cos spots seront. 
calculees avec la distribution p n . 



2.2 Rappels sur le bruit poissonnien du systeme 

Nous tenons compte id des bilans photometriques expliques dans le document du simulateur (Georges Le 
Goualher, Gregoire Malandain, Aymeric Perchant), ainsi que dans le bilan de fin d'annee de Magalie Surivet. 
Les notions sont les memes qu pour le document du simulateur. Ce bilan peut se resumer ainsi. 




Bilan photometrique 


Parametres 


Id 


Defaut | Domaine 


Puissance dans le guide d'image 




770//.W 


[0.01, 5]mW 


Rapport cyclique 


T,„ 


4 




Coefficient de reflexion parasite 


r p 


0.55 10" 4 




Coefficient de retrodiffusion maximale 


R Tn „r 


0.001 




Coefficient de reflexion au point considere 


R(x) 




[0,017; 0,04] 


Coefficient d'attenuation de la retrodiffusion 




0.103447 




Energie d'un photon 


E1A(A) 


2.89 10-' 9 





Le nombre de photons qui vont etre compte par le photomultiplicateur est au total de : 
N(R(x)) = N rp + N utUe = N guide (r p + a r * R max * JS(as)) , 



N Tp = N guide *r p , 
NutMe = N guide * a r * R max * R(x), 

PguideT e mission T puUe . 



N < 



T e . m .is 



Le signal N(R(x)) est. alors entache d'un bruit, poissonnien. Le signal bruite est. note N(R(x)). Ce signal 
comporte un offset par rapport a la partie utile que nous retirons, puis nous appliquons une correction du 
gain aver. : 

_ N(R(x)) - N gvMe *r p + k*VN 

a r * Rmax 



R.(X) = - 



(1) 



afin de retrouver un estime note R(x) de R(x). Le facteur entier k regie un offset supplemental exprime 
en terme d'ecart-type du bruit. Ceci permet de ne pas saturer le signal si le bruit est, important. 
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2.3 Mesures bruitees 

Un pixel est. rer.onsit.u6 a partir en moyenne de 8 mesures de spots pour les parametres par defaut. Nous 
noterons N spot / pixe i ce nombre. Chacun de ces spots a une amplitude modulee par le systeme guide d'image 
+ trou de filtrage. Nous utiliserons la distribution de probabilite p a pour simuler des repartitions. 
Nous notons (a^j = |a ?; , i G {l..N spot / p « e ;}} l'ensemble des amplitudes (non bruitees) des N spol/vixeI spot 
du pixel considers Ces amplitudes sont uniquement dues aux passages successifs de la lumiere dans le guide 
d'image, puis dans le trou de filtrage. De plus celles-ci sont par defaut des coefficient d'attenuation entre 
dans [0; 1], ou bien lorsque cela est. precise, les a 4 sont normalises dans [0; 200] afin de respecter la dynamique 
de 8 bits du systeme (cf. figures) . 



2.3.1 Mesures sur tissu 

Pour un ti ss u, nous - s upposons -qu o -nous avons j?,, 
>n iehse par li - >t 'h ut>1 i nrflw n \ 
compte de la variation d'ampli idea du system* les equations pra dentes d< iennent 



ations de la qual-i-te de couplago, 



, e (r p + a r -*-R mnx ^R{x^ r - 



avec l'indice i qui represente le ieme spot qui intersecte dans le pixel. Cette amplitude theorique est perturbee 
par le bruit de poisson pour produire une variable aleatoire poissonnienne Ni(R{x)). N etant grand, on peut 
approcher cette variable par la loi normal e : 



Ni(R(x)) - N , 



* N m 



Ur P + a, T * R max * R(x)), sfii * N guidK (r p + a 



*R(x))). 



Nous noterons avec Texposant (j) les realisations de cette variable Ny'(R(x)), j £ {1, N N ). Par comodite 
ces realisations sont notees simplement N-^ lorsqu'il n'y a pas d'ambiguite. 

Ces realisations sont ensuite normalises lineairement, et. quantifiees sur 8 bits. En pratique les varient, 
entre 0 et iV fl ,„ :rfe * (r p + a r * R max ). La normalisation doit prendre en compte deux phenomenes. D'abord, la 
distribution des coefficient a* module directement l'amplitude des : cela implique que toute la dynamique 
de ceux-ci est une dynamique utile. Ensuite, pour un o.i fixe la dynamique reelle a utiliser n'est pas [0;a» * 
Ngui,ie(r p + a r * R max )} mais [m * N guide r p ; o.i * N guide (r p + a r * R m ax)}- 

On ne peut cependant pas connaitre a priori quel est. la valeur du coefficient de couplage a,i pour une mesure 
nV\ et il faut done ne prendre en compte que le premier phenomene pour la correction de la. dynamique, 
et done y integrer la composante parasite qui sera retiree avec les procedures de calibration en amplitude 

Plusieurs methodes sont possibles : 

- methode des quartiles (ler quart jusqu'a 0.04, 2nd quart 0.125 , 3eme quart 0.305), 

- normalisation max/min, 

max = N guide * (r p + a r * R mnx ) 



= 0 



- normalisation moyenne (ou mediane) ; point fixe. Moyenne = 0.2, mediane = 0.125, point fixe=???. 

Idee du moment... : prendre un certain pourcentage a partir de 0 car les signaux les plus faibles sont aussi les 
plus nombreux et done les plus importa.nts pour une reconstruction des pixels. Les trois premiers quartiles 
representent 30.5% de la dynamique maximale si on ne prend pas en compte le coefficent de retrodiffusion 
R.(x), soit en nombre de photons : 0.305 *N guide *(r p + a T *R mnx ) = 15461. Si on prend 90% de la dynamique, 
on a alors 0.54 * N auide * (r p + o. r * R ma x) = 27373 photons a compter. 

Voici un tableau de correspondance entre le pourcentage de la dynamique et le pourcentage de mesure que 



% de mesures 


10 


20 


30 


40 


50 


60 


70 


80 


90 


95 


98 


99 


100 


% de la dynamique 


1 


3 


5.5 


8.5 


12.5 


18 


25.5 


36.5 


54 


69 


83.5 


90 


100 



En saturant a 54% de la dynamique les amplitudes (tres improbables) les plus grandes (done on en garde 
90% bien calibrees) la normalisation devient : 



n^(R(x)) = 



0.54 * N guide * (r p + a r * R max ) ' 



avec les nf'(R(x)) € [0; 1], les valeurs depassant les bornes etant alors saturees sur ces bornes. 
Avec cette normalisation, les n\ j \B.(x)) suivent la loi normale : 



n, Wx))~A{ 0 .54*(r p +a r *R max ) ' 0.54 * N guide * (r p + a r * R max ) > 
i, en saturant a la moitie de la dynamique maximum, on diminue et on homogeneise le bruit. 



Nous pouvons realiser les memes mesures sur un miroir; R. max = 1 et R(x) = 1.1 
re- G~r En -teha nt-compl-e-de-la- va riafri ori-d^m 



s notations sont similaires ; 



C i (R(x))=a i *N guide {r v + a r ), 



avec l'indice i qui represents le ieme spot qui intersecte dans le pixel. Cette amplitude theorique est perturbee 
par le bruit de poisson pour produire une variable aleatoire poissonnienne C*. C etant grand, on peut 
approcher cette variable par la loi normale 



Ci - N[ai * N guid , (r P + a r ), y/m * N guide (r p + a r )) . 



Nous noterons avec l'exposant (J) les realisations de cette variable C\ 3 \ j € {l,...,Nc}. Ces realisations 
sont normalises avec : 

C ' N guide (r p + a r y 
Nous avons done c\ j) qui est une realisation d'une variable aleatoire 



^( ai ^N gu Jr p+ ^)- 



2.3.3 Probleme 1 du moyennage 
Le probleme peut alors se formuler ainsi : 



Etant donnes un premier jeu de mesures acquises sur un 
miroir et un second jeu de mesures acquises sur un 

tissu (n^ {R(x))*) , quel est le meilleur estimateur R{x)de 



On peut, a priori distinguer trois methodes generiques : 

Moyennage temporel, puis sur le pixel : la calibration (c^) doit permettre d'estimer chaque a,;. Si 
on eonnait a,, alors on peut estimer plus farilement R.(x) a partir de (n\ j ^ (R(x))^ . 

Moyennage sur le pixel puis temporel : le pixel est reconstruit pour les donnees de calibration, et pour 
les donnees sur tissu, et ce a chaque instant de mesure. Puis la succession de pixels dans le temps est 
moyenne pour attenuer le bruit. 



Moyennage tempore] et sur le pixel simultanement : le moyennage utilise s'effectue simultannement 

sur le pixel et dans le temps. 
Un point commun des methode est, que Ton cherche a s'afranchir des ai dans l'estimation des N t afin de 
pouvoir retrouver R.{x). On cherche done a estimer : 

b . , n r v + a r *R max *R(x) 
= r p + a T *R mnx ■ 

ou R(x) est un estimateur biaise de R(x). La calibration ne doit, cependant pas s'arreter la car il reste un 
offset a soustraire, puis il faut retablir la dynamique. On realise done l'operation suivante pour oter le biais 
de R(x) : 

-3 Moyennage temporel, puis sur le pixel 



La somme du numerateur correspond a 1 'estimateur du maximum de vraissemblance qui permet. d'estimer 
le coefficient, a,;. On divise ensuite chaque terme de la somme du numerateur par l'estime de a,;. La somme 
du numerateur est alors elle aussi un estimateur du maximum de vraissemblance de R(x) |s 
Une fois que l'on a un estime de la valeur reelle de chaque spot, il faut les integrer afin de creer un pixel de 
l'image reconstruite. Nous proposons plusieurs methodes simples : 

R(x) = —^ (3) 

M spot/pixel t 

A(x) = Median 1 {A < (a)} (4) 
R(x) = OWAj {&(»)} (5) 
Avec OWA,; l'operateur de moyenne ponderee ordonnee. Soient (x,i) une suite de N valeurs reelles et (j/,:) la 
suite des valeurs (x.i) mais ordonnee par ordre selon une relation d'ordre predefinie. Alors : 

L'OWA permetde privilegier les valeurs qui sont plus grandes que les autres par rapport a la relation d'ordre. 
Dans notre cas precis la relation d'ordre est equivalents a la relation d'ordre "<" des amplitudes estimees a,; 
a partir des . 

4 Moyennage de spots pour reconstruire un pixel 



R u) (x) = - 



II faut remarquer ici qu'en 1'absence du bruit, de poisson, et. done en remplac.ant les variables n/' et. c/' 
leurs moyennes respoectives a,;i?((x)) et a* l'operation serait. : 
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Cet operateur est un operateur do moyenne ponderee en fonct.ion de l'importance des o.j. L'equation 6 
represents un operateur de moyenne pondere par l'estime Cette premiere operation pormet de tenir 
direct.ernent, compt.e de la fiabilite de ehaque mesure. cette fiabilite etant oxprimee par l'amplit.ude a,. Cet. 
amplitude regie le niveau du signal recu, et 1'importance du bruit poissonnien associe. L'equation 6 est, une 
moyenne ponderee ou les mesures les moins bruitees ont un poids plus important que les mesures bruitees. 
Nous disposons done de mesures du pixel a ehaque instant d'acquisition. II reste a les moyenner t.emporelle- 
ment. afin de reduire le reste du bruit. Nous proposons alors les operateurs vu precedement : 




5 Quantification 

| Tous ces cakuls sont. a verifier par un expert, du signal ? | 

5.1 Rapport signal sur bruit de la quantification 

On peut modeller la distribution des amplitudes des a* par une loi exponentielle decroissante. On a alors : 

avec a as jf 222 -, et. N ma x est le niveau maximum de la quantification (255 en 8 bits), le niveau min etant. 0 
par defaut. On divise l'intervalle [0; N ma x] en 2 N intervalles de longueurs identiques q — •^W*. Le rapport, 
signal sur bruit est donne en dB par : 

SW? dB = 101og 10 (^), 

avec P s la puissance du signal et P& la puissance du bruit. On a alors pour le bruit de quantification : 
b 12 12 2 2/v 

et, pour signal : 

P»=(J~Jae-°*dx)-(l~xae-™ dxf 

a 2 W 
_ _1_ 
~ a 2 

Le rapport signal sur bruit vaut done : 

SNR dB = lOio.9,0 ^ + 10/o 9l0 2 2/v = W - 1.25 

Soit. 46.75rIB sur 8 bits. 



5.2 Rapport signal sur bruit du au bruit de poisson 

Si on suppose que le bruit, de poisson est le bruit predominant, la puissance du signal non bruite vaut A 2 , A 
etant la valeur exacte du signal et done le parametre du bruit de poisson. La puissance du bruit de poisson 
vaut A (moment centre d'ordre 2). Le rapport signal sur bruit s'obtient done facilement avec : 

SNB. dB = 10/oflioA. 

Le bruit de poisson est un bruit de comptage du signal, et A represents bien ici le nombre de photon 
comptes dans le photomultiplicateur (PM). On ne prend pas en compt.e ici uniquement, le signal utile mais 
le signal entier car on ne peut pas differencier les deux avant l'etape de calibration informatique. II faut 
done tout, quantifier. Avec les parametres par defaut on compte environ 56000 photons dans le PM, et done 
SNB. dn = 47dB. 



— BiLan- 

On constate facilement que le niveau du bruit de poisson et le niveau de la quantification sont. equivalents. 
D d ax cho es Fu . soit i i ci it au pci un 1 el on suppo qui la c li l ration sen i arfaiti et. qui i >n 
am >as de pert oit on ,(1! quelqw bits pour etre ct iif'c i T i 1 l< 

Dans le premier cas (on a cru au pere noel mail il n'est. pas venu) il vaut. mieux saturer les niveaux hauts et 
done bien quantifier les niveaux bas car ils sont. plus nombreux. Dans le second cas, la question est. : combien 
faut-il rajouter de bits pour etre confortable? Ma reponse intuitive est qu'il faut 10 ou 12 bits au total. 
Derniere question : quel gain peut.-on esperer par bit. rajoute ? 



6 Resultats 

Les principaux resultats attendus sont : 

- comportement de la quantification sur 8 bits en fonction du contrfile du gain, 

- forme et comportement de la calibration (ce point est etroitement lie a la forme du moyennage), 

- comportement des moyennages en fonction de : 

- N N et C/v, 

- la fonction de moyennage (moyennage tempore! d'abord, ou pixel d'abord), 

- du nombre de spots dans un pixel, 

- du coefficient de retrodiffusion observe. 

Le nombre de parametres est assez important, mais nous pouvons tout de meme degager un certain nombre 
de resultats. 



6.1 Comparaison des deux types de methodes de moyennage 

On observe R(x) = 0.5. Les parametres utilises dans ce paragraphe : 

R(x) = 0.5 -4 e'est, le coefficient, normalise dans [0, 1], soit. une valeur reelle de 128 ici. 

Nc = 5 -4 nombre d'acquisitions pour la calibration 

A^^ = 5 -> nombre d'acquisitions pour la mesure sur tissu 

N spots /pixel = 3 ou 8 — > nombre moyen de spots pour un pixel dans le cas reel 
Statistiques sur 1000 mesures -> les statistiques sont. significatives 
La calibration ut.ilisee est. la calibration optimale par defaut. 

On simule la reconstruction du pixel avec les deux types de methodes. On repete l'experience 1000 fois afin 
de faire des statistiques. Les histogrammes des resultats obtenus sont represent.es sur les figures 3 et 4. Le 
tableau suivant donne les moyennes, les medianes et. 1'ecart. typde des mesures obtenues pour chaque methode 

{N spots/pixel = 8). 



Met.hode 


Moyenne 


Mediane 


Ecart-type 


Moy. temp, puis moy. pixel 


130.5 


128 


7.2 


Moy. temp, puis median pixel 


127.4 


127.5 


1.4 


Moy. temp, puis OWA pixel 


128.2 


127.7 


2.2 


Moy. pixel puis moy. temp. 


127.6 


127.6 


0.9 


Moy. pixel puis median temp. 


127.6 


128.0 


1.0 


Moy. pixel puis OWA temp. 


127.8 


127.9 


0.9 



On remarque que les moyennes temporelles sont bien moins performantes que ^integration dans un pixel. 
Cela peut s'expliquer par le fait que l'integration dans un pixel est implicitement realise par une moyenne 
ponderee par rimportance intrinscque du bruit, poissonien de la mesure. Le moyennage temporel s'effectue 
lui sur des donnees forcement mal calibrees et sn revele moins bon. 



Le tableau suivant donne les moyennes, les medianos et. l'ecart typde des mesures obtenues pour chaque 
met.hode avec N spoU / pixel = 3 (les figures correspondanl.es sont les figures 5 et. 6). 





Methode 
-Meyr-temp^ puis moy. pixel - 
Moy. temp, puis median pixel 


Moyenne 
128.2 


M . ; i , ■ 
128.0 


Ecart-type 
8.1 






Moy. .temp, puis OWA pixel 
Moy. pixel puis moy. temp. 
Moy. pixel puis median temp. 
Moy. pixel puis OWA temp. 


129.3 
— 127.8 
127.8 
128.3 


— 127.7 
127.6 
128.0 
128.0 


7^2 

2.5 " 

3.3 

2.7 





On remarque ici que les resultat enter les deux families de methodes sont plus proches. L'integration dans le 
pixel semble encore meilleure, au moins pour la moyenne. Ce test ne change pas a conclusion, mais montre 
bien la sensibilite de la reconstruction par rapport au facteur de sur-echantillonnage. 



6.2 Estimation des parametre "optimaux" 

A FAIRE : 

- determiner le nombre d'images necessaires a la calibration, 

- determiner le nomnre d'image necessaire pour les mesures sur tissu, 

- quantifier I'influence du decalage lors de la calibration, 

- etude theorique analytique de la calibration, 

- influence de la repartition des amplitudes des spots, 



7 Conclusion sur le moyennage 

La premiere conclusion qui s'impose est le type de moyennage a effectuer : les donnees de calibrations doivent 
etre globalement. moyennees avec : 

No N, pMlpi „, 

puis, utiliser cette "calibration" pour integrer les mesures de chaque acquisition dans le pixel. Le moyennage 
temporel est effectue en dernier en moyennant les pixels des images consecutives. Ce dernier type de moyen- 
nage importe pen car les resultats obtenus dans la partie 6.1 son similaires pour la moyenne arithmetique, 
la medianne, et l'OWA. Le plus simple est done d'utiliser la moyenne arithmetique qui represente le cout de 
calcul le plus faible. 
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histogramme des mesures calibrees reconstruites dans un pixel 
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Fig. 3 - Histogramme des mesures calibrees et. reconstruites dans un pixel (N spots /p ixe i = 8) pour les 
parametres de la partie 6.1. On effectue d'abord le moyennage temporel, puis Integration dans le pixel a.vec 
les trois methodes differentes. 



11 



histogramme des mesures calibrees reconstruites dans un pixel 
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Fig. 4 - Histogramme des mesures calibrees et reconstruites dans un pixel {N svois / vixe i = 8) pour les 
parametres de la partie 6.1. On effectue d'abord Integration dans le pixel, puis le moyennage temporel du 
pixel avec les trois methodes differentes. 
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histogramme des mesures calibrees reconstruites dans un pixel 
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FlG. 5 - Histogramme des mesures calibrees et. reconstruites dans un pixel pour les parametres de la partie 
6.1 (N spots / pixe j = 3). On effectue d'abord le moyennage tempore!, puis 1'integration dans le pixel aver, les 
trois methodes differentes. 
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histogramme des mesures calibrees reconstruites dans un pixel 
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Fig. 6 - Histogramme des mesures calibrees et. reconstruites dans un pixel pour les parametres de la partie 
6.1 {N apoU / V i X ei = 3). On effectue d'abord l'integration dans le pixel, puis le moyennage tempore! du pixel 
avec les trois methodes differentes. 



14 



histogramme des mesures calibrees reconstruites dans un pixel 
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Fig. 7 - Histogramme des mesures calibrees et. reconstruites dans un pixel pour les para.met.res de la partie 
6.1 (N spols / pixel = 8). On effectue d'abord le moyennage temporal, puis Integration dans le pixel avec les 
trois methodes differentes. Les mesures sont, sur-normalisees (=t,rop faibles). 
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histogramme des mesures calibrees reconstruites dans un pixel 
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FlG. 8 - Histogramme des mesures calibrees et reconstruites dans un pixel pour les parametres de la partie 
6-1 {N spots / pixel = 8)- On effectue d'abord Integration dans le pixel, puis le moyennage tempore] du pixel 
avec les trois methodes differentes. Les mesures sont sur-normalisees (=trop faibles). 
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histogramme des mesures calibrees reconstruites dans un pixel 
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Fig. 9 - Histogramme des mesures calibrees et. reconstruites dans un pixel pour les parametres de la pa.rt.ie 
6.1 (N apotll - 8). On effectue d'abord le moyennage temporal, puis 1 'integration dans le pixel avec les 
trois methodes differentes. Les mesures sont sous-normalisees (=sat.urees). 
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histogramme des mesures calibrees reconstruites dans un pixel 
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Fig. 10 - Histogramme des mesures calibrees et reconstruites dans un pixel pour les parametres de la partie 
6-1 (N S p 0 t s /pixel = 8). On effectue d'abord l'integration dans le pixel, puis le moyennage tempore] du pixel 
avec les trois methodes difFerent.es. Les mesures sont, sous-normalisees (=saturees). 
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8 Restauration sans calibration 

Cette partie propose une methode de restauration des images de mesures n'utilisant pas le calibrage en 
amplitude. Le but est de s'affranchir des phenomene d'echantillonnage sur le guide d'images. Les deux 
principaux phenomenes sont : 
le probleme d'injection non uniforme : 

aliasing : du a une frequence d'echantillonnage inferieure a la frequence de Nyquist. 
8.1 Analyse dans le plan de fourrier 

Supposons que le spot balaye continument une seule fibre, sur toute son diametre. Mors la fonction d'injection 
est (pour un spot du meme diametre que la fibre) : 



en supposant la fibre centree en x = 0, alors x indique la position du spot par rapport a la fibre, et rf ihrf: 
TF qui est aussi une gaussienne : 

TF( ai )(u) =r fibre ^ er 





Fig. 11 - (a) est la fonction d'injection dans une fibre en fonction de la distance entre le centre du spot laser 
et le centre de la fibre, (b) est la TF de (a). 

Soit v c la frequence de coupure approchee. La portion de gaussiene repliee (aliasing) a la forme d'une queue 
de gaussienne c,rfc{^ On pent done approcher une frequence de coupure a v c — 3a„ — y^,,-^, = 
5.40 10 5 . On en deduit la frequence d'echantillonnage optimale (spatiale) a ^^.^ . ce qui donne une periode 
spatiale «optimale» de 0.93/im pour r Iibre = 1.25//. m. Si on prend v € = 2r7„, alors la periode spatiale est de 
1.38//, m. 

La fibre etant circulaire, on peut supposer que cette frequence spatiale doit etre respect.ee en x et en y. 
En pratique la periode temporelle d'echantillonnage est 120n.s. Si on suppose une vitesse constante du spot 
d'environ llm.s" 1 , cela donne une periode spatiale de 1.32/im en x. En y la. periode spatiale est variable en 
fonction de la vitesse du miroir galva. Les parametres par defaut donnent une frequence de raffraichissement 
de 15 images par seconde, avec une periode ligne de 125/z.s. En supposant un balyage sur 800/xm la periode 
spatiale en y est de 1.5//.m. 
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Ces chiffres indiquent que Ton va faire apparaitre de l'alisasing sur les images, avec une frequence en x 
differente de la frequence en y. L 'aliasing sera plus important suivant la direction y. 

La figure 12-(a) represents la TF de l'image 12(b) qui est une image simulee, avec le guide d'image tourne 
de 10 degres, qui image une surface uniforme, et avec un echa.ntillonna.ge instable a 1.5//. m pres pour chaque 
mesure. 




(a) (b) 
Fig. 12 - (a) Image, (b) FFT de (a). 

La forme centrale est due a la forme du guide. Les 6 pics les plus importants sont le signal utile melange aux 
pics representant le maillage hexagonal des fibres. Les pics secondaires ressemblant aux precedents sont de 
l'aliasing. Le reste du plan de fourrier ressemble a du bruit blanc. 

On remarque que le pic representant la frequence quasi verticale des fibres est tres proche de 1, ce qui 
indique que Ton echantillonne avec une periode proche de celle des fibres. II faut done augmenter la frequence 
d'echantillonnage verticale. Si on double cette frequence, on diminue tres nettement l'aliasing (figure 13). 
Celui-ci reste cependant nettement visible sur le TF, mais n'apparalt pas sur l'image. Ceci permet de justifier 
les calculs precedents... 

8.2 Technique de masquage de frequence 

8.3 Resultats 

9 Modele parametrique de restauration 

10 Stabilisation d'image 

i.e. recalage simple. 
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Fig. 13 - (a) Image, (b) FFT de (a). 

10.1 Recalage ID pour translation 

10.2 Calibration des projections 

10.3 Correlation croisee normalisee 

10.4 Resultats 
11 Conclusion 
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